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In  this  study,  we  deal  with  the  analysis  and  synthesis  of 
two  classes  of  multivariable  feedback  systems,  high  gain  feedback 
systems  and  variable  structure  systems  (VSS),  subject  to  parameter 
variations  and  disturbances. 

We  first  examine  the  insensitivity  property  of  high  gain 
state  feedback  systems.  The  idea  of  incorporating  high  gain  feedback 
loops  in  observers  is  then  introduced.  We  then  study  the  behavior  of 
the  resulting  "two-high-gain- loops"  system.  A design  procedure  which 
enhances  the  insensitivity  of  the  closed  loop  system  for  the  high  gain 
feedback  matrices  is  developed.  Second,  we  establish  the  relationship 
between  high  gain  feedback  systems  and  VSS.  We  reveal  the  insensitivity 
property  of  VSS.  Then,  observers  with  variable  structure  are  developed 
and  we  examine  the  behavior  of  the  resulting  VSS  with  variable  structure 
observers.  Third,  we  investigate  the  robustness  property  of  these  two 
classes  of  feedback  systems  with  respect  to  model  reduction.  Neglection 
of  acuator  and  sensor  dynamics  is  considered.  Finally,  we  exploit  the 
insensitivity  of  these  two  classes  of  feedback  control  to  enhance  the 
reliability  of  interconnected  systems.  Design  procedures  of  high  gain 
feedback  and  variable  structure  controllers  are  illustrated  in  the  designs 
for  a distillation  column,  a synchronous  machine  and  the  longitudinal 
motions  of  an  aircraft. 
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CHAPTER  1 
INTRODUCTION 

High  gain  feedback  has  been  a classical  tool  for  reduction 
of  effects  of  disturbances,  parameter  variations  and  distortions.  Al- 
though limited  to  single  input-single  output  feedback  systems,  the 
early  investigations  of  structures  permitting  high  gains  [1],  the  rules 
for  root  locus  asymtotes  [2]  and  the  results  on  sensitivity  and  return 
difference  [3,4]  has  greatly  deepened  the  intuition  of  control  engineers 
in  the  1950's.  Recent  developments  in  the  multivariable  system  theory 
have  revived  the  interest  in  high  gain  feedback  systems.  First,  recent 
studies  on  multivariable  transmission  zeros  [5,6]  and  the  work  on  multi- 
variable  root  locus  asymptotes  [7]  generalize  the  notion  of  zeros  of 
transfer  functions  in  single  input-single  output  systems  and  consider  the 
asymptotic  behavior  of  multivariable  high  gain  feedback  systems  when  gain 
tends  to  infinity.  Second,  works  on  disturbance  rejection  [S],  parameter 
uncertainty  [9]  and  decoupling  of  large  scale  systems  [10]  either  pur- 
posely introduce  high  gain  feedback  in  the  problem  statement  or  they 
implicitly  appear  in  the  resulting  feedback  structures.  Third,  feedback 
implementations  of  linear  optimal  controls  when  only  small  penalties  are 
made  on  the  control  variables  (the  so-called  "cheap"  control  problem) 
result  in  loops  with  high  gain  [11,12,13]. 


Another  class  of  feedback  controls  capable  of  reducing  parameter 
sensitivities  and  rejecting  disturbances  is  the  so-called  variable  struc- 
ture control.  Basically,  it  is  a feedback  control  discontinuous  on  some 


switching  surface  defined  on  the  state  space.  This  cl  tss  of  feedback 
svsteos  is  called  variable  structure  systems  (VSS’'.  It  has  been  developed 
in  the  USSR  in  the  last  fifteen  years  [ 14 , 1 5 , 1 b , 1 7, 18 ] . For  a recent 
survey,  see  [19].  A salient  feature  of  VSS  is  the  so-called  sliding  mode 
occurs  on  the  switching  surface.  It  is  in  sliding  mode  that  VSS  becomes 
s. *' s s sensitive  to  parameter  variations  and  disturbances. 

In  this  study,  uv  deal  with  the  analysis  and  synthesis  of  these 
two  classes  of  mil t ivar table  feedback  systems,  high  gain  feedback  systems 
and  VSS,  subject  to  parameter  variations  and  disturbances . First,  we 
examine  the  insensitivity  property  of  high  gain  state  feedback  systems. 
Then,  in  the  case  when  only  some  of  the  states  variables  are  accessible, 
we  introduce  the  idea  of  incorporat i ng  high  gain  feedback  loops  In  ob 
servers.  By  allowing  high  gain  feedback  in  the  observer  structure,  the 
observation  error  dynamics  enjoy  the  same  insensitivity  property  inherent 
in  high  gain  feedback  systems.  In  place  of  high  gain  state  feedback,  the 
"high  gain"  observer  states  are  fed  back  through  the  main  high  gain  feed- 


back loops,  '.v'o  study  the  behavior  of  this  "t vo-htgh-ga in- loops"  feed- 
back system.  A practical  alternative  to  examining  the  sensitivity  of  all 
the  states  of  high  gain  feedback  systems  is  to  investigate  only  the 
sensitivity  o*’  the  variables  which  are  critical  to  performance  degradation 
or  have  to  meet  certain  design  specif icat ions . Wo  call  these  variables, 
the  regulated  variables.  We  discover  that  there  remain  some  degrees  of 
freedom  in  the  design  of  the  high  gain  feedback  matrices  which  can  be  util- 
ized to  enhance  the  insensitivity  of  high  gain  feedback  systems.  The  num- 
ber of  degrees  of  freedom  depends  on  the  number  of  control  variables, 
available  measurements , disturbance  Inputs  and  regulated  variables.  \ 


J 


] 


procedure  Is  developed  for  the  design  of  the  feedback  matrices  which 


exercises  the  available  degrees  of  freedom. 


Second,  we  establish  the  relationship  between  high  gain  feed- 


back systems  and  VSS.  We  find  that  in  sliding  mode,  VSS  enjoy  the  same 


insensitivity  property  in  high  gain  feedback  systems.  This  encourages 


us  to  develop  observers  with  variable  structure  feedback  in  parallel  to 


high  gain  observers.  We  then  proceed  to  examine  the  behavior  of  variable 


structure  feedback  svstems  with  variable  structure  observers.  The  same 


degrees  of  freedom  in  the  design  of  variable  structure  feedback  is 


discovered  as  in  high  gain  systems. 


Third,  we  address  an  important  consideration  in  application: 


the  robustness  property  of  these  classes  of  feedback  systems  with  respect 


to  model  reduction.  As  a step  in  this  direction,  we  consider  the  most 


common  model  reduction  in  practice,  the  neglection  of  actuator  and  sensor 


dynamics.  This  robustness  property  differentiates  VSS  and  high  gain 


feedback  systems.  As  a reward  for  the  added  complexity,  VSS  are  more 


robust  with  respect  to  the  neglected  small  time  constants.  This  is 


due  to  the  fact  that  variable  structure  control  does  not  force  the  motion 


to  be  fast. 


Finally,  we  exploit  the  insensitivity  property  in  high  gain 


feedback  svstems  and  VSS  to  enhance  the  reliability  of  interconnected 


systems  in  the  spirit  of  [20].  The  design  procedures  of  high  gain  feed- 


back and  variable  structure  controllers  are  illustrated  by  the  design 


of  control  systems  for  a distillation  column,  a synchronous  machine  and 


the  longitudinal  motions  of  an  aircraft. 
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CHAPTER  2 

SINGULAR  PERTURBATION  AND  HIGH  GAIN  FEEDBACK  SYSTEMS 

The  high  gain  linear  time-invariant  feedback  system  we  consider 
is  of  the  form 


x = A x + B u 
o o 


(2.1) 

(2.2) 

where  g is  a large  positive  scalar,  the  state  x is  an  n-vector  and  control  u 


u = gC  x 
o 


is  an  m-vector.  We  assume  that  the  pair  (A  ,B  ) is  controllable, 

o o 

matrices  A , B and  C are  known  and  rank  B = rank  C = m.  In  [21],  it 
o o o o o 

is  observed  that  a dynamic  loop  whose  gain  tends  to  infinity  causes  change 
in  the  system  order,  characteristic  of  singular  perturbation.  On  the 
other  hand,  a typical  singularly  perturbed  linear  time -invariant  feedback 
s vs  tern  is 


X1  = (F11  + FlVXl  + (F12  + GiH2)x2 

Ux2  = (F21  + G?H1)x1  + (F,2  + G9H2)x0 
where  y is  a small  positive  scalar.  Clearly  by  defining 


u = -[F21  + «2l'l)xl  + F22  + G2H2)x2] 


(2.3) 

(2.4) 


(2.5) 


to  be  an  artificial  high  gain  feedback,  the  singularly  perturbed  system 
(2.3),  (2.4)  is  equivalent  to  a high  gain  feedback  system  (2.3),  (2.5)  and 

x9  = u (2.6) 

where  the  gain  factor  g = --  and  g * ” as  )J  •*  0.  Motivated  by  this  analogy, 
we  use  the  singular  perturbation  methodology  to  examine  the  behavior  of 


5 

of  high  gain  feedback  systems.  In  particular,  we  focus  on  the  class  of 
high  gain  feedback  systems  which  possesses  a two  time  scale  property. 

In  the  case  when  all  the  states  are  accessible,  it  is  possible  to  design 
the  feedback  matrix  K in  the  high  gain  feedback  loop 

u - gKx  . (2.  7) 

A procedure  for  separate  placement  of  fast  and  slow  modes  and  a de- 
composition of  near  optimal  high  gain  regulator  problems  are  presented. 

2.1.  Fast  and  Slow  Modes 

The  high  gain  feedback  system  (2.1)- (2. 2)  with  y H g 1 can  be 

written  as 

UX  - ( UA  + B C )x  . (2.8) 

o o o 

In  order  to  apply  singular  perturbation  methods,  we  transform  (S'!  into 
the  form  of  system  (2.3),  (2. A).  Let  M be  such  that 


3 ~ 

(2.9) 

o 

2_ 

where  is  a m x m nonsingular  matrix.  Such  a matrix  M exists  and  is  the 

T 

product  of  elementary  row  transformations.  Let  x*  - [x|  x^]  = Mx  and  in 

the  new  coordinates,  (2.8)  becomes 


A x * + A v * 

n 1 12  2 


Ux2  “ ;iA21  +b:Ci'Ix[  + <>A22  + B2CVX2 


(2.10) 

(2.11) 


where  C^,  C.,  are  m x (n-m)  and  m x m matrix  respectively  defined  by 


and  matrices  A^j  are  defined  in  Appendix  A.  System  (2.10),  I 1 1 is  new 

in  the  same  form  as  (2.3),  (2.4).  Nevertheless,  since  it  is  required  In 

singular  perturbation  analysis  r22,23]  that  t lie  matrix  F.,.,  t C0H.,  in  (2.4) 

be  a stable  matrix,  we  nuke  the  assumption  that  u A°,,  4 b!^C,  is  a stable 

matrix.  For  u sufficiently  small,  that  is,  the  gain  v sufficiently  large, 

this  assumption  is  equivalent  to  Re\(fC’c.,)<  0.  Since  is  nonsingular, 

this  requires  C.,  to  be  nonsingular.  From  the  identity  C B “ C 11°,  this 
- o o 2 2 

assumption  becomes 


Re  ).(<’  B ) < 0 . 
v o o 


(2.1 3) 


The  following  lemma  expresses  assumption  (-’.1  11  in  terms  of  the  matrix 


B C in  (2.8) 
o o 


Lemma  2.1: 


Assumption  (2.13)  Is  equivalent  to 


rank  B C = rank  C 11  “ m 

o o o o 


(2.14) 


and  the  nonzero  eigenvalues  of  11  C have  negative  real  parts, 

o o 1 


Re  \ (B  C.o)  • 0,  i“l m . 


(2.13) 


Proof:  Let  \ be  an  eigenvalue  of  C 11  . then  there  exists  a vector  v 
o o 

such  that 


C B v “ \v  . 
o o 


(2.16) 


Multiplying  to  the  left  with  B on  both  sides  of  (2.16),  we  have 

B C w ■ \w 
o o 

where  w s B v.  Thus  \ is  a nonzero  eigenvalue  of  B 0 . Thus  (2.1')  is 
O 0 o v 

equivalent  to  (2.14)  and  (2.15). 


- — 


pi  p p 


'"vr~'~Tr 


Under  the  assumptions  (2.14),  (2.15),  it  will  now  be  shown  that 
the  motion  of  (2.8)  consists  of  a fast  transient  to  an  OCy)  neighborhood 


of  the  subspace  Cqx  = 0,  followed  by  a slow  motion  in  this  neighborhood. 


Let  x ■ [x  x2!  = Tx  where 


n-m 


LC1 


-2J 


M a Tm 


(2.17) 


We  note  that  x^  = x|  and  x2  ■ Cqx.  It  is  shown  in  Appendix  A that  the 
meaning  of  T is  that  it  decompress  the  state  space  of  x into  the  null 


space  of  Cq,  V (Cq),  and  the  range  of  Bq,  R (Bq).  Then  x can  be  written  as 


1 ~ 


x = N x,  + B (C  B ) x0 
1 o o o 2 


(2.18) 


where  the  columns  of  the  n * (n  - r)  matrix  N span  the  null  space  of  Cq, 


that  is. 


C N * 0 
o 


(2.19) 


and 


x^  = ML^x 

where  matrix  is  the  first  n-m  rows  of  M,  that  is 


(2.20) 


M,B  = 0 . 
1 o 


(2.21) 


Such  decompositions  also  appear  in  [6,7].  The  system  (2.8)  is  thus 
transformed  into 


X1  = Fllxl  + F12X2 


(2.22) 


= pH1x1  + (p^  + CqBo)x2 


(2.23) 


where  H^,  are  defined  in  Appendix  A.  The  block  triangularization 


[24,25]  which  is  simpler  than  Jordan  transformation  used  in  [6,7] 
is  now  applied  to  system  (2.22),  (2.23)  to  exhibit  the  two  time 


8 


scale  property  of  high  gain  feedback  system  (2.8).  We  introduce  the  "fast" 
variable 

= x9  + UI.x1  (2.24) 

where 

L = Lq  + UG  (2.25) 

Lo  2 (CoBo  +MH2)’1h!  <2*26) 

and  a recursive  formula  for  calculating  G is  given  in  Appendix  B.  The 
resulting  separation  of  fast  and  slow  modes  is  now  summarized  in  the 
following  theorem. 


Theorem  2.1: 


Let  the  upper  bound  for  P (that  is  a lower  bound  g^  for  g) 
be  given  in  (B3).  If  the  high  gain  feedback  system  (2.8)  satisfies 
assumptions  (2.14)  and  (2.15)  and  if 


0 * P £ p 


(2.27) 


then  (2.8)  is  equivalent  to  the  block  triangular  system 


*1  ' LFll-^i2Ljxl  +F12’ 


(2.28) 


^ = "CoBo+UH2  + U LF12r 


(2.29) 


and  hence  its  eigenvalues  are  approximately. 


Tf  -L'W  + °o>'!  » 


(2. 30) 


= Xj(Fll>  + °^>‘ 


j=l , • • • , n-m  . 


(2.31) 


More,  for  sufficiently  high  gain,  that  is,  for  p sufficiently  small,  the 
fast  subsystem  (2.29)  is  asymtoticallv  stable. 


9 


Proof:  The  proof  follows  from  (2 . 17) - (2 .26)  and  the  existence  of  G is 
established  in  Appendix  B.  The  asymototic  stability  of  (2.29)  follows  from 
Assumption  (2.13). 


The  two  time  scale  property  of  high  gain  feedback  system  (2.8) 

is  exemplified  by  the  presence  of  O(-jj)  large  eigenvalues  and  the 

eigenvalues  \ which  are  0(1).  The  fast  variable  T)  decavs_exponentially  in 

-8Us 

the  "stretched"  time  scale  t/U  and  for  t 2 tg,  it  is  0(e  ) where 

X ■ | max  Re  X(CqBo)| . Applying  lemma  Cl,  the  slow  variable  can  be 


approximated  by 


- F11C~ 

xx(t)  - e Xl(0)  + 0(u)  . 


(2.32) 


Hence,  the  decomposition  (2.18)  of  the  original  state  space  corresponds 
to  the  separation  of  time  scales  and  x(t)  is  approximated  by 

Ft  C B ~ 

x(t)  = Ne  11  xx(0)  + Bo(CQBo)'1e  0 ° UT)  (0)  + 0(u)  (2.33) 

which  proves  our  earlier  assertion  that,  in  the  limit,  the  fast  transient 

occurs  in  R(B  1 and  the  slow  motion  is  confined  to  N(C  ).  For  t 2 t . 

o o s 


x2(t>  = C^x(  t)  - liLox1(t)  + 0(p  ) . 


(2.34) 


If  z = Cqx  is  defined  as  the  regulated  variables  of  system  (1>  to  be 

forced  to  zero,  that  is,  a "zeroing  the  output"  problem  [26]  is  considered 

where  z is  the  output,  then  from  (2.34),  we  see  that  z(t)  is  reduced  to  an 

2 

0(u  ) quantity  after  t 2 t . 

We  have  shown  that  under  the  assumption  (2.141  and  (2.151,  the 
high  gain  feedback  system  (2.81  possesses  a two  time  scale  property. 

The  recent  study  of  multivariable  root  locus  asymptotes  [7  ] indicates  that 
for  more  general  structures  of  Bq  and  C , instability  of  the  closed  loop 


LO 


system  can  result  from  high  gain  feedback.  For  example,  in  single  input 

system,  (2.14)  means  that  the  inner  product  of  the  input  and  output 

vector c be  nonzero.  In  [27],  singular  perturbation  methods  are  used 

to  analyze  the  behavior  of  the  single  high  gain  loop  system.  It  is 
T 

shown  that  when  c ^b^  = 0,  the  closed  loop  system  contains  high  frequency 

T i 

oscillation  modes.  It,  in  addition,  c A b = 0 for  some  j > 0,  then  the 

’ o o o J 

system  becomes  unstable.  Thus,  we  focus  on  high  gain  feedback  systems 
which  satisfy  assumptions  (2.14)  and  (2.15).  Clearly,  if  all  the  states 
are  accessible,  then  there  always  exists  feedback  gain  matrices  such 
that  (2.14)  and  (2.15)  holds. 

2.2.  Root  l.ocus  Asymptotes  and  Transmission  Zeros 

In  Theorem  2.1,  we  have  revealed  the  asymptotic  behavior  of  the 
eigenvalues  of  the  high  gain  feedback  system  (2.8)  as  gain  tends  to 
infinity.  For  single  input  system,  the  asymptotic  behavior  of  the  eigen- 
values can  be  determined  by  the  well  known  root  locus  method  [5]  which 
provides  graphical  estimates  of  the  closed  loop  eigenvalues  as  a function 
of  the  loop  gain  if  the  open  loop  transfer  function  is  known.  The 
finite  eigenvalues  of  the  closed  loop  system  coincide  with  the  zeros  of 
the  open  loop  transfer  function  while  the  large  eigenvalues  approach 
the  root  locus  asymptotes  in  the  limit  as  gain  tends  to  infinity.  In 
[7],  some  results  of  the  multivariable  root  locus  asymptotes  have  been 
obtained.  The  relationship  of  the  eigenvalues  of  a high  gain  feedback 
system  with  the  zeros  of  the  open  loop  system  is  established  in  [5  ]. 
These  zeros  are  specifically  termed  transmission  zeros  since  in  general 
they  are  quite  different  from  the  zeros  of  the  various  transfer  functions 
relating  the  inputs  and  outputs  of  the  system.  Later,  we  will  interpret 


I 


(2.30)  and  (2.31)  in  these  terras.  First,  we  define  transmission  zeros. 


Definition  1:  [ 5 ] 


Given  the  system  (2.1)  with  input  u and  output 


y - C x , 
o 


(2.35) 


where  y is  an  r-vector,  the  transmission  zeros  of  (2.1),  (2.35)  (or  in 

short  the  triple  (C  ,A  ,B  ))  are  defined  to  be  the  set  of  complex  number  C 
o o o r 

which  satisfy  the  following  inequality 


A - 01  B 
o o 


< n + min(r,m)  . 


(2.36) 


This  is  one  of  the  many  definitions  of  transmission  zeros  given  in  [ 5 ] . 
An  equivalent  definition  for  the  case  when  m*r  is  given  as  follows. 


Definition  2:  [5,48] 

Suppose  the  pair  (Ao>Bq)  is  controllable,  the  pair  (A^,Co)  is 

obserable  and  m=r,  then  the  transmission  zeros  of  the  triplet  (C  ,A  ,B  ) 

o o o 

are  the  zeros  of  the  polyominal 


det(sI-A  )det[ C (sI-A  ) B ] - 0 . 
o o o o 


(2.37) 


We  note  that  if  m = r-  1,  then  the  zeros  of  (2.37)  are  in  fact  the  zeros 
of  the  open  loop  transfer  function. 

From  (2.30),  it  can  be  seen  that  as  g-*“,  u -*  0,  the  eigenvalues 
tend  to  infinity  along  the  asymptotes  defined  by  the  directions  of 
X^(C^B^)  which  are  the  root  locus  asymptotes  obtained  in  " 7j.  For  large 
but  finite  values  of  gain  g,  the  eigenvalues  can  be  computed  from 
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neglecting  the  0(y“)  terms  in  2.29). 

To  connect  the  slow  eigenvalues  in  (2.31)  with  transmission 
zeros,  we  need  the  next  theorem  proved  in  [ 5 ]. 


Theorem 


Let  m = r and  high  gain  output  feedback  u = gy  be  applied  to 


(2.1),  then  the  finite  eigenvalues  of  the  matrix  A + gB  C coincide  with 

o o o 

the  transmission  zeros  of  the  tripie  (C  ,A  , B ). 

o o o 

Thus,  the  limits  'j(F^)  of  the  eigenvalues  of  (2.31)  are  the  transmis- 
sion zeros  of  the  triple  (C^,A  ,B  ).  For  large  but  finite  values  of  gain  g, 

O 

the  eigenvalues  'j  can  he  approximated  to  any  order  of  \i  from  the  matrix 

•y 

fjj  - U F^0G.  Hence,  we  conclude  that  for  the  high  gain  feedback  system 
(2.1),  (2.2)  which  sat isf ies  condi tions  (2.14)  and  (2.15),  Theorem  2.1  encom- 
passes and  extends  the  results  on  root  locus  asymptotes  and  transmission  zeros 
of  multivariable  systems  in  [5,6,7  ].  Moreover,  Theorem  2.1  and 
Appendix  A provide  the  following  procedure  for  calculation  of  transmission 


Find  the  matrix  Ml  = [ M*  in  (2.9)  and  its  inverse  M 


[Sj  S.,  1 where  M0  and  S.,  are  n < m matrices.  This  inverse  is  easy  to  com- 
pute since  M is  a product  of  elementary  row  transformations.  Then  form 
C^M  * = [Cj.  C0  ] as  in  (2.12).  Compute  and  the  eigenvalues  of 


Fn  = M1Ao(S1  - S2C^  ACX) 


(2.38) 


which  are  the  transmission  zeros  of  the  triple  (C  ,A  ,B  ). 

o o o 

The  derivation  of  (2.38)  is  given  in  Appendix  A where  it  is  also 
shown  that  the  columns  of  the  matrix  - S,C.,  ^ span  the  null  space  of 


C , Chat  is, 
o 

R (N)  - R (St  - S2C;1C1)  (2.39) 

where  N is  defined  in  (2.19).  Thus,  our  matrix  coincides  with  the 

"zero  matrix"  given  in  r6  , 7,9]  for  the  specific  choice  of  the  basis  for 

)J(C  ).  The  calculation  of  the  transmission  zeros  from  matrices  of  the 
o 

form  (2.38)  avoid  the  ill-conditioning  due  to  the  presence  of  the  largo 
gain  factor  g in  the  direct  computation  of  the  eigenvalues  of  the  matrix 
A i + g which  is  the  approach  used  in  [ 5,23J- 

2.3.  Separate  Pole  Placement 

In  the  case  when  all  the  states  are  accessible,  we  consider  the 
synthesis  of  high  gain  state  feedback  control  u ° gKx  (2.7)  such  that  the 
closed  loop  system  (2.1),  (2.7)  possesses  desirable  properties.  In  this 
section,  we  are  concerned  with  the  synthesis  of  matrix  K such  that 

rank  KB  " rank  B K " m (2.40) 

o o 

and  the  fast  and  slow  eigenvalues  of  the  high  gain  state  feedback  system 
(2.1),  (2.7)  are  placed  at  desirable  locations.  By  substituting  K for  C} 
in  (2.2),  Theorem  2.1  shows  that  the  fast  and  slow  eigenvalues  are  approxi- 
mated by 

j~  C\ i (K  Bo)  + 0(u)]  i-1,  ....  m (2.41) 


and 
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det[K  (q  I -A  )'  B ] = 0, 
L z'^j  o o ’ 


j*l n-m  . 


(2.47) 


or  alternatively  the  linear  equations  [ 3 l] 


det 


A - q . I B 
o n j o 


0,  j«l , . . . , n-m 


(2.48) 


are  satisfied.  We  note  that  (2.48)  defines  q^  as  the  transmission  zeros 
of  (Kt,Ao>B^)  as  in  (2.36).  The  placement  of  the  n-m  transmission  zeros 
implies  that  K^R  is  nonsingular  [ 6].  Thus,  K^  satisfies  (2.40)  and  the 


-1 


fast  eigenvalues  in  the  limit  are  \.(K  B ).  Since  (KB)  exists,  there 

i z o z o 


exists  matrix  K.  such  that 
i 


^i^KiKzBo^  =Pi’  1=1  ’ •••’  m * 


(2.49) 


The  transmission  zeros  of  (K  ,A  ,B  ) and  (K.K  ,A  ,B  ) are  the  same  since 

Z O O x zoo 


rank 


’ A -crl  B 

A -cl  B 

o o 

o 0 

° rank 

K K 0 

K 0 

i z 

z 

(2.50) 


This  proves  the  lemma. 

A procedure  to  place  the  slow  and  fast  eigenvalues  is  developed 
in  the  proof  of  this  lemma.  The  design  consists  of  two  stages,  first  the 
slow  eigenvalues  are  placed  by  solving  for  K^  in  (2.47)  or  (2.48).  Then, 
the  fast  eigenvalues  are  placed  by  solving  for  K^in  (2.49).  This  design 
procedure  encompasses  the  design  in  [29]  where  it  is  demanded  that  K3q 
be  negative  definite  which  is  equivalent  to  restricting  p.  in  (2.44)  to  be 


real  numbers. 
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We  note  chat  the  placement  of  the  slow  eigenvalues  in  this 

design  is  achieved  by  finding  an  appropriate  "output"  matrix  K so  that  the 

triple  (K,A  ,B  ) has  transmission  zeros  coincide  with  q..  Since  we  have 
o o n j 

shown  that  the  transmission  zeros  coincide  with  (A^  - ) , it  is 

reasonable  to  expect  that  the  placement  of  the  slow  eigenvalues  can  be 
achieved  by  solving  a pole  placement  problem.  To  show  that  this  is  indeed 
the  case,  let  M be  the  transformation  defined  by  (2.9)  and  x = [xj  x9 ' s Mx. 
Then,  (2.1)  is  transformed  into 


X1  = AllXi  + A12X2 
X2  ~ A 2 1 x 1 + A22X2  + u 


(2.51) 

(2.52) 


Lemma  2 . 3 : 

If  the  pair  (A^.B^)  is  controllable  ( stabil izab le ) , then  the  pair 

o o 

(AU,A  ,)  is  controllable  (stabilizable) . 


Proof 


The  controllability  of  (Aq,Bo)  implies  [31] 


rank 


cl  - A 


11  12 


r-  T - \°  R° 

z 1 w 1 A22  °2-l 


(2.53) 


and,  since  rank  B9  = m, 


rank  r c I - A 


11 

for  all  complex  number  c.  Thus,  (A?, ,A°,)  is  controllable.  When  (A  ,R  ) 

111-^  o o 

is  only  stabilizable,  consider  u = v - B9  A9^x,.  Since  \;A99,B9)  is 
controllable,  rank  cI-A^  A°9_  < n_,n  antl  uncontrollable  but  stable 


A12]  = n - m 


(2.54) 
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eigenvalues  of  Aq  are  eigenvalues  of  A°^.  Hence  (A^.A^)  is  stabilizable. 

Since  bv  assumption,  (A  ,B  ) is  controllable  from  Lemma  2.3,  there  exists 

o o 

an  (n-m)xrn  matrix  K such  that 
v s 

(A 1 1 " Ai2Ks)  = qj’  J"1'  •***  n-m  ' (2.55) 

By  defining  Kg  * we  have  shown  that  the  slow  eigenvalues  can  be 

placed  by  solving  a (n-m)-dimensional  pole  placement  problem.  By  observing 
that 

\i(KBo)  = ^(KjBj)  * \i(B°K2)  , (2.56) 

the  fast  eigenvalues  are  placed  by  solving  a ra-dimensional  pole  placement 
problem.  By  defining  and  letting  = K^Kg,  we  recognize  that 

the  matrix  Kg  and  which  place  the  slow  and  fast  eigenvalues  respectively 
are  designed  independently. 


Theorem  2.3: 

Let  (A  ,B  ) be  a controllable  pair  and  let  the  (n-m)  x m 
o o 

matrix  Kg  and  m x m nonsingular  matrix  K.  be  such  that  (2.55)  and 


Xi(KfB2)  “ Pi* 


i-1. 


(2.57) 


holds.  Then,  the  use  of  the  high  gain  feedback 


u * g[  KfKg  xj  + Kf  x^] 


(2. 53) 


places  the  eigenvalues  of  the  system  (2.51),  (2.52)  to  q . + 0(— ) and 

J 8 

SlP  t + 0(1)]  . 

Proof : Substituting  (2.58')  in  (2.51),  (2.52)  and  applying  the  transfor- 
mation x = F x1  where  T is  as  in  (2.17)  yields 





— gggj 


wmm—mmm 
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(A11  ' A12Ks)XL  + AL2Kf  X2 


(2.59) 


Ux  2 “uKf  lKsA11  " A22Ks  " KsA12Ks  + A12  xl 

-1 


+ KfB°[l  + MB°  (Ks +A22Kf1)^X2  * 


(2.60) 

The  proof  follows  from  Lemma  2.1  and  Theorem  2.1  and  (2.56). 

From  the  above  theorem,  a procedure  for  completely  separated 
placement  of  the  slow  and  fast  eigenvalues  can  be  summarized:  design 


o o 

K.  to  place  the  eigenvalue,  of  A,,  -A.„K  and  design  Kr  to  place  the  eigen- 
s ° 1112s  f * 


values  of  R#>  then  form  the  high  gain  feedback  control  as  in  (2,58), 
As  a simple  illustration  consider  a system  in  the  form  of 
(2.51),  (2.52), 


3 

1 

' — 

1 

r 

0 

0 

-6 

1 

0 

X + 

l 

1 

0 

0 

3 

0 

1 

—1 

(2.61) 


First,  suppose  that  a high  gain  feedback  control  is  to  be  found  in  place  of 


the  slow  eigenvalue  near  = -3  and  the  fast  eigenvalues  near 


2,3 


g(-l  + jl).  We  solve  the  two  lower  order  pole  placement  problems 


(2.55),  (2.57)  and  obtain 


~6" 

0 f 

__ 

. K = 

’ f 

_0_ 

_p2  0_ 

(2.62) 


The  high  gain  feedback  (2.58)  is 


r o 


-12 


-2 


(2.63) 


1 


L9 


Second,  consider  the  high  gain  feedback  system  (2.61),  (2.63)  with  (2.63) 

written  as  u ■ g C^x  and  apply  Theorem  2.1.  The  closed  loop  eigenvalues 

2 

are  computed  approximately  by  neglecting  the  y terms  in  (2. 23),  (2.29), 
that  is, 


X1  " X(Fn"uF12Lo)  “ X('3  + 30ij)  " "3  + 30  u 


(2.64) 


X2,3  ' ^ X<CoBo  + V * Z 


+ .!] +u 


- (-1  + 5u  + j Jl  + 16 u - 4)i 2 ) 


(2.65) 


The  upper  bound  in  (2.27)  is  ■ 0.0074  corresponding  to  g^  * 135.6.  The 

exact  eigenvalues  are  computed  for  comparison  purpose.  Denoting  the  error 

between  the  approximate  and  exact  eigenvalues  as  e and  between  the 

app 

specified  and  exact  eigenvalues  as  *SpC»  the  results  are  summarised  by  the 
following  table. 


~1 


2.-4.  Near  Optimum  !1  i gh  Cain  Rogul  a cars 

High  gain  feedback  systems  can  also  result  for  the  optimization 
of  system  (2.1)  with  respect  to  a quadratic  performance  index  having  small 
penalty  on  the  control  variables  u, 


1 ■ T *>  t n 

J = ~ „ [x  Qox+y“u  Ru^dt 

o 


(2.66) 


where  is  symmetric  positive  semide finite , R is  symmetric  positive  definite 
and  p is  a small  positive  scalar  parameter.  This  class  of  problems  belong  to 
the  so  called  "cheap  control"  problems  which  have  been  considered  in 
^11,12,13,32*33].  Detailed  results  exist  [ll]  for  the  case  when 


B Q B > 0 . 
o o o 


(2.67) 


Under  this  condition,  the  resulting  optimal  high  gain  feedback  system  satisfies 
the  assumption  (2.1s).  Thus  the  analysis  of  section  2.1  is  applicable  to 
reveal  the  two  time  scale  property  of  the  optimal  state  regulator  developed 
in  [ll].  We  assume  that  the  transformation  M has  been  applied  to  (2.1) 
and  we  consider  the  open  loop  system  (2.51),  (2.52).  Correspondingly,  the 

.1  T _ 1 

weighting  matrix  on  the  states  becomes  Q = (M  ) Q M and  we  denote 

o 


Q11  Q12 


^12  ^22 


(2.68) 


where  and  Q9,,  are  (n-m)  * (n-m)  and  m x m matrices  respectively.  The 
condition  (2.67)  becomes 


B2  Q22  B2  > ° • 


(2.60) 


Since  B„  is  nonsingular,  (2.69)  implies  Q.,.,  is  nonsingular.  We  now  summarize 
the  results  in  r ll]  with  the  next  theorem. 
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Theorem _2  . : 

Suppose  condition  (2.69)  holds,  the  pair  (A^  - A^Qo^’Qp  ,A°,) 
is  stabilizable  and  the  pair  (P,A^  - A^Q^Qp ) is  detectable  where 


0TV 


<11 


" ^12Q22Q12 


(2.70) 


then  the  optimal  Riccati  solution  P possesses  a power  series  expansion 
tor  t>  *=  0, 


P 


‘ l 


i _ 

M P12iU 


1H 

,t  a+i 


*122."  22£P 


H+l 

2+1 


(2.71) 


where  and  P1>(,  are  (n-m)  x (n-ra)  and  in  x m matrices  respectively,  and 

the  zeroth  order  term  P q ts  the  unique  symmetric  semidefinite  solution  of 
the  matrix  Riccati  equation. 


o -1„T 


1..T 


110(AirA!2  Q22Q12}  + (An-Ai.-Ai^>X->)  p 


11  12  12s22x12 


110 


+ V V - P110A12A12^22A12  P110’ 


P.,9q  is  a symmetric  positive  definite  matrix. 


-1/2  -l/7 

P220  - S lS 


Q,->s 


e-l/2,l/2g-l/2 


with 


S 


(2. 72) 


(2.73) 


(2. 74) 


and  ,q  satisfies 


120 


^P110A12+^12)P220S  1 


(2. 75) 


Moreover,  the  optimal  feedback  control  is  given  by 


-u  R-1BoT  (I  [pTnx*  + 

i=o 


(2.76) 


From  Theorem  2.1,  we  see  that  the  eigenvalues  of  the  optimal  high  gain 
feedback  system  (51),  (52),  (76)  are  approximately 


X!  “ y [V"K  le2Tp220B2)  + 0(m)]’  1 = 1 


(2.77) 


Xj  = Xj(All  ' A12P220P120)  + 0(U)*  J c 1 


, . . . , n-m 


(2.78) 


Since  the  matrices  R and  P9 ,)(j  are  positive  definite  and  K1^  is  non- 
singular,  the  matrix  -R  is  negative  definite.  Hence,  the  root 

locus  asymptotes  form  an  angle  of  130°  with  the  positive  real  axis  in 
the  complex  plane. 


From  the  fast-slow  separation  of  the  pole  placement  design,  it 
can  be  expected  that  a similar  decomposition  is  feasible  in  the  near 
optimum  high  gain  regulator  design.  Instead  of  the  full  order  problem, 
the  following  two  reduced  order  regulator  problems  are  solved. 


Problem  "s": 


(2.80) 


Optimize  system  (2.71)  with  respect  to  the  performance  index 


, VJ  , , Cl 

x = A, , x + A, „u 
3 11  s 12  s 

0CO 

JS  " 2 j (xsqilxs  + 2xsQ12Us  + V'22us)dt- 


(2.79) 


(2.30) 


23 


Problem  "f": 

Optimize  system  (2.81)  with  respect  to  the  performance 
index  (2.82) 


Xf  = B°uf 

(2.81) 

00 

Jf  " 1 \ (xf^22Xf  + ufRuf>dt- 

(2.82) 

o 


In  (2.79),  (2.81),  matrices  A°^,  A°2  and  B°  are  as  in  (2.51),  (2.52) 
and  xg  is  an  (n-m)-vector , x^.u^  and  ug  are  m-vectors.  In  (2.80), 
(2.82),  matrices  are  as  in  (2.68)  and  matrix  R is  as  in  (2.66). 


Lemma  2.4: 


If  the  pair  (A  ,B  ) is  stabilizable  and 
o o 


(An,P)  is  detectable 


(2.83) 

where  V is  defined  in  (2.70)  and  if  (2.69)  is  satisfied,  then  there 


exists  a unique  stabilizing  solution  P (which  is  symmetric  positive 

s 

semidef inite)  if  the  matrix  Riccati  equation 


Ps(All  " A12Q22Q12)  ~ (A11  “ A12Q22Q12)  Ps 


, p a°  p 0^0 

s 12Q22  12  s 


(2.84) 


and  the  optimal  control  for  Problem  "s"  is 


“s  ■ - <>22  «I2  + - G.V 


(2.85) 


Proof:  By  Lemma  2.3,  the  pair  (A°^,  A°?)  is  stabilizable  and  hence 
(A°j  - A^Q^Q^.A^)  is  stabilizable.  Together  with  (2.83),  the 
existence  and  uniqueness  of  Pg  is  guaranteed  [ 34  J . 

Since  B°  and  are  nonsingular,  the  optimal  control  for  problem  "f"  is 


- . • _ 


lx  i 


TL3 


readily  obtained  as 


“'‘“IVf 2 Vi 


(2.86) 


where 


P . 
t 


•1/2  1/2  1/2  l/'1  -1/2 


0. 


(2.87) 


A near  optimal  control  for  the  system  (2.51),  (2.52)  is  now  comprised  of 
the  controls  us  and  u^.  in  the  form  of  the  control  law  (2.58) 


u.-  ■ *<cfv; + V:> 


(2.88) 


It  is  crucial  to  note  that  the  feedback  gain  matrices  are  operating  on  the 
state  variables  x’  and  x.'  of  the  system  (2.51),  (2.52)  rather  than  on  x 

t 4.  S 

and  xf.  The  feedback  system  (2.51),  (2.52),  (2.88)  satisfies  Theorem  2.3 

with  G and  G,  playing  the  roles  of  K and  K,,  respectively.  Thus,  J 
s r si  s 

and  J^.  can  be  interpreted  as  the  performance  specifications  for  the  slow 
and  fast  modes. 


Theorem  2^5 : 

Under  the  conditions  of  Lemma  2.4,  the  performance  J of  the 
teedbaek  system  (2.  51),  (2.52),  (2.88)  is  0(1.1*')  close  to  its  optimum 
performance . 

Proof : Since  (2.84)  and  (2.S7)  coincide  with  (2.72),  (2.73),  the  leading 
terms  and  P,,,0  of  the  optimal  Riccati  solution  expansion  (2.71)  are 

identical  to  P^  and  P , respectively.  From  the  identity  v. 2 . 75 ) which  is 


1 


equivalent  to 


i 


it  follows  that  u is  0(u)  close  to  the  optimal  control  (2.76),  that  is, 
c 


uc  - u + 0(u)  (2.90) 

and  the  states  x of  the  feedback  system  (2.51),  (2.52),  (2.88)  is  0(u) 
close  to  the  states  of  the  optimal  feedback  system  (2.51),  (2.52),  (2.76). 
By  Theorem  3 in  [35],  the  performance  of  (2.51),  (2.52),  (2.88)  is 

O 

0(u‘‘)  near  optimal. 

The  results  on  cheap  control  [11]  are  thus  recast  in  terms  of 

2 

two  separate  problems  which  give  an  0(y  ) near  optimal  solution.  In 

view  of  Lemma  2.4,  we  note  that  the  stabilizability  assumption  in  Theorem 

2.4  is  equivalent  to  the  assumption  that^A^.B^)  be  stabilizable  which  is 

the  standard  assumption  made  in  linear  quadratic  regulator  problems.  The 

T 

standard  assumption  on  the  detectabilitv  of  (C,A  ) with  C C ■ Q is, 

o o 

however,  replaced  by  detectability  assumptions  (2.83)  and  (2.69)  in  cheap 
control  problems  which  can  be  interpreted  as  the  detectability  of  the 
slow  and  fast  subsystems  respectively. 

If  instead  of  a state  regulator  problem,  a output  regulator 
problem  is  to  be  solved,  that  is,  (2.66)  is  replaced  by 


J 


00 


T*  1 T 

[z  z + y“u  Ru]dt 


(2.91) 


where  the  output  is  denoted  bv  an  m-vector,  z = C x.  If  rank  C = m, 

oo 

T 

and  C B is  nonsingular,  then  bv  making  the  substitution  Q * C C in 
oo  ooo 

2 

(2.66),  an  O(u^)  near  optimal  control  can  be  found  by  the  preceding 
decomposition  procedure.  The  matrix  Q of  (2.68)  has  the  form 


rpw — ’ 


i. 


o o 

'll  '12 

Q 2> 


2b 


“ ,T  ,,T  , ' 

I C1  CjC, 

i _T  , _T 

! *-  ■•>  ^ ^ 2 ^ i 


(2.92) 


where  C , C,  .ire  as  in  (2.12).  For  Problem  "s"  and  "f",  introduce 


y >=  c x 

s 1 s 


yf  - C2xf 


(2.93) 

(2.9/,) 


as  the  output  variables  of  system  (2.79)  and  (2. SI),  respectively. 

Then,  these  two  problems  become  output  regulator  problems.  For  the 
matrix  Q given  by  (2.92),  V * 0 in  (2.70).  Thus,  we  have  the  following 
lemma . 


Lemma  2.5: 


If  the  pair  (A  ,K  ) is  stabilicable  and 
o o 


/ . o 


R°X(A11  “ A12C21(V  " 0 


(2.95) 


then  the  unique  stabilizing  solution  P of  the  matrix  Rlceati  equation 
(2.89)  is 


P - 0 
s 


(2.96) 


and  the  optimal  control  for  Problem  "s"  is 


„-l„ 


u = - C,  C.x  = G x . 
s 2 1s  s s 


(2.97) 


Proo f : This  lemma  is  a consequence  of  Theorem  13  in  [36]. 
The  composite  control  u in  this  case  becomes 


e = g(R‘13t>t.c;1C,x’  + K~V>.x', 

^ -l  - II  4-tJ 


(2.98) 


and  the  slow  eigenvalues  of  the  optimal  and  near  optimal  feedback  systems 
can  be  approximated  by 

Xj  " Xj(All  ' A12C21C1)  + °(U)’  j = (2.99) 

From  Theorem  2.2,  it  is  shown  that  (A°^  - A°?^2^C1^  are  C^e  transmissi°n 
zeros  of  the  triple  (CQ,  AQ,  BQ) . Hence,  we  conclude  that  if  all  the 
transmission  zeros  of  (CQ,  AQ,  BQ)  lie  in  the  open  left  half  complex 
plane,  then  the  finite  eigenvalues  of  the  optimal  feedback  system  tend  to 
these  transmission  zeros  as  the  small  penalty  on  the  control  variables 
U 0.  The  condition  (2.95)  is  analogous  to  the  condition  on  transmis- 
sion zeros  for  achieving  "maximal  accuracy"  in  [12]. 
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CHAPTER  d 

HIGH  GAIN  FEEDBACK.  FOR  SYSTEMS  WITH 
UNCERTAINTIES  AND  INACCFSSIBI E STATES 

It  is  wo l 1 known  from  classical  control  theory  that  high 
gain  loops  arc  capable  of  reducing  sensitivity  of  single  input  control 
systems  with  respect  to  parameter  variations  and  disturbances.  In  the 
same  spirit,  the  properties  of  multivariable  linear  systems  with  high 
gain  loops  are  under  investigation  in  this  chapter.  The  system  con- 
sidered is  of  the  form 


x « (A  + <SA)  x + (H  f 
o o 

5B) u + E f 

0 

id. 

i) 

v * (0  + £C)x 

o 

id. 

:) 

2 * D X 

o 

id. 

d) 

where  the  state  x € Rn,  the  control 

input  u t N , the  me a 

sured  out 

put 

v £ K\ 

the  disturbance  f c Rk  and 

the  regulated  variable 

s /.  c RC. 

1 1 is 

assumed 

: that  the  nominal  values  of 

the  system  parameters 

Ao * V 0 

o 

and  E 

0 

are  known  and  the  matrix  P 

0 

is  also  known.  Kurthe 

rmore,  the 

d i s t urh 

ance  f is  not  measurable  and 

the  onlv  information 

about  it 

is 

t hat 

f k hj  and  f - h,  for  some  positive  constants  h^  and  h,. 

1'he  following  problems  are  addressed.  Given  the  measurements 
vit).  under  what  conditions  does  a feedback  control  exist  such  that 
desirable  attenuation  of  the  effects  of  disturbances  and  parameter 
variations  on  the  regulated  variables  /.it)  is  achieved?  What  are  the 
svnthesls  procedures  for  such  a feedback  control.' 


w 
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The  class  of  feedback  controls  which  we  conjecture  to  have  the 
desirable  properties  is  high  gain  observer  feedback 

u = gKS  (3.4) 

where  g is  a large  scalar.  This  calls  for  a full  order  observer  which 
has  the  property  that  the  estimation  error 

e ' = x - ft  (3.5) 

between  the  states  x and  observer  states  $ can  be  reduced  to  zero 
asvmtotically  in  the  presence  of  parameter  variations  and  disturbances. 

A prospective  structure  for  an  observer  with  such  capability  is  one  that 
itself  contains  high  gain  loops,  that  is, 


A x + B u - Hv 

(3.6) 

O 0 

g (y  - cox) 

(3.7) 

where  g is  a large  scalar.  The  closed  loop  system  (3.1)-  (3.4),  (3.6)  - 
(3.7)  contains  two  groups  of  high  gain  loops  and  the  design  parameters 
are  the  matrices  K and  H. 

The  organization  of  this  chapter  is  as  follows.  First,  we 
analyze  the  properties  of  system  (3.1)  with  high  gain  state  feedback 

u = gKx  (3.8) 

whre  K is  a given  matrix  appropriately  designed,  for  example,  by  the 
methods  in  Sections  2.3,  2.4.  A condition  for  x(t)  to  have  arbitrarily 
small  sensitivities  with  respect  to  parameter  variations  and  disturbances 
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is  obtained.  Similarly,  the  analysis  is  then  applied  to  the  high  gain 
observer  system  (3.6),  (3.7)  for  a given  matrix  H and  a condition  for 
e\t)  to  have  arbitrarily  small  sensitivities  is  obtained.  Thirdly, 
conditions  arc  derived  for  the  regulated  variables  z(t)  to  have  small 
sensitivities  when  full  state  high  gain  feedback  (3.8)  is  permissible. 
Then  we  develop  conditions  when  high  gain  observer  feedback.  (3.4)  is 
applied.  Finally,  we  look  into  the  synthesis  problems  of  matrices  K and 
H that  satisfy  the  above  conditions  and  at  the  same  time  place  the 
eigenvalues  of  the  main  loops  and  the  observer  loops. 


3 . 1 High  Gain  State  Feedback 

The  behavior  of  system  (3.1)  when  high  gain  state  feedback 
(3.4)  is  applied  is  examined.  We  assume  that  the  feedback  matrix  K 
is  designed  such  that  the  closed  loop  undisturbed  nominal  system,  that 
is,  (3.1),  (3.4)  with  5A,  OB,  5C  and  Eq  equal  to  null  matrices  of 
appropriate  dimensions,  possesses  a two  time  scale  property  and  desirable 
characteristics.  As  it  is  shown  in  the  last  chapter,  we  can  assume  that 
the  transformation  M of  (2.9)  has  been  applied  to  (3.1)  and  in  the  new 
coordinates  (3.1)  is  represented  by 


x’  = (A°n  + 6An)x;J  + (A°n  + 5A10)x^  + OB x u + eJ  f (3.9) 

x'  = (A°x  + 5An)x:  + (A°,  + 6A^)x^  + (S°  + C'B^u  + E°f  (3.: 


T t T 

where  x = [x|  x9 ] , x^  is  an  (n-m)  vector,  x’  is  an  m-vector  and, 

B^  is  a mxm  nonsingular  matrix.  We  adopt  the  convention  that  matrices 
with  superscript  o and  prefix  5 denote  nominal  values  and  parameter 
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variations,  respectively.  Let  the  high  gain  state  feedback  (3.4)  in 
these  coordinates  be 


u - g[KfKsx|  + Kfx^; 


as  in  (2. 53).  Our  assumption  on  K then  becomes 


Re  A(A?-A°  K ) < 0 
11  12  s 


and 


Re  \(KB  ) < 0. 
o 


(3.11) 


(3.12) 


(3.13) 


Introducing  the  variable 


U)  = K,K  x^  + K.x' 
f s 1 f 2 


(3.14) 


the  condition  for  x(t)  to  have  arbitrarily  small  sensitivities  with 
respect  to  parameter  variations  and  disturbances  is  now  stated  in  the 
following  theorem. 

Theorem  3.1: 

The  solution  of  the  closed  loop  system  (3.9)-(3.11)  can  be 
approximated  by 


(A11  ‘ A12Ks)c  -1 

xL(t)  = e 1 x1(0)  + 0 (g  i) 

gK(B  + 5B)t  _ 

oo(t)  = e 0 co(0)  + 0(g-i) 


(3.15) 

(3.16) 


if  and  onlv  if 


R(E)  q R(B), 
o o 


(3.17) 





1 


R(6 A)  C R(Bq) 


R(6B)  c r(b  ) 
o 


(3.18) 


(3.19) 


Re  A [ K ( B + 6B)  ] < 0. 
o 


(3.20) 


Moreover,  let  CTg  = max  ReA[K(BQ  + 5B)  ] and  tg  = jaj  u^np.  Then  for  all 
t ^ tg,  u(t)  is  reduced  to  0(g  l) , that  is,  joj(t)|  = 0(g_1). 

Proof:  In  terms  of  the  new  variable  u),  the  closed  loop  system  (3.9)— (3. 11) 
has  the  form 


x'  = (Z°  + &Z1 ) + (Z°  + 6Z2  + goB^oj  + E°f  (3.21) 

J>  = (F°  + 6F  ) x'  + [F°  + 5F  + gK(B  + oB)  ]gj  + RE  f (3.22) 

1 L X Z Z O O 


where 


z°i  ■ Aii  - AnKs 


7°  - 1°  V-1 

z2  “ A12Kf 


F°  = W°n  + KfA°x  - F°KfKg 


F2  = (KfA22  + KfKsA12)Kf1 


(3.23) 


(3.24) 


(3.25) 


(3.26) 


and  the  matrices  6 Z^,  5F^  are  defined  as  in  (3.23),  (3.24)  and  (3.25), 
(3.26)  with  replaced  by  SA„ . We  first  prove  sufficiency.  Condi- 
tions (3.17)  — (3. 19 ) imply  that  6B^  = 0,  E°  = 0,  SA^  = 0 and  5A°0  = 0. 
Consequently,  5Z^  = 0 and  6Z0  = 0.  Using  Lemma  C2 , we  obtain  (3.15). 


T 
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Applying  the  variation  of  parameters  formula  on  (3.22),  we  have 


K(B  + SB)— 
o U 


,C  K(B  + 6B) (t-r) 


ui(t)  « e 


u)(0)  + 


M KE  f (t) dT  + 0(u).  (3.27) 


Integrating  the  integral  by  parts,  (3.27)  becomes 


K(B  + SB)— 
o U 

= e cj(0)  + 


(3.28) 


K(B  +SB) 


(t-T) 


[K(B  + (SB)  ]_1e 
o 


KE  f(T)dT  + 0(u) 
o 


Since  |f|  - h^  by  assumption  (3.28)  becomes  (3.16).  Necessity  can  be 
proved  by  contradiction.  Suppose  condition  (3.17)  is  violated,  then 


^ 0 and 

x^(t)  = e x^(0)  + 


r'  z°<t-x)  0 

e 1 E^f (t ) dx  + 0 (p) . 


(3.29) 


Clearly,  the  integral  in  (3.29)  is  not  0(u).  Now,  suppose  conditions 

(3. 18) (3. 19)  are  violated,  then  SB,  # 0,  5Z  / 0 and  5Z,,  # 0.  Lemma  C3 

ii  - 

gives 


x^(t)  - e 


(Z,  + SZ,  + SB  L)  t 


[xl  (0)  - SB  [K(B  + SB)  ] u)(o  ) ] 
1 1 o 


+ SB,  [K(B  + SB)  ] co(t) 
1 o 


r (Z°  + SZL  + SB,L)(t-r). 
e 


SB,  [K(B +5B)  ]_1KE  f(T)dT+0(u) 
1 o o 


(3.30) 


where  L = [K(B  + SB)]  ^ ( F"?  + 5F,).  This  completes  the  proof, 
o 11 
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The  following  corollary  is  evident  from  the  proof. 


Corol larv  3.1: 


The  solution  of  the  systems  (3.9)— (3. 11)  can  be  approximated 


(Z  + 6Z  + 5B  L)t 

x!(t)  = o 1 [x!(0)  - 5B,(K(B  + 5B))_1a)(0)] 

1 1 1 o 


+ 5B,  (K(B  + 5B) ) U)(t)  + 0(U) 

1 o 


(3.31) 


and  (3.16)  if  and  only  if 


Rl(I  - 6 B ( K ( B + 5B))_1K)E  ] c R(B  ) 
o o—o 


(3.32) 


For  sufficiently  small  y,  this  system  is  asvmtotieally  stable  if  and 
on 1 v if 


Re  A [Z^  + 6ZX  + 5B1L]  < 0 


(3. 33) 


and  condition  (3.20)  holds  where  the  eigenvalues  of  the  matrix  + 6Zj 

+ 53,1.  are  the  transmission  zeros  of  the  triple  (K,  A + 5A,  B + 5B) . 

1 o o 

The  above  theorem  and  its  corollary  show  that  the  sensitivities 

of  the  high  gain  state  feedback  system  (3.9)— (3. 11)  with  respect  to 

parameter  variations  5A,  5B  and  disturbances  f depend  on  the  ranges  of 

5A,  5B  and  E . Wien  the  ranges  of  these  matrices  all  lie  inside  the 
o 

range  of  Bq,  the  motion  x|(t)  has  0(y)  sensitivities  and  the  influences 

of  disturbance  and  parameter  variations  to  o)(t)  for  t 2 t are  also 

s 

0(u).  The  stability  of  the  fast  transient  in  ui(t)  in  the  "stretched" 


time  scale  gt  depends  on  5B.  By  an  earlier  assumption  (3.13),  there 
exists  5B  with  || 5 B 1 1 sufficiently  small  such  that  the  fast  transient 


. 


35 


decays  exponentially.  Thus,  we  say,  the  stability  of  the  fast  motion 
u)(t)  is  robust  in  the  sense  of  (3.20).  If  the  range  conditions  (3.17)- 
(3.19)  are  not  satisfied  but  condition  (3.32)  holds,  then  the  sensitivity 
of  the  motion  x^(t)  with  respect  to  f is  0(p).  From  (3.31),  the  motion 
xj(t)  is  no  longer  "slow"  and  is  influenced  by  the  system  parameter  vari- 
ations 6A  and  5B.  Furthermore,  the  stability  of  system  (3.9)— (3. 11) 
is  only  robust  in  the  sense  of  (3.20)  and  (3.33). 

As  an  illustration,  consider  the  3rd  order  system  given  by 
(2.61)  in  the  Section  2.3.  Suppose  this  is  the  nominal  undisturbed  sys- 
tem. If  the  columns  of  the  disturbance  matrix  E^  and  the  parameter  vari- 

T 

ation  matrix  5A,  5B  are  of  the  form  [0  5^  + 6.,  6^]  for  any  values 

of  6^  and  , then  the  range  conditions  ( 3. 17)- ( 3. 19)  are  satisfied.  If 
the  first  element  of  any  of  the  columns  is  nonzero,  then  ( 3. 1 7)- ( 3 . 19) 
do  not  hold. 

3. 2 High  Gain  Observer 

The  objective  of  the  high  gain  observer  (3. 6), (3. 7)  is  to 
reduce  the  estimation  error  e'  to  zero.  Expressing  (3. 6), (3. 7)  in  terms 
of  e ' , that  is. 


e'  = A e'  + Hv  + iSAx  + 6Bu  + E f 
o o 


v » g(C  ^e ' + 6Cx) 


(3.34) 

(3. 35) 


it  is  clear  that  (3.34),  (3.35)  represent  a high  gain  feedback  system  with 
disturbances  x,  u and  f.  The  design  parameter  is  the  "input"  matrix  H 
while  the  feedback  matrix  is  specified.  Thus,  the  high  gain  observer 
design  becomes  the  design  of  matrix  H such  that  the  system  (3. 34) ,( 3. 35) 


I 


r.  ' 


is  asymptot ical ly  stab le . This  design  problem  can  be  viewed  as  the  dual 


problem  of  the  high  gain  state  feedback  design  problem  solved  in  Section 


2.3.  This  duality  is  set  forth  by  the  following  theorem. 


Theorem  3.2: 


Let  (A^.C^)  be  an  observable  pair,  then  there  exists  a matrix 


H such  that  the  slow  and  fast  eigenvalues  of  the  high  gain  feedback  svs- 


e'  = (A  + gHC  )e' 
o o 


(3.36) 


are  placed  to  + 0(g  *) , j = l,...,n-r  and  g[pA  + 0(g  1) j i = l,...,r 


where  q^  and  gp^  are  the  prescribed  locations  of  the  slow  and  fast  eigen- 


values, respectively. 


T T 

Proof:  Bv  the  Duality  Theorem  of  Kalman,  (A  , C ) is  a controllable  pair. 
' o o 


T /%  T T 

Since  A(Aq  + gC^H  ) = \(A^  + gHC  ),  the  pole  placement  problem  for  (3.36) 


T T T 

is  equivalent  to  finding  the  feedback  matrix  H for  the  pair  (A^,  C ) . 


By  Theorem  2.3,  the  theorem  is  proved. 


In  the  rest  of  this  section,  we  assume  that  the  design  of  matrix 


H by  the  methods  presented  in  Section  2.3  has  been  carried  out.  Since 


system  (3.34),  (3.35)  is  similar  to  the  high  gain  state  feedback  system 


( 3. 9)- ( 3. 11) , it  is  to  be  expected  that  the  conditions  for  e' (t)  to  be 


insensitive  to  the  disturbances  involve  the  ranges  of  5A,  6B,  E and  H. 

o 


Before  we  state  the  theorem,  it  is  convenient  to  represent  the  system 


(3.34),  (3.35)  in  the  new  coordinates  defined  by 


r —i 

ei 


(3.37) 
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where  e^  is  an  (n-r)-vector , e.,  is  an  r-vector  and  ft  Is  defined  by 


and  H,  is  an  r x r nonsingular  matrix.  Thus  we  have 


(3.38) 


•1  ' ^iiei  + lh°2  * SV  + SV  + E°‘ 


(3.39) 


e,,  * e^  + A°0e0  + H^v  + iSA.,x  + 5B.,u  + E^f 


(3.40) 


v “ + C°e2]  + g6Cx 


(3.41) 


where  C“  is  an  r x r nonsingular  matrix.  Our  assumption  on  H is  then 


Re\  a°n  - A°2C°  C J)  < 0 


(3.42) 


ReA  (C  H)  < 0 
o 


(3.43) 


The  theorem  is  stated  in  terms  of  the  variable  n 


„o  , _o 
H * C1e1  + C^e0. 


(3.44) 


Theorem  3. 3: 


The  solution  of  the  closed  loop  system  (3.37),  (3.3S)  can  be 


approximated  by 


/a°  a°  r°  r°U 

^11  Ai'>C->  ^i-  c ^-X 

eL(t)  - e 11  L~  1 eL(0)  + 0(g  l) 


(3.45) 


gC  Ht  , 

n(t)  - e L n(°)  + o(g  ) 


(3.46) 
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if  and  onlv  if 


R(SA)  c K(H) 
R(5B)  c R (H) 


R(E  ) C R(H) 
o - 


5C  = 0. 


(3.47) 


(3.48) 


(3.49) 


(3.50) 


Moreover,  let  3,  = max  Re\(C  H)  and  t , = | CT  | ^jjinu  with  u = g 
b o » » 

Then,  for  t ^ £ , rj ( t ) is  reduced  to  0(g  *),  that  is,  |n(t)|  = 0(g  ^) 

Proof : The  proof  is  analogous  to  the  proof  for  Theorem  3.1,  details  are 
omitted. 


3.3.  Regulation  by  High  Gain  State  Feedback 

In  Section  3.1,  we  have  derived  conditions  for  the  states  x to 
have  arbitrarily  small  sensitivities  with  respect  to  parameter  varia- 
tions and  disturbances.  We  recognize  that,  in  reality,  very  often  we  are 
only  interested  in  the  regulation  of  some  of  the  state  variables.  In 
general,  these  variables  can  be  expressed  as  linear  combinations  of  the 
states  which  we  have  denoted  by  z in  (3.3).  We  assume  that  there  are 
l of  these  variables  to  be  regulated,  that  is,  z ~ D^x  is  an  12,-vector  and 

Z < n.  Without  loss  of  generalitv,  it  is  assumed  that  rank  D = l.  In 

o 

this  section,  we  derive  the  conditions  on  D for  the  existence  of  a high 

o 

gain  state  feedback  control  which  regulates  the  variables  z.  Suppose  that 
the  feedback  gain  matrix  R is  designed  such  that  condition  (3.13),  (3.13) 


are  satisfied,  we  have  the  following  theorem. 
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Theoreir  3.4: 

Let  t be  the  time  instant  as  defined  in  Theorem  3.1.  Suppose 

that  the  system  (3.9)— (3. 11)  when  f * 0 is  asymtotically  stable,  then  the 

motion  z(t)  is  reduced  to  an  0(g  quantity,  that  is,  |z(t)|  = 0(g 

for  all  t * t if  and  only  if 
s 


R(DT)  C R(KT). 
o 


Proof : Sufficiency.  Analogous  to  (2.18),  x can  be  written  as 


(3.51) 


x = Qx'  + B (KB  ) a) 
1 oo 


(3.52) 


where  co  is  as  in  (3.14)  and  the  columns  of  the  nx(n-m)  matrix  Q span 


the  null  space  of  K,  that  is. 


KQ  = 0. 


It  follows  that 


(3.53) 


z = D Qx’  + D B (KB  )_1w. 
O 1 0 0 o 


(3.54) 


Theorem  3.1  gives  | co ( t ) | = 0(g  ^)  for  all  t ^ t^,  hence  |z(t)|  = 0(g  ^) 
if  D^Q  - 0 which  is  equivalent  to  (3.51).  Necessity  is  proved  by  contra- 
diction. Suppose  (3.51)  is  violated,  then  DqQ  ^ 0.  From  (3.54),  since 
|x'(t)|  t 0(g  ^)  for  all  t £ t . The  theorem  is  proved. 

T 

In  the  above  theorem,  it  is  shown  that  if  the  range  of  D lies 

o 

T 

in  the  range  of  K , then  z(t)  can  be  reduced  to  zero  as  rapidly  as  desired 
since  t __  is  0(g  ^).  Another  interpretation  of  (3.51)  is  that  the  varia- 
bles z are  linear  combinations  of  the  fast  variables  a).  Since  u>(t)  decays 


exponentially  in  the  time  scale  gt , z(t)  is  reduced  to  zero  in  the  same 
time  scale.  The  following  theorem  states  the  condition  for  z(t)  to 
be  zeroed  in  a finite  time  interval  using  high  gain  state  feedback. 

Theorem  3.5: 

Let  t^  be  a finite  positive  constant.  Under  the  same  assump- 
tion as  in  Theorem  3.4,  the  motion  z(t)  is  reduced  to  an  0(g  quantity, 
that  is,  jz(t)|  = 0(g  *)  for  all  t 2 t if  and  only  if 

< Z°  + <5Z  + 5B  L | R[M  (I  - 6B(K(B  + 6B))_1K)E  ] ~ W(D  Q)  (3.55) 

1111  o o — o 

where  is  as  in  (2.21)  and  the  notations  in  Corollary  3.1  are  used. 

Proof : From  the  proof  of  the  previous  theorem,  we  recognize  that  we  only 

need  to  prove  that  (3.55)  is  necessary  and  sufficient  for  |D^Qx|(t)|  = 
0(g  1)  for  all  t 1 t..  Lemma  C3  shows  that  there  exists  t.  > 0 such  that 

r * IZ°  + SZL  + OB1L](t  - t) 

D Qx.'(t)  = D Q e [E?  - 5B.(K(B  + (SB))-1!^  ]f(r)d 

o i o 1 i o o 

0 

+ 0 (g_1)  (3.56) 

for  all  t - t^..  In  [33],  it  is  shown  that  (3.55)  is  necessary  and 
sufficient  for  the  integral  term  in  (3.56)  to  vanish.  This  proves  the 
theorem. 

J. 

Let  F be  an  n x n and  G an  n x m matrix,  then  <F|R(G)>  denotes 
the  subspace  spanned  by  the  columns  of  the  matrices  G,  FG,  ....  F^  _ ^G 
where  kin. 


'-*■ ---  ■: 


* 


1 


rL. 


t 
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Given  the  variables  to  be  regulated,  condition  (3.55)  character- 
izes the  ranges  of  5A,  5B  and  Eq  which  allow  regulation  to  occur  by  apply- 
ing high  gain  state  feedback.  Nevertheless,  it  is  difficult  to  deter- 
mine from  (3.55)  the  ranges  of  <$A,  6B  and  Eq  that  satisfy  it.  The  follow- 
ing corollary  relates  (3.51)  with  (3.55). 


Corollary  3.2: 


If  the  pair  (Z°  + 5Z.  + $B.L,  E°  - 6B  [K(B  + 6B) ]_1KE  ) is 
l 1 l 1 1 o o 

controllable,  then  condition  (3.55)  becomes  (3.51). 

We  note  that  (3.51)  is  a sufficient  condition  (3.55).  When  there  is  no  in- 
formation available  about  the  ranges  of  fiA,  5B,  then  to  guarantee  regulation 
of  z(t),  we  have  to  take  into  account  all  possible  parameter  variations.  In 
this  case,  we  recognize  from  this  corollary  that  condition  (3.55)  should  be 
replaced  by  the  condition  (3.51).  Thus  (3.51)  becomes  both  necessary 
and  sufficient  condition  for  regulation  of  z(t).  In  view  of  Theorem 
3.4,  z(t)  can  be  reduced  to  zero  as  rapidly  as  desired.  (3.51)  can  also 
be  interpreted  as  a specification  to  the  synthesis  of  matrix  K.  This 
aspect  will  be  brought  up  again  in  Section  3.6. 


3.4.  High  Gain  Observer  Feedback 

We  now  begin  to  examine  the  behavior  of  the  system  (3.1),  (3.2) 
when  high  gain  feedback  is  applied  through  the  high  gain  observer  (3.6), 
(3.7),  that  is,  the  high  gain  feedback  control  is  of  the  form  u = gKx  as 
it  is  given  in  (3.4).  It  is  recognized  that  this  closed  loop  system  con- 
tains two  groups  of  high  gain  loops.  To  facilitate  our  analysis,  we 
express  this  system  in  terms  of  the  state  x of  (3.1)  and  error  e'  between 
x and  observer  state  x. 


— — - 


■w  w: 


! 

I 
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x = [A  + 6A  + g ( B + 6B)K]x  - g(B  + <5B)Ke ' + E f (3.57) 

o o o o 

e’  = (5A  + gH5C  + g5BK)x  + (A^  + gHCQ  - g5BK)e'  + E f (3.58) 

At  this  stage,  we  assume  that  the  gains  of  the  main  loops  and  the 
observer  loops  are  of  the  same  magnitudes,  that  is,  g = Ctg  when  a is  some 
positive  scalar  which  is  0(1).  The  feedback  gain  matrix  K of  (3.4)  is 
designed  such  that  (3.12),  (3.13)  are  satisfied  and  the  matrix  H of  the 
high  gain  observer  satisfies  ( 3 . 42) , ( 3. 43) . Under  these  assumptions,  we 
first  state  the  conditions  for  x(t)  to  have  arbitrarily  small  sensitivi- 
ties with  respect  to  disturbances  f.  Then  we  establish  conditions  for 
arbitrarily  small  sensitivities  with  respect  to  parameter  variations 
<$A,  (SB  and  5C  and  disturbances  f. 

Theorem  3.6: 

Let  S be  the  matrix  defined  in  (D5).  Using  the  notations  in 
Appendix  D,  let 


"B2<K2 

+ ?.,) 

A 

"V 

— 

2 

f 

(3. 

.59) 

_a_1H2 

(C2  - 

<5c  2 ) 

ct-1H 

c ® 1 

T2  j 

'B2<K1 

+ Rx) 

-b9k 

1 

(3. 

, 60) 

-a_1H2 

(Ci- 

5CX) 

a-1H 

2Cl 

r 

All 

0 

AT2 

— 

0 

_ 

1 

* li 

i 

- 

1 z 

l 

) 

(3. 

61) 

_Fn 

A° 

n 

_F12 

A° 

^12 

5 

2 ?1 

Then  the  solution  of  system  (3.57,  (3.58)  can  be  approximated  hv 


>(tl  V R (0)' 


t1<t) 


Lci(0) 


+ 0(g-1) 


(3.62) 


x"  (t)|  g^6.t  fx"  <oTI  -1  R <0) 

2 =e  2 2 + L L 

[e2(t)J  Le2(0)J  21  lfi(0). 


-1  1 (t)  . 

- L L + o(g_1) 

62  61  Lc^oJ 


(3.63) 


if  and  only  if 


R[  (I  - S)E  ] C R(H) 
o — 


R(E  ) c R(B  + 5B) 
o — o 


(3.64) 


(3.65) 


ReX(^  ) < 0. 


(3.66) 


Proof : As  it  is  shown  in  Appendix  D,  system  (3.57),  (3.58)  is  equivalent 
to  the  system  (D8)  - (Dll).  Lemma  C2  gives  (3.63)  and 


x"  (t)  Zit  xi'  (0)  Z r ( t— x)  E . 

=e°  1 + e°  _ f(T)dr  + 0(g_i). 

fl(t)_,  L£1(0)J  LE1. 

0 


(3.67) 


From  (D3) , (D4) , (D5)  and  (D14) , E.^  = 0,  E = 0 providing  that 
(3.64),  (3.65)  holds.  This  concludes  the  proof. 


I 


The  following  Corollary  establishes  the  stability  of  system 

(3. 57), (3. 58). 

Coro  I La  ry  3_._3: 

For  g sufficiently  large,  the  system  (3 . 57) , ( 3. 58)  is 
asymptotically  stable  if  and  only  if  ( 3. 64)- (3 . 66)  and 


ReX(Z.)  < 0 


(3.68) 


holds.  Furthermore,  X(Z^)  are  the  transmission  zeros  of  the  triple 


A + 5A 
o 


5 A 


B + SB 
o 


5B 


K -K 


6C 


and 


X(L  ) = X 


KB 


-a'1KH 


-1, 


6CB  + C 5B  a C H 


(3.69) 


Proof : X(Z^)  and  gX(£~  ) are  the  0(1)  and  0(g)  eigenvalues  of  the  systo 

°2 

( 3. 57) , (3. 58)  , in  the  limit  when  g -►  °°.  The  remainder  of  the  corollary 
follows  from  the  results  in  Sections  2.1  and  2.2.  Details  are  omitted. 


From  Theorem  3.6  and  its  corollarv,  we  conclude  that  if  E satisfies 

o 


(3.64),  (3.65),  then  the  effects  of  disturbances  f on  the  motion  of  system 

-1 


(3.57),  (3.58)  are  attenuated  by  a factor  of  g . The  stability  of  this 
high  gain  feedback  system  depends  on  the  parameter  variations  5a,  cB  and 
5C.  When  there  are  no  parameter  variations,  that  is,5A  =0,  5B  = 0 and 
5C  = 0,  the  fast  and  slow  eigenvalues  as  g 00  are  the  eigenvalues  of  the 


mammmam 


i 1 


matrices 


-1 

KB  -Cl  KH 

o 

o d-1C  H 

o 


(3.70) 


A11  - A12^s  ' a"^fC2  # 


A°  - A°  C°_1C° 
11  12°2  1 


(3.71) 


respectively  where  Kf_,  K^  and  as  defined  in  (3.74)  this  can  be  shown  by 

A 

letting  M.,  = 0 and  Mr  = I in  (Dl)  and  assuming  that  the  x-svstem  has  the 

o o n 

form  ( 3.9) , ( 3. 10) . By  the  assumptions  on  K and  H,  (3.12),  (3.13)  (3.42), 
and  (3.43),  the  eigenvalues  of  the  matrices  in  (3.70),  (3.71)  have 
negative  real  parts.  By  continuity  argument,  \(Z^)  and  A (J^  ) have  nega- 
tive real  parts  providing  that  the  parameter  variations  ||  6A  j|  , | j SB  j | and 
||  5C||  are  sufficiently  small.  Hence,  we  say  that  the  stability  of  high 
gain  observer  feedback  system  is  robust  with  respect  to  5A,  5B  and  5C  in 


the  sense  of  (3.66)  and  (3.63)  when  the  range  of  Eq  satisfies  (3.64), 


(3.65) 


We  now  state  the  theorem  regarding  the  conditions  for  x(t)  to 


have  arbitrarily  small  sensitivities  with  respect  to  parameter  variations 
5A,  5B , i$C  and  disturbances  f. 


Theorem  3.7: 


Let  Z be  the  matrix  defined  in  (3.71)  and  denote 


KB°  + 6B,)Kf 


5B.,Kf 


(Bl,  + 6B2)KfKs  -(b\  + 6B9)KfKg 


(B0  + 5B„)K 


_ 1 O /V  A 

a H-C"  - 6 B , K 
2 o 1 t 


6B0K£K 
2 f s 


a H_C  - iSB-K.K 
2 o 2 f s 


(3.72) 


(3.73) 


where 


-v  '_]  ■'_i  ~ 

K = K M , K = K,.M  , 6B=  / E MiSB. 

9 9 ‘ ‘ SB, 


(3.74) 


Let  the  variable  4*  Be  defined 


* " KfKsXl  + KfX2  ' Kf‘vsCl  " S‘C 2 


(3.75) 


where  x^,  x.'  and  e , e.0  are  the  variables  defined  in  (3.9),  (3.10)  and 
(3.39),  (3.40),  respectively.  Then  the  solution  of  the  high  gain  feed- 
back system  (3.57),  (3.58)  can  be  approximated  by 


x (t) 

1 Z°t 

* e 

e,  (t) 

— * 


x-(°n 


+ D(g-1) 


ei(0)J 


X2<e)  gV°t  fX2(0) 


-1  x,’  (0) 

+ n y° 
l2  “1 


e,(0). 


eL(0) 


(3.7b) 


-r2  r; 


<;(t)~i 


e,(t)  i 


+ o(g'1) 


(3. 77) 


p r 
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if  and  only  if 

R(E  ) c R(B  ),  R(E  ) C R(H)  (3. 78) 

o~o  o — 

R(6A)  c R(B  ),  R(SA)  cr  R(H)  (3.79) 

— o — 

R(6B)  c R(B  ),  R(6B)  c R(H)  (3.80) 

— o — 

SC  - 0 (3.81) 

and 

ReA(£°)  < 0.  (3.82) 


Moreover,  let  o , » max  Re\(T°)  and  t , = la  J ^g  Vlng  ^ . Then  for  all 
of  2 ot  ' of 

t - t £,  4>(t)  is  reduced  to  an  0(g  *)  quantity,  that  is,  |4>(t)|  * 0(g  *)• 

Proof : Sufficiency.  If  SB,  SC  satisfy  conditions  (3.79),  (3.  80),  let 

Mc  * I , ft.  « 0 in  the  transformat  ion  matrix  M defined  in  (Dl)  and  assume 

0 n 5 

that  the  x-svstem  (3.1)  has  the  form  (3.9),  (3.10).  Then  by  definition 

X1  ’ X1  x2  " X2 ’ C1  " V C2  ’ e2  and  Aij  " Aij  + 6Aij*  B2  “ B2  + 6B2’ 

E » E°,  E » f^,  K , C and  SC  are  null  matrices  and  K.  * K,.K^, 

1 11  ill  1.  1 t s 

K,  =Kj.  Condition  (3.78)  implies  5A^  and  SF^  are  null  matrices  and  F.^, 

F.j  are  null  if  (3.77)  holds.  Lemma  C2  gives  (3.  76),  (3.  77).  Necessity 
is  proved  by  contradiction.  Details  are  omitted.  This  concludes  the  proof. 
We  have  the  following  corollary  concerning  stability. 

Corolla rv  3.4 : 

Given  that  the  parameter  variations  fiA,  SB,  SC  and  disturbances  f satis- 
fying (3. 78)-(3.32)  then  for  g sufficiently  large,  the  system  (3.57),  (3.58) 


■ 1 


is  asymptotically  stable. 


Furthermore , 


and  \(Z°)  are  the  transmission  zeros  of  the  triple 


f 

A + <5A 
o 

0 

~B  + 5B 
o 

0 

r~ 

K 

■*) 

\ 

SA 

tmm 

A 

o. 

f 

SB 

H 

, SC 

Co -/. 

(3.83) 


Proof : By  assumption  K and  H are  designed  such  that  (3.12),  (3.13)  and 
(3.42),  (3.43)  holds,  hence  ReA(Z°)  < 0.  The  proof  is  then  analogous  to 
that  of  Corollary  3.3,  details  are  omitted. 

In  Theorem  3.6  and  Theorem  3.7,  it  is  shown  that  the  sensitivity 
of  x(t)  when  high  gain  feedback  control  u = gKx  is  applied  to  the  system 
depends  on  the  ranges  of  the  parameter  variations  matrices  <5A,  SB,  SC, 
the  disturbance  matrix  Eq  and  the  "input"  matrix  H in  the  high  gain 
observer  system  (3.6),  (3.7).  If  the  ranges  of  SA,  SB  and  E^  lie  in  the 
range  of  B^  and  H simultaneously  and  there  is  no  parameter  uncertainty  in 
the  measurement  matrix  C,  then  the  motion  x^(t)  has  0(u)  sensitivity  with 
respect  to  parameter  variations  and  disturbances  and  the  influence  of  f 
and  parameter  variations  to  the  motion  x^(t)  for  t > t^  is  also  0(y). 

The  stability  of  the  fast  transients  in  x^(t)  depends  on  <5B  and  is  robust 
in  the  sense  of  (3.83).  If  the  range  conditions  (3. 78)— (3. SO ) do  not 
hold  and  <5C  ^ 0,  but  conditions  (3.64),  (3.65)  hold,  then  the  sensitivity 
of  x^'(t)  with  respect  to  f is  0(u).  Furthermore,  the  stability  of  system 


49 


(3.57),  (3.58)  is  only  robust  in  the  sense  of  (3.66)  and  (3.68).  From 
(3.67),  we  see  that  if  the  range  of  Eq  lies  outside  the  range  of  H,  then 
e^(t)  is  sensitive  to  disturbance  f and  consequently  x|(t)  is  disturbed 
by  f. 

If  the  gains  of  the  main  loops  and  the  observer  loops  are 

different  in  order  of  magnitudes,  that  is  the  products  gg  1 or  gg  ^ tend 

to  zero  as  g and  g both  tend  to  infinity,  then  instead  of  the  two  time 

scales,  t and  gt  we  have  shown  in  Theorem  3.6,  there  are  three  time 

scales,  t,  gt  and  gt.  We  now  consider  the  case  when  lim  gg  * = 0,  that 

g,g  00 

is,  roughly  speaking,  g » g. 

Theorem  3.8: 

Suppose  that  lim  gg  ^ = 0,  then  the  solution  of  the  system 
g,g  -*■  00 

(3.57),  (3.58)  can  be  approximated  by 

n(t)  = egCoHtn(0)  + (Kg-1)  (3.84) 

gK(B  + 6B)t  A A — 1 

x'(t)  - e ° [x'  (0)  + K~1(KfK  - K,C°  C7)e.  (0)  + K x.'(0)] 

2 2 t t s f 2 1 i si 

-Kf-1(KfKs-KfC°  Cj)ex(t)  - Kgx^(t)  + 0(g_1)  (3.85) 


\(cf 

o 

rt 

= e 

i 

e1(0) 


“I 


+ 0 (g  x) 


(3.86) 


ci(0)J 


if  and  only  if  the  range  conditions  (3. 77) — ( 3 . 30)  and  the  condition  (3.20) 


are  satisfied. 


Proof : Proof  is  analogous  to  that  of  Theorem  3.1,  details  are  omitted. 

The  three  time  scales  property  of  the  high  gain  feedback 
system  (3.57),  (3.58)  when  g » g is  evident  from  ( 3. 84)- ( 3. 86) . The 
motions  n(t),  x9(t),  x^(t)  and  e^(t)  can  be  characterized  by  eigen- 
values of  different  orders  of  magnitudes,  g,\(C  H) , gX(K(B  + 6B))  and 

o o 

X(Z°).  If  the  range  conditions  (3. 77)- ( 3. 80)  are  not  satisfied  but 
conditions  (3.64),  (3.65)  hold,  then  analogous  to  Theorem  3.8,  we  have 
the  following  Theorem. 


Theorem  3.9: 


Using  the  same  notations  as  in  Appendix  D,  let 


T,  = [B9(K9  + k9)  - k2(cqh)  h2(c,  - SC7) 


(3.87) 


\ = - K2(CoH)_1H2C°  K1  + K - K2(CoH)"1H2(C1  - 6C  ) 


F„  + A°  (C  H)_1H9(C.  - 6C  ) 
11  Id  o d l 1 


(3.88) 


F12  + Ai2(CoH)  H:(C2"  6C2l 


p—  1 p 

d 1 1 ’ 


(3.89) 


Then  the  sensitivities  of  the  motions  x(t)  and  e(t)  with  respect  to  f 
are  of  0(g  S if  and  only  if  conditions  (3.64)  and  (3.65)  hold.  Further- 
more, x(t)  and  e(t)  can  be  characterized  bv  the  eignevalues  g\(C  H) , 

o 

g\(r2)  and  \(Z1). 
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Proof:  The  proof  is  analogous  to  that  of  Theorem  3.8,  details  are 
omitted. 

By  comparing  the  solutions  (3.76),  (3.77)  and  (3. 84)- (3. 86) , 

we  recognize  that  for  the  class  of  parameter  sanctions  5A,  SB  and  SC 

and  disturbance  matrix  Eq  satisfying  conditions  (3. 78 ) — ( 3 . 81) , the 

effect  of  having  g » g is  that  the  fast  eigenvalues,  X (£°)  of  (3.72), 

which  are  0(g)  decompose  into  two  groups  of  eigenvalues  of  different 

magnitudes,  gX(C  H)  and  g\(K(Bo  + SB)).  Moreover,  from  Theorems  3.1  and 
o 

3.3,  these  eigenvalues,  gX(CoH)  and  gX(K(Bo  + 5B)),  are  the  "fast"  eigen- 
values of  the  high  gain  observer  system  (3. 39)-(3.40)  and  the  high  gain 
state  feedback  system  (3.9)-(3.11) , respectively.  Furthermore,  the  0(1) 
eigenvalues  of  the  high  gain  observer  feedback  system  (3.57) , (3.58)  are 
the  union  of  the  sets  of  0(1)  eigenvalues  X(A°  - A°  K ) of  (3.9)— (3. 11) 

IX  X 1Z  S 

and  X(A°^  - A^C^  ^1  ^ (3. 39)-(3. 41)  which  is  evident  from  (3.71). 

Thus  we  conclude  that  the  separation  property  of  regular  gain  observer 
feedback  systems  holds  if  conditions  (3. 78)- (3. 81)  are  satisfied  and 
g » g.  By  comparing  the  results  of  Theorem  3.6  and  3.9,  we  recognize 
that  this  separation  property  does  not  hold  even  when  the  high  gain 
observer  feedback  system  (3.57),  (3.50)  is  0(p)  sensitive  to  disturbance  f. 

We  will  not  consider  the  case  when  g » § since  this  implies 
the  observer  states  are  fed  back  even  before  the  fast  transients  in  the 
estimation  error  are  over.  Such  considerations  are  in  conflict  with  the 
concept  of  observer  feedback. 
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3.5.  Regulation  by  High  Gain  Observer  Feedback 

Under  the  assumption  that  K and  H are  designed  such  that  (3.12), 

(3.13)  and  (3.42),  (3.43),  respectively,  we  now  establish  the  conditions 

for  the  existence  of  a high  gain  observer  feedback  control  that  solves 

the  partial  states  regulation  problem  described  in  Section  3.3,  that 

is,  the  regulation  of  the  variables  z = D x.  Before  we  state  the 

o 

theorem,  it  is  essential  to  decompose  x and  e into  slow  and  fast  variables 
analogous  to  (3.52). 

Lemma  3.1: 

Assume  that  the  matrix  of  (3.59)  is  nonsingular.  Then  the 

°2 

variables  x and  e can  be  expressed  in  terms  of  fast  and  slow  variables 
as 


KB  -a  1KH' 
o 


X 

x’1 

B O" 

1 

0 

= y 

+ 

e 

-i- 

0 H 

6CB 


a 1C  H 

o _ 


(3.90) 


where  $ and  f]  are  the  fast  variables  approximated  by 


+ 0(g-1) 


and  x " and  are  the  slow  variables  approximated  by 


~Ut) 

“•(of 

= gK282t 

fl(  t)_ 

_fK°)_ 

(3.91) 


’x"  (t)" 

__  1 

rx”  (on 

1 

0^-  A2K2 

1 

i 

e1(t) 

- e 

+ 

(A1  ’ M2~\)(c"t) 


+ 0(g_1) 


^ f (x)dT 


(3.92) 


and  the  matrices  A^,  K^,  and  8 2 are  as  defined  in  (D19),  (D23)- 


(D25).  In  (3.90),  the  matrix  4*  is  as  defined  in  (D30) , that  is. 


-K  K 


5C  C 


¥ - 0 


(m+r)  x (2n-m-r), 


(3.93) 


Proof : See  Appendix  D. 


The  theorem  on  regulation  of  z(t)  is  as  follows. 


Theorem  3.10: 


*T  T 

Let  be  the  matrix  whose  columns  span  the  null  space  of  H , 


that  is. 


MjH  = 0. 


(3.94) 


Let  4*^,  4^  be  the  n x (2n-m-r)  submatrices  of  41, 


(3.95) 


T T 

be  the  matrix  whose  columns  span  the  null  space  of  (Bq  + 5B)  , that 


M.  (B  + 6B)  - 0 
61  ° 


(3.96) 


Let  tf  be  a finite  positive  real  number  and  matrix  S be  defined  as  in 


(D5). 


Suppose  that  high  gain  observer  feedback  control  (3.4)  is 


applied  to  system  (3.1).  Then  there  exists  E such  that  for  all  t > t 
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the  regulated  variables  z(t)  * D^x(t)  are  reduced  to  0(g  ^),  that  is 
| z ( c ) | * 0(g  1)  if  and  only  if  system  (3.57),  (3.58)  is  asymptotically 
stable  and 


<A1  - A2K2-1K1|R( 


Mr 


Mjd  - S) 


E ) W(D  4L). 
o — o i 


(3.97) 


Proof : Lemma  3.1  establishes  the  decomposition  of  the  variables  x and  e 

into  fast  and  slow  variables.  Then  the  proof  follows  closely  that  of 
Theorem  3.5,  details  are  omitted. 


Corollarv  3.5: 


Let 


Mr 


M1(I  - S) 


F . 

r 


If  the  pair  (A  - A K.y  \ E ) is  controllable  then  condition 
1 Z Z 1 9 ^ 


C R 


KT  iSCT 


T T 
-K1  C , 

o J, 


(3.98) 
(3.97)  is 


(3.99) 


Proof : If  the  pair  (A  - A?<2  ^Ki»^r)is  controllable,  then 

<Ai  A2^2  ^ m_r.  (3.97)  holds  in  this  case  if  and  only  if 

D V,  = 0 or  equivalently  [D  0]T  = 0.  Since  columns  of  f span 
o 1 o 

M([*  “J  ]),  (3.99)  is  obtained. 

oL  LO 

We  note  from  their  definitions  that  the  matrices  A , A,,  K..,  depend  on 

all  parameter  variations  5A,  5B  and  5c  while  the  matrices  M,.  , S depend 
only  on  5B.  Hence  in  order  to  guarantee  regulation  of  z(t)  for  all 


possible  i$A,  5B  and  5C,  condition  (3.97)  in  Theorem  3.10  is  replaced  by 


r 

J 

l 

l 

i 


l 

i 

l 

l 


l 

l 

l 

l 


<j 

it 
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(3.99) .  The  following  lemma  explores  the  implications  of  condition 

(3.99)  on  the  structure  of  matrix  D . 

o 


Lemma  3.2: 


A necessary  condition  for  (3.99)  to  hold  is  that 


R(DT)  c R([CT  + 6CT]) 
o — o 


Proof : (3.99)  implies  the  existence  of  a matrix  P 


CT  + 6CT 
o 


-r 


such  that 


(3.100) 


This  lemma  shows  that  when  no  restriction  is  imposed  on  parameter  variations, 
then  regulation  of  z(t)  can  be  maintained  only  if  the  regulated  variables  z 
appears  as  linear  combinations  of  the  measured  variables  y.  We  note  that 
so  far  we  assume  that  high  gain  feedback  is  made  of  an  observer  whose  sensi- 
tivity with  respect  to  parameter  variations  and  disturbances  is  not  reduced 
to  0(g  ■*")  by  the  observer  high  gain  loops.  By  assuming  that  the  distur- 
bances and  parameter  variations  are  "rejected"  by  the  high  gain  observer, 
we  then  have  the  following  conditions  on  regulation  of  z(t). 


Theorem  3.11: 

Suppose  that  the  matrices  5A,  6B,  6C  and  satisfy  (3. 47)-(3. 50) 
and  the  high  gain  observer  has  the  properties  of  Theorem  3.3.  Let  the 
condition  (3.82)  be  satisfied  and  tf  be  a finite  positive  constant.  Then 
there  exists  t such  that  for  all  t > t , z(t)  is  reduced  to  an  0(g  ^) 
quantity,  that  is,  |z(t)|  = 0(g  ^)  if  and  only  if 


i 


• . 1 ' ■ 


l 
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<A11  - A12K21kiIR(M51Eo)>Q  W(DoQ) 


(3.101) 


and 


ReX(A u - A12K"1Ki)  < 0 


(3.102) 


where  Q is  the  matrix  in  (3.53),  >1^  as  in  (3.96)  and  A^ , are  as 


defined  in  (D12)-(D16). 


Proof:  If  6A,  5B,  5C  and  satisfy  conditions  (3.47)— (3. 50 ) let 


>1^  = 0 in  the  transformation  matrix  M of  (Dl).  Then  = e^  and  c0  *»  e., 

and  K.  , £ . , <5C.  and  5F..  are  null  matrices.  Since  SC  = 0,  we  can  take 
i i i ii 


H1  in  (3.93)  to  be 


Q N 
0 N 


(3.103) 


where  N is  as  in  (3.19),  that  is,  it  satisfied  C N = 0*  Then  (3.90) 


becomes 


x = Qx^’  + Ne^  + B (KB  )_1  $ + KH(CoH)-1n 


(3.104) 


If  (3.82)  holds,  then 


(A  - A K^K  )C  -1 

xj  (t)  = e ‘ x”  (0)  + Ar(K;  K.,C°  C°  - 1C"  K1)e1(t) 


,-l, 


(A  - A K XK  )(t  - T) 
e E1f(T)dT  + 0(g  L) 


(3.105) 


(A°1  - Z°  C° 

e: (t)  = e ' ~ ex(0)  + 0(g  l) 


(3.106) 


and  3(t),  n(t)  are  given  in  (D26).  The  remaining  of  the  proof  is 
analogous  to  that  of  Theorem  3.5. 


.■^■1  - - ■ _ . • 2; ..  . 
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Corollary  3.6: 


If  Che  pair  (A,.  - A K K.  , M.  E ) is  controllable,  then 
r 11  12  2 1 o 


(3.101)  is  replaced  by 


T T 

R(D*)  C R(KX) 


(3.107) 


which  is  the  same  condition  as  (3.51). 


Condition  (3.107) implies  that  l £ m,  that  is,  the  number  of  regulated 
variables  should  not  be  more  than  the  number  of  control  variables. 

The  coincidence  of  condition  (3.107)  and  (3.51)  agrees  with  intuition. 
Since  the  high  gain  observer  rejects  the  parameter  variations  and 
disturbances  and  thus  provides  estimated  states  x which  tracks  x 
asymptotically,  the  condition  for  regulation  of  z(t)  should  be  the 
same  whether  high  gain  state  feedback  or  high  gain  observer 
feedback  is  employed.  However,  we  note  that  if  high  gain  state 
feedback  is  used,  the  transient  responses  are  faster,  that  is,  z(t)  is 
reduced  to  0(g  ^)  after  the  time  tg  which  is  0(g  ^) . Since  the  matrices 
A , K and  M,  depends  on  parameter  variations  <5A,  SB,  in  order  to 

lj  1 

guarantee  regulation  of  z(t)  for  all  possible  parameter  variations,  »/e 
will  focus  on  condition  (3.107)  when  we  consider  the  synthesis  of  matrix 
K in  the  next  section. 

If  C •=  D , that  is,  the  measured  variables  are  also  the 
o o 

regulated  variables,  and  m = r = i,  then  the  following  corollary  of 
Theorem  3.11  is  obtained. 


Corollary  3.7: 

Let  C = D and  m * r = i.  Suppose  it  is  required  to  regulate 
o o 

z(t)  for  all  possible  parameter  variations  5A,  SB  and  disturbance  matrix 
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E^,  then  under  the  same  assumptions  of  Theorem  3.11,  the  conditions 
(3.101),  (3.102)  are  replaced  by 


T T 

R(CX)  = ROC1) 
o 


(3.108) 


and  the  condition  that  the  transmission  zeros  of  the  triple  (C  ,A  + 5A, 

oo 

B^  + 6B)  all  lie  in  the  open  left  half  plane. 

Proof:  IfC  = D , ra  = r = l,  then  (3.107)  is  (3.108).  Consequently, 
o o 

X(A,,  - A,.K„^K, ) are  the  transmission  zeros  of  the  triple  (C  , A + A, 

11  12  l 1 o o 

B + SB)  from  the  results  of  Section  2.2. 
o 

We  see  that  when  the  number  of  control  variables  is  the  same 

as  the  number  of  measured  variables  and  the  measured  variables  are  the 

variables  to  be  regulsted,  then  in  spite  of  the  use  of  a "disturbance 

rejecting"  high  gain  observer  in  the  high  gain  feedback  loops,  the 

finite  eigenvalues  of  system  (3.51),  (3.58)  are  determined  bv  the  system 

parameters  A + SA,  B + SB  and  C . In  fact,  Theorem  2.2  shows  that  these 
r o o o 

finite  eigenvalues  coincide  with  the  finite  eigenvalues  of  the  system 
(3.1)  when  high  gain  output  feedback  u = gC^x  is  applied.  If  m > we 
recognize  that  there  exists  some  degrees  of  freedom  which  allows  the 
placement  of  the  finite  eigenvalues.  These  aspects  of  the  synthesis 
of  the  matrices  K and  H are  the  subject  of  the  next  section. 


3 . 6 Synthesis  Problems 

So  far  in  this  chapter,  we  have  demonstrated  that  high  gain  feed- 
back is  capable  of  reducing  the  sensitivities  of  control  systems  with 
respect  to  parameter  variations  and  disturbances  to  as  desired  providing 
that  the  ranges  of  SA,  SB  and  lie  in  the  proper  ranges.  The  range 
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conditions  that  involve  the  range  of  B and  C are  dictated  by  the 

o o 

input  output  structures  of  the  open  loop  system.  The  ones  that  involve 
range  of  K and  H can  be  manipulated  since  matrices  K and  H are  design 
parameters.  However,  we  recall  that  these  matrices  also  determine  the 
dynamical  characteristics  of  the  high  gain  feedback  system.  Hence, 
we  are  faced  with  synthesis  problems  with  two  objectives,  that  is  the 
matrices  K and  H are  designed  such  that  sensitivity  reduction  by  high 
gain  feedback  is  achievable  while  the  high  gain  feedback  system  possesses 
desirable  characteristics.  We  recognize  that  in  multivariable  systems, 
the  solution  to  some  feedback  gains  design  problems  are  nonunique  and 
our  success  in  achieving  the  design  goal  depends  on  the  degrees  of 
freedom  available  to  satisfy  both  objectives. 

The  range  conditions  that  we  want  to  manipulate  are  as  follows 


R(5A)  E R(H),  R(Sb)  c R(H),R(E  ) E R(H) 

o 


R(dJ)  C r(kt). 


(3.47)— (3.49) 
(3.107) 


For  design  purposes,  there  is  no  loss  of  generality  to 
assume  that  R(5A)  C R(E),  R(5B)  c R(E)  and  R(Eq)  E R(E) . Thus,  in  place 
of  (3. 4 7)  — ( 3 . 49) , we  consider  the  condition 


R(E)  <=  R(H) . 


(3.109) 


The  design  problems  regarding  matrices  K and  H can  be  defined 
formally  as  Problem  K and  Problem  H. 


i 
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Problem  K: 

Design  matrix  K such  that 

U)  R(DT)=R(KT)  (3.107) 

o 

(ii)  ReA(A°x  - A°2Ks)  < 0 (3.12) 

(ill)  ReA(KB  ) < 0.  (3.13) 

o 

Problem  H: 

Design  matrix  H such  that 


(i) 

R(E)  c R(H) 

A /V  ”1 

(3.109) 

(ii) 

ReX(A°x  - A°2C°  C°)  < 0 

(3.42) 

(iii) 

Re\(C  H)  < 0. 
o 

(3.43) 

Although  the  specifications  (ii)  and  (iii)  demand  only  asymptotic  stability 
of  the  fast  and  slow  modes  of  the  high  gain  feedback  systems,  feasibil- 
ities for  slow  and  fast  pole  placement  will  also  be  investigated.  We 
have  not  taken  into  account  the  effects  of  parameter  variations  on 
stability  in  these  specifications  since  these  are  not  the  problems  of 
concern  for  the  design  of  K and  H.  We  recognize  that,  from  duality,  the 
design  procedures  to  solve  problems  K and  H are  the  same.  Therefore,  we 
derive  solvability  conditions  only  for  Problem  K.  First,  we  state  the 
constraints  of  requirements  (ii)  and  (iii)  on  the  structure  of  matrix  K. 


Lemma  3.3 


The  matrix  KB  is  nonsingular  if  and  only  if 


Proof 


such  that  R(B  q) C N(K).  But  this  implies  rank  (KB  ) < m which  is  a con- 


tradiction. To  prove  sufficiency,  (3.110)  implies  there  exists  no  nonzero 


vector  q such  that  KB  q = 0.  It  follows  that  KB  is  nonsingular 


Corollarv  3.8 


Condition  (3.110)  is  equivalent  to 


Lemma  3.3  and  its  Corollary  shows  that  (iii)  restricts  the  range  of  K 


This  observation  leads  to  the  first 


solvability  condition  of  Problem  K 


Lemma  3.4:  A necessary  condition  for  the  solution  of  Problem  K to 


exist  is  that 


Proof 


hence  (3.112)  is  necessary  to  satisfy  both  (i)  and  (iii) 


Thus,  we  assume  that  matrix  D satisfies  (3.112).  This  assump 


M x where 


tion  allows  us  to  introduce  the  transformation  x 
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M = I I (3.113) 

is  as  defined  in  (2.10)  and  the  m x n matrix  is  chosen  such  that 

(3.114) 


Lemma  3.5: 


Let 


and 


R(dt)c  R(m'T). 
O ” z 


KM,_1  = [K'k'  K'] 

t s x 


U)'  = KM'-1x'. 
Then  the  solution  of  the  system 


(3.115) 

(3.116) 


x = (Aq  + gBoK)x 


can  be  approximated  by 

(A’°  - A'^K' )t 

• x^(t)  = e 11  12  x|(0)  • 

gK’B°  . 

u)’  (t)  = e V (0)  + 0(g~l). 

where 


0(g  l) 


M'A  M' 
o 


,"1 


A’° 

ii 

12 

A’° 

A’° 

_ 21 

12 

(3.117) 


(3.118) 

(3.119) 


(3.120) 


and  B0  is  as  defined  in  (2.9). 
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Proof : The  proof  is  analogous  to  that  of  Theorem  2.1,  details  are 

omitted. 

Clearly,  (ii)  and  (iii)  are  equivalent  to 

(ii')  ReX(A’°  - A'°K’)  < 0 
11  12  s 

(iii')  ReX(K^B°)  < 0 

respectively.  The  motivation  behind  the  use  of  the  transformation  M'  is 
clear  from  the  following  lemma. 


Lemma  3.6: 

Let  T be  an  £ x a matrix  satisfying 

D*  = m’2  rT  (3.121) 

then 


D M'_1 
o 


^£x(n-m) 


n 


(3.122) 


Proof:  Direct  computation  gives  (3.122). 

We  thus  have  the  following  solvability  condition  for  Problem  K. 


Theorem  3.12: 

Let  II  be  an  m x (m-£)  matrix  whose  columns  span  the  null  space 
of  f defined  in  (3.121),  that  is, 


rn  = o 


(3.123) 
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Suppose  that  satisfies  (3.112),  then  a necessary  and  sufficient  condi- 
tion for  the  existence  of  a matrix  K which  fulfills  the  condition 
(3.107),  ( i i ' ) and  (iii')  is  that  the  pair  (A^°,  A^°JI)  be  stabilizable. 

Proof : Condition  (3.107)  holds  if  and  only  if  there  exists  an  m x £, 
matrix  P such  that 

K’TKlTP  = 0 (3.124) 

S I 

and 

K^TP  *=  rT.  (3.125) 

(iii')  means  that  K'^  ^ exists,  hence  (3.124)  becomes 

HC'  = 0.  (3.126) 

s 

Thus,  (3.107)  is  equivalent  to  the  matrix  has  the  form 

K'  * UK”  (3.127) 

s s 

where  K^'  is  a matrix  of  appropriate  dimensions.  Substituting  (3.127) 
into  (ii' ),  it  becomes 

ReA(A'°  - A'°IIk"  ) < 0.  (3.128) 

11  12  s 

Clearly,  there  exists  matrix  k”  such  that  (3.128)  is  satisfied  if  and 
only  if  the  pair  (A^,  A^II)  is  stabilizable.  Since  (iii')  implies 
is  nonsingular,  there  exists  k'^  such  that  (iii'  ) is  fulfilled.  This 
concludes  the  proof  of  this  theorem. 

The  condition  for  the  existence  of  matrix  K that  places  the  eigenvalues 
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1 


subjected  to  the  structural  constraint  (3.107)  is  evident  from  the 
proof.  We  have  the  following  theorem. 

Theorem  3.13K: 

Let  pi , i = 1, . . . ,n-m  and  gq^ , j = 1, ... ,m  be  the  desired 
locations  for  the  slow  and  fast  eigenvalues  of  the  system  (3.117).  Then 
under  the  same  assumptions  of  Theorem  3.12,  the  necessary  and  sufficient 
condition  for  the  existence  of  a matrix  K such  that  condition  (3.107) 
is  fulfilled  and 


XK  - 


(3.129) 


*(KjB°)  = q. 


(3.130) 


is  that  the  pair  (A^°,  A^JI)  be  controllable. 


If  H = m,  then  clearly  T = 1^  and  II 


= 0.  The  following  corollary 


is  obtained. 


Corollarv  3.9K: 


If  i 8 m,  then  the  necessary  and  sufficient  condition  in 


Theorem  3.12  becomes 


ReX(A,11)  < 0 


(3.131) 


that  is,  the  transmission  zeros  of  the  triple  (D  , A , B ) lie  in  the 

r o o o 

open  left  half  complex  plane. 


Proof:  If  Z = m,  then  F = I and  II  = 0 and  (3.128)  becomes  (3.131). 
m 

From  the  results  of  Section  2.2,  X(A|°)  are  the  transmission  zeros  of 

the  triple  (D  , A , B ). 
r o o o 

In  Theorem  3.12  and  3.13K,  the  solvability  of  Problem  K is 
stated  in  terms  of  II  which  depends  on  the  transformation  matrix  M^.  We 
now  show  that  the  controllability  of  the  pair  (A|°,A^°n)  dependent 
only  on  system  matrices  Bq  and  Dq. 

Theorem  3.14K: 

Let  Q be  an  n * (m-£)  matrix  with  rank  = m-Z  such  that 

R(f2)  H [R(D*XjN(B*)]  = 0 (3.132) 

Then  the  controllability  of  the  pair  (A^°,  Aj°n)  e<lu:'-valent:  t0 
controllability  of  the  pair  (M^AqS^ .M^A^fi)  where  the  inverse  of  M?  of 
(3.113)  is  defined 

M'"1  = [Sx  S^\  (3.133) 

ra 


Proof : By  definition  (3.120) 


A11  = MlAoSl  ’ A12  “ MlAoS2 


To  prove  the  equivalence  statement,  we  need  only  to  show  ^(M^A^S^H)  = 


R(M,A  f2) . Let  M'  be  chosen  such  that 
1 o 2. 


M 


(3.134) 


where  G is  a (m-2.)  * n matrix  such  that  rank  Ml  = m. 
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Denote 


S2  = [S1,  s22] 


l m-£ 


Since  MIS.  = I , 
l l ra 


°oS2  = °rM) 


GS2  = Vi. 


For  the  choice  of  Ml,  in  (3.134), 


s2n  - S9 


It  follows  that 


R(SJI)  c W(D  ) 

£■  O 


R(s2rr)  n w(g)  = 0 


(3.135) 


(3.136) 


(3.137) 


(3.138) 


(3.139) 


(3.140) 


Since  R(GT) , R(D*)  and  N(B*)  span  Rn,  (3.139),  (3.140)  leads  to 
R(S,II)  n [R(DT)UW(BT)]  = 0 

L O O 


This  concludes  the  proof. 


(3.141) 


We  now  summarize  the  solvability  of  Problem  H in  the  following 


Theorem. 


Theorem  3.1311: 


Assume  that 


R(E)  n W(C  ) = 0. 
o 


(3.142) 


Let  M be  as  n ^ n matrix  defined  bv 


(3.143) 


where  H0  is  an  m * n matrix 


T 0 T ° 

- T-  -1  11  12 

MA  M = 

0 7 ° T o 

A,,l  "22 


Let  r be  an  r * p matrix  satisfying 


(3.144) 


(3.145) 


and  let T bo  an  r x (r-p)  matrix  whose  columns  span  the  null  space 


of  F,  that  is. 


m = o 


(3.146) 


Then  a necessary  and  sufficient  condition  for  the  existence  of  a matrix 

A A “ 1 

H such  that  (3.109)  is  satisfied  and  \(A?.  - A°  C°  0°)  and  \ (C  H)  are 

ll  12-1  o 

positioned  at  desired  locations  is  that  the  pair  (A^°,  A ^JI)  be  con- 


trollable. 


I t T 

Proof:  Bv  let t inn  K *-*  H , B '-*•  C , D *-*•  E , A •*-*•  A1  and  applvins 

■ oooooo  • ° 


Theorem  3.13K,  this  theorem  is  proved. 


The  dual  theorem  of  Theorem  3.14K  is  stated  below. 


Theorem  3.14H: 

Let  £2  be  an  n * (r-p)  matrix  with  rank  £2  = r-p 

such  that 


R(ft)  n [R(E  ) U N(C  )]  = 0. 


(3.147) 


Then  the  controllability  of  the  pair  (A^°,  A^°  i-s  e<?uiva^-ent  to  c^le 

— T—  — T—  — 

controllability  of  the  pair  (M,  A S.,M  A £2)  where  the  inverse  of  M of 

1 o I J.  o 

(3.143)  is  defined 


M = [S  S2J. 


(3.148) 


If  p = r,  then  T = I and  II  = 0.  The  dual  corollary  of  Corollary  3.9K  is 
as  follows. 


Corollary  3.9H: 

• If  p = r,  then  the  necessary  and  sufficient  condition  in 
Theorem  3.13H  becomes 

ReX(A  °)  < 0 (3.149) 


that  is,  the  transmission  zeros  of  the  triple  (Cq,  A^,  E)  all  lie  in 
the  open  left  half  complex  plane. 
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CHAPTER  4 

VARIABLE  STRUCTURE  FEEDBACK  SYSTEMS 
WITH  SLIDING  MODE 


Another  class  of  feedback  controls  capable  of  reducing 
parameter  sensitivities  and  rejecting  disturbances  is  the  so-called 
variable  structure  control.  Basically,  it  is  a feedback  control  dis- 
continuous on  a switching  surface  s (x)  = 0, 


ui(x) 


i£(x)  s.(x)  > 0 
ui(x)  S (x)  < 0 


(4.1) 


where  u^  and  s^  denote  the  i component  of  the  m-vectors  u and  s,  re- 
spectively. The  class  of  control  systems  with  discontinuous  feedback 
control  (4.1)  is  called  variable  structure  systems  (VSS)  and  has  been 
developed  in  the  USSR  in  the  last  fifteen  years  [14,15,16,17,18].  It  has 
found  applications  in  control  of  a wide  range  of  processes  in  steel, 
power,  chemical  and  aerospace  industries.  For  a recent  survey, 
see  [19].  The  salient  feature  of  VSS  is  that  the  so-called  sliding 
mode  occurs  on  the  switching  surface  s(x)  = 0.  While  in  sliding  mode, 
the  feedback  system  may  be  insensitive  to  parameter  variations  and 
disturbances . 

The  organization  of  this  chapter  is  as  follows.  We 
first  review  the  properties  of  sliding  mode  in  VSS.  Then  we  establish 
the  relationship  between  high  gain  feedback  systems  and  VSS  with  sliding 
mode.  We  reveal  the  insensitivity  property  of  sliding  mode  by  drawing 
an  analogy  with  high  gain  feedback  systems  which  we  have  investigated  in 
the  previous  chapters.  First,  the  case  when  all  the  states  are  accessible 
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is  considered.  We  then  examine  the  behavior  of  VSS  with  Luenburger 
observer  feedback.  A variable  structure  observer  which  is  closely 
related  to  the  high  gain  observer  of  Section  3.2  is  developed.  Finally, 
we  look  into  VSS  with  variable  structure  observer  feedback. 


4.1.  Properties  of  Sliding  Mode 

To  illustrate  some  of  the  fundamental  concepts  in  the  theory 
of  variable  structure  systems  with  sliding  mode,  we  consider  a second 
order  linear  time-invariant  system, 


X1  " x2 


X2  - ax2  - bu. 


a,  b > 0 


let  the  variable  structure  control  law  be 
'ax.^  , x^s  > 0 

l-ou^  , x^s  < 0 


where  the  switching  line  is 


j 

(4.2) 


(4.3) 


s * cx^  + x.,  = 0,  c > 0.  (4.4) 

The  phase  plane  trajectories  of  the  closed  loop  system  (4.2), (4.3)  are 
illustrated  in  Figure  4.1.  They  are  a combination  of  phase  trajectories 
of  two  linear  feedback  systems  with  u = ctx^  and  u = -ctx^.  One  of  these 
systems  is  marginally  stable  and  the  other  is  unstable.  By  selecting  a 
such  that  s and  s have  opposite  signs  in  the  neighborhood  of  s = 0,  that 
is. 
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lim  s > 0 and  lim  s < 0,  (4.5) 

s-*— 0 s-M-0 


the  phase  trajectories  will  continue  along  the  switching  line  (4.4).  We 
call  the  motion  along  the  switching  line  sliding  mode.  We  note  that 
sliding  mode  is  not  a part  of  trajectories  which  belong  to  either  of 
the  two  feedback  systems.  When  system  (4.2),  (4.3)  is  in  sliding  mode, 
its  trajectories  satisfy  s = cx^(t)  + X£(t)  = 0,  that  is,  x^(t)  is 
described  by  the  first  order  system 


- cx^. 


(4.6) 


We  note  that  there  is  a reduction  in  system  order  when  sliding  mode 
occurs.  The  reduced  order  system  is  obtained  by  the  equivalent 
control  method  proposed  by  Utkin  [38].  According  to  his  method,  u is 
solved  from  the  algebraic  equation 


CX2  + ax2  - bu 


(4.7) 


The  unique  solution  of  (4.7)  called  the  equivalent  control  ug^  is  then 
substituted  for  u in  (4.2)  and  (4.6)  is  obtained  by  restricting  x^,  x,, 
to  satisfy  s = 0.  We  call  (4.6)  the  equivalent  control  system.  It  is 
dependent  only  on  the  parameter  c of  the  switching  plane  s = 0 of 
(4.4).  This  insensitivity  property  in  sliding  mode  depends  on  the 
system  structure  (4.2).  In  general,  the  equivalent  control  system  which 
describes  the  dynamical  behavior  of  the  VSS  in  sliding  mode  depends  both 
on  system  parameters  and  switching  plane  parameters.  The  design  of  VSS 
with  sliding  mode  thus  involves  the  Synthesis  of  the  switching  planes 
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such  that  the  equivalent  control  system  is  asymptotically  stable.  In  this 
example,  by  choosing  c to  be  positive  guarantees  the  asymptotic  stability 
of  (4.6)  and  in  fact  -c  is  the  eigenvalue  of  (4.6). 

Another  important  consideration  in  the  design  of  VSS  in  the 
condition  for  the  system  trajectory  to  reach  the  switching  planes. 

For  the  linear  system  (4.2),  it  happens  that  by  choosing  the  vari- 
able structure  control  as  in  (4.3),  all  its  trajectories  reach  the 
switching  line.  For  illustration  purposes,  we  consider  the  system 
(4.2)  subjected  to  parameter  variations  and  disturbance,  that  is. 


u = 


(4.10) 


a9sgn(b) , x2s  > 0 


-oc,sgn(b),  x9s  < 0, 


(4.11) 


and  a_  > 0, 


awsgn(b),  ws  > 0 


-awsgn(b) , ws  < 0 


(4.12) 


and  a >0,  then  condition  (4.9)  is  satisfied  if 
w 


a-  > 

2 m„ 


a > — 
w m- 


(4.13) 


We  note  that  the  equivalent  control  system  for  the  VSS  (4.8),  (4.10)  is 


(4.6)  which  is  the  same  for  the  VSS  (4.2),  (4.3).  This  clearly  exhibits 
the  insensitivity  with  respect  to  parameter  variations  and  the  disturbance 


rejection  capability  of  VSS  with  sliding  mode. 

If  the  disturbance  w(t)  is  not  measurable,  then  by  letting  the 
variable  structure  control  be 


u - Y2x2  + H'r 


(4.14) 


where  V,,  is  as  in  (4.11)  and 


arsgn(b),  s > 0 


-arsgn(b) , s < 0 


(4.15) 


and  a >0,  the  condition  (4.9)  is  satisfied  if 
r 


4 
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a. 


c + m2 


(4.16) 


The  difference  between  the  variable  structure  controls  (4.10)  and  (4.14) 
is  the  component  associated  with  w(t).  When  the  disturbance  w(t)  is  accessi- 
ble. then  it  is  'fed  forward"  with  a switching  gain  as  in  (4.10).  If  w(t) 
cannot  be  measured,  then  a relay  function  is  adopted  as  in  (4.14).  This 
determination  of  the  variable  structure  control  law  according  to  the  ac- 
cessibility of  the  disturbances  carries  over  in  the  multivariable  case. 

We  now  summarize  some  of  the  results  in  [19,37,38]  which  we  need 
later.  Consider  the  linear  multivariable  system 

• 

x = Ax  + Bu  + E^f  (4.17) 

where  as  in  (3.1),  the  state  x t 8°,  the  control  u £ R™  and  the  disturb- 

p 

ance  f £ R . Assuming  that  f is  not  measurable,  we  let  the  variable 
structure  control  law  be  of  the  form 

u(x)  = 4»°x  + i>b  (4.18) 

where 


Y = 

Tdiag( 

sgn  x , . . . 

» sgn  X ) 
n 

(4. 

.19) 

H 

a 

(<[>,.••  • 

,*> 

m 

(4. 

, 20) 

1 ^ 

. si(x) 

> 0 

♦i- 

, i * 1, . . . ,m,b 

(4. 

21) 

K 

» si(x) 

< 0 
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and  s^(x)  = 0 are  the  switching  planes.  We  define  an  m-vector  s(x) 
whose  iC^  component  is  denoted  as  s^(x)  by 


s(x)  = [s,(x),...,s  (x) ] = Kx  = 0. 

i m 


(4.22) 


In  (4.20),  (4.21),  and  i = l,...,m,  are  n-vectors  and  and  r^ 
are  m-vectors.  We  note  that  u^,  the  iC^  component  of  u,  is  discontinuous 
on  s^(x)  = 0 and  is  a linear  function  in  x if  x is  not  on  any  of  the 
switching  planes. 


Lemma  4.1: 


Assuming  that  the  matrix  KB  is  nonsingular  and  sliding  mode 


exists  on  the  intersections  of  the  switching  planes  s(x)  = 0 (4.22), 


then  an  equivalent  control  ug^  exists  and  is  unique, 

u = (KB)_1K[Ax  + E f], 
eq  o 


(4.23) 


Furthermore,  when  sliding  mode  occurs,  the  VSS  (4.17),  (4.18)  is  governed 


x = (I  - B(KB)  K) [Ax  + E f ] 

o 


Kx  = 0. 


(4.24) 


(4.25) 


Proof:  Following  the  equivalent  control  method  [38],  (4.23)  is  the 


unique  solution  of  u in  the  algebraic  equation 


s = KAx  + KBu  + KE  f = 0. 

o 


(4.26) 


Substituting  (4.23)  into  (4.17),  (4.24)  is  obtained  and  the  constraint 
(4.25)  is  the  sliding  mode  condition  s(x)  = 0. 


We  call  (4.24),  (4.25),  the  equations  of  sliding  mode.  In 
the  design  of  VSS  with  sliding  mode,  we  consider  the  synthesis  of  the 
switching  planes,  that  is,  the  choice  of  a matrix  K such  that  the 
dynamical  behavior  of  the  VSS  in  sliding  mode  is  desirable.  Some 
synthesis  procedures  are  developed  in  [15,17,37].  As  it  is  illustrated 
in  the  example,  the  second  design  problem  (the  so-called  reaching 
problem)  in  VSS  is  to  guarantee  that  the  trajectories  move  towards  the 
intersections  of  the  switching  planes  and  that  sliding  mode  exists.  In 
[19,37],  the  reaching  conditions  and  the  existence  conditions  are  derived 
from  the  conditions  for  the  asymptotic  stability  in  the  large  and  in  the 
small,  respectively,  of  the  subspace  s(x)  = 0 for  the  discontinuous  feed- 
back system  (4.17),  (418).  To  derive  these  conditions,  it  is  necessary 
to  construct  Lyapunov  functions  for  the  discontinuous  system 

s = KAx  + KBu  4-  KE  f (4.27) 

o 

Since  there  is  no  standard  method  to  construct  Lyapunov  functions  for 
arbitrary  nonlinear  systems,  the  existence  and  reaching  conditions  can  be 
derived  only  when  appropriate  Lyapnov  functions  are  found.  Two  cases 
when  the  Lyapunov  functions  exist  are  summarized  in  the  next  lemma  [37]. 

Lemma  4.2; 

If  the  matrices  K and  B are  such  that  either 
(i)  there  exist  a positive  definite  symmetric  matrix  W and  a 
matrix  L which  is  diagonally  dominant  such  that 
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or 


(ii)  KB  is  symmetric, 

then  there  exists  Lyapunov  functions  which  are  quadratic  forms  of  s. 

Alternatively,  the  existence  and  reaching  conditions  can  be 
found  by  utilizing  the  conditions  for  VSS  with  scalar  control  where 
there  is  only  one  switching  plane  a = 0.  If 

oa  < 0,  (4.29) 

that  is,  a and  o have  opposite  signs,  then  the  trajectories  move  toward 
0=0  and  sliding  mode  exists.  Two  cases  when  it  is  possible  to 
decouple  the  multivariable  problem  into  m scalar  problems  is  given  in 
the  next  lemma  [37]. 

Lemma  4.3: 

If  the  matrices  K and  B are  such  that  KB  is  diagonally  dom- 
inant then  the  conditions  on  r^  of  (4.21)  such  that  the  trajectories 
reach  the  switching  plane  s^^  = 0 and  sliding  mode  exists  on  s^,  = 0 can 
be  derived  from 

sisi  < 0 * 1 “ 1. •••.“•  (4.30) 

From  the  last  two  lemmas,  we  see  that  in  order  to  derive  the  existence 
and  reaching  conditions,  it  is  necessary  that  the  matrix  KB  is  of  specific 
structures.  The  following  lemma  [37]  enlarges  the  class  of  matrix  KB. 
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Lemma  4.4: 

The  equations  of  sliding  mode  (4.24),  (4.25)  are  invariant 
with  respect  to  switching  plane  transformations 

s = H s (4.31) 

s 

and  input  transformations 

u*  = H u (4.32) 

u 

where  H and  H are  m x m nonsingular  matrices, 
s u 

If  the  matrix  B is  known,  then  there  exists  matrix  (Hs) 

such  that  with  respect  to  the  new  inputs  (new  switching  planes),  the  new 

KBq  belongs  to  either  cases  in  Lemma  4.2  and  4.3. 

We  note  that,  however,  in  the  presence  of  parameter  variations 

in  B,  there  exists  no  transformation  H or  H that  reduces  to  the  cases 

s u 

of  Lemma  4.2,  4.3.  On  the  contrary,  the  hierarchy  of  controls  method 
which  is'also  given  in  [19,37],  does  not  require  KB  to  have  any  particular 
structure.  Hence,  this  method  is  suitable  for  our  purposes  since  all 
the  VSS  we  considered  later  have  parameter  variations.  To  introduce 
this  method,  we  use  the  following  definitions. 

Definition  4.1: 

Let  u^  be  discontinuous  on  s^  = 0.  If  u^  is  designed  so  that 
sliding  mode  occurs  on  s^  = 0 and  assuming  s^  = 0,  u?  in  designed  so 
that  sliding  mode  occurs  on  s9  = 0,  then  the  hierarchy  of  controls  is 


U1  u2’ 


Definition  4.2: 


If  the  hierarchy  of  controls  is  u^  + then  the  hierarchy 
of  switching  planes  is  s^  -*■  S2,  that  is,  sliding  mode  first  occurs  on 
s^  ■ 0 and  then  on  the  intersections  of  s^  = 0 and  S2  =*  0. 

The  hierarchy  of  controls  method  consists  of  the  following  steps. 
Step  1.  Suppose  a hierarchy  of  controls  is  given, 

12  m 

Step  2.  We  begin  with  the  bottom  switching  plane  s = 0 in 

m 

the  hierarchy  by  letting  i = m. 

Step  3.  Suppose  sliding  mode  occurs  on  the  first  i-1 

+ 

switching  planes,  that  is,  s^  ■=  0,  j = l,...,i-l.  Solve  for  the 

equivalent  control  u*  ^ of  the  variable  structure  controls 

eq 

i-1  _ , , T r . - . i+1  _ . 

u = (u^,...,u^  ^)  as  a function  of  u^  and  u = (u^+^,..., 

T • 

) from  the  algebraic  equations  s =0,  j * l,...,i-l. 
i+1  m j 

Since  (4.17)  is  a linear  svstem,  u^  ^ is  linear  in  x,  u*+\  f and  u., 

eq  i 

that  is,' 


P.  .x  + Q.  .u  + T.  .f  + ct.u. 
i-1  l-l  i-l  i 1 


(4.33) 


where  matrices  P^_^*  an<^  Ti-i  an<*  scalar  depend  on  the  first 

i-1  rows  of  matrices  KA,  KB  and  KE  . We  note  that  ui+1  is  known 

o 

since  r^,r^,  j=i+l,...,ra  have  been  determined  previously. 

Step  4.  For  the  i switching  plane  s^  = 0,  solve  for  r^  such 


SiSi  < °* 


(4.34) 


For  i = 1,  we  do  not  assume  sliding  mode  occurs  on  any  switch- 


ing plane. 
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The  resulting  conditions  on  are 

Qi  trilxl+  + rbl  ~ " min  [pIx  + cif  + <bf  + a*  ^ 0 (4.35) 

i+1  , 
u ,f 

where  r*  denotes  the  iC^  component  of  r^  and  the  vectors  p^,  t^  and 
depend  on  the  matrices  ^ and  scalar  a^,  respectively 

and  all  of  them  depend  on  the  i*^  row  of  the  matrices  KA,  KB  and  KEq.  In 

(4.35) ,  we  assume  that  the  bounds  on  the  components  of  the  disturbance 
f are  known.  Since  u^"+^  is  of  the  form  (4.18),  the  minimizati  . in 

(4.35)  can  be  performed  component  wise,  that  is, 


* u 

Vij  - - min  [(pij  t £ qikrkj)sgnxj]’j  = 1****»n 


(4.36) 


r . ..... ,r  . k =i+l 
1+1.J  m,J 


* i < 

Vb  - 


, i+1  m 

t,rb, . . . ,rb 


k=i+l 


(4.37) 


where  p and  r„  are  the  j component  of  p_^  and  r^,  respectively. 

Step  5.  Let  i = i-1.  If  i > 0,  go  to  Step  3,  else  stop. 
Because  of  the  inequality  conditions  (4.36),  (4.37)  in  the  above  proced- 
ures, the  hierarchy  of  controls  method  is  suitable  for  systems  with 
parameter  variations.  In  this  case,  the  maxima  of  the  arguments  on  the 
right  hand  sides  of  (4.36),  (4.37)  with  respect  to  the  uncertainty 
parameters  is  minimized.  It  is  sufficient  to  know  only  the  bounds 


T T 

Let  x = (x1,...,xn)  , then  |x|  = ( | | , . . . , | | ) . 


of  the  variations  in  order  to  determine  the  switching  gains  r^. 

For  the  hierarchy  of  controls  specified  in  Step  1,  the  solution  of  r^ 
assume  the  worst  case  for  r , r^,  j *=  i+l,...,m  and  in  particular, 
r^  assures  sliding  mode  on  s^  = 0 assuming  that  sliding  mode  does  not 
occur  on  Sj  » 0,  j = 2,...,m.  We  will  employ  this  method  exclusively  for 
the  design  of  VSS  with  sliding  mode  in  this  study. 

Thus  far,  we  have  considered  the  ideal  behavior  of  VSS  with 
sliding  mode,  that  is,  we  assume  that  the  switching  of  the  variable 
structure  control  when  the  state  x is  on  the  switching  planes  is  infi- 
nitely fast  and  consequently,  the  trajectories  lie  in  the  intersection 
of  the  switching  planes  when  sliding  mode  occurs.  In  reality,  due  to 
switching  delays,  neglected  time  constants,  hysteresis,  etc.,  the  control 
is  switched  at  a finite  frequency  and  the  corresponding  trajectories  are 
in  some  vicinity  of  the  switching  planes.  We  refer  to  these  trajectories 
as  nonideal  sliding  mode.  It  is  shown  in  [38]  that  nonideal  sliding  mode 
is  "close"  to  the  ideal  sliding  mode  in  the  sense  described  in  the  next 
theorem  [38]. 

Theorem  4. 1 : 

Let  the  behavior  of  the  VSS  (4.17),  (4.18)  with  nonideal  switch- 
ing devices  be  described  by  the  system 

x = Ax  + Bu  + E f (4.38) 

o 

where  the  nonidealities  are  taken  into  account  by  the  control  u.  Suppose 
that  u is  such  that  x(t),  the  motion  of  (4.38),  exists  and  is  in  a 
A-neighborhood  of  s(x)  ■ 0,  that  is, 


5 


I 


84 


| s (x ( t ) ) I s A,  t € (4.39) 

and  A is  a parameter  depending  on  the  type  and  size  of  the  switching  non- 
idealities. Then  there  exists  a positive  number  y such  that 

|x(t)  - x ( t ) | * yA  (4.40) 

for  all  t € [t^,“)  where  x(t)  is  the  motion  of  the  VSS  (4.17),  (4.18)  in 
sliding  mode  characterized  by  the  equations  of  sliding  mode  (4.24),  (4.25). 

If  the  trajectory  moves  to  some  A-neighborhood  of  s (x)  = 0, 
then  this  theorem  guarantees  that  it  will  remain  in  some  vicinity  of 
s(x)  = 0.  In  other  words,  if  the  reaching  conditions  hold  in  the  presence 
of  switching  nonidealities,  then  the  motion  is  close  to  the  ideal  sliding 
mode.  Furthermore,  (4.40)  implies  that  if  ideal  sliding  mode  x(t)  is  asymptoti- 
cally stable  in  the  subspace  s(x)  = 0,  then  x(t)  is  at  least  stable  in  the 
sense  of  Lyapunov. 

In  general,  when  the  VSS  (4.17),  (4.18)  is  in  sliding  mode, 

its  motion  governed  by  the  equations  of  sliding  mode  (4.24),  (4.25)  is 

influenced  by  the  disturbance  f.  If  there  are  uncertainties  in  the  system 

papameters  A,  B,  these  parameter  variations  also  affect  the  behavior  of 

sliding  mode.  Let  A = A + 5A  and  B = B + SB  when  A , B denote  the 

o o o o 

nominal  values,  then  providing  that  K(Bq  + 5B)  is  nonsingular,  the  equiva- 
lent control  is 

u = [K(B  + (SB)  ]_1K[  (A  + 6A)  x + E f]  (4.41) 

eq  o o o 

and  (4.24)  becomes 

x = (I  - (B  + <5B)[K(B  + 6B)]-1K)[(A  + 6A)x  + E f ] . (4.42) 

0 0 o o 


In  [39| , the  following  conditions  of  invariancy  of  sliding  mode  with  re- 
spect to  plant  parameter  variations^  and  disturbances  is  given. 

Lemma  4.5: 

If 


rank[B  + 6B  ; E ] = rank  (B  + SB)  (4.43) 

o o o 

and 

rank[B  + SB*  5A]  = rank  (B  + SB)  (4.44) 

o • o 

then  the  equations  of  sliding  mode  (4.42),  (4.25)  do  not  depend  on 
matrix  5A  and  the  disturbance  f. 

Conditions  (4.43),  (4.44)  imply  there  exist  matrices  P and  Q such  that 

6A  = (B  + SB)P  and  E = (B  + 5B)R,  then  (4.42)  becomes 
o o o 

x = (I  - (B  + SB) [K(B  + 5B)]-1K)A  x (4.45) 

0 0 o 

which  is  clearly  invariant  with  respect  to  plant  parameter  variations  SA 
and  disturbance  f.  We  note  that  conditions  (4.43),  (4.44)  resemble  the 
conditions  of  arbitrary  small  sensitivities  (3.17),  (3.18)  using  high 
gain  feedback.  The  relationships  of  VSS  with  sliding  mode  and  high  gain 
feedback  systems  are  explored  in  details  in  the  following  section. 

4.2.  Sliding  Mode  and  High  Gain  Feedback 

A type  of  nonidealities  in  switching  devices  is  the  finite  slope 
switching.  In  [40,41],  it  is  demonstrated  that  an  ideal  switching  device 
can  be  represented  by  an  amplifier  with  an  infinitely  large  gain  and 
saturation.  Thus,  it  is  reasonable  to  conjecture  that  the  behavior  of  VSS 


This  means  variations  in  matrix  A only. 
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in  the  neighborhood  of  the  switching  planes  is  close  to  that  of  high  gain 
feedback  systems.  This  conjecture  is  verified  in  the  following  theorem. 


Theorem  4.2: 


In  the  VSS  given  by  (4.17),  (4.18),  let  A = A^,  B = Bq  and 
f = 0.  If  sliding  mode  occurs  on  s = Kx  = 0 (4.22),  then  the  motion  of 
system  (4.17),  (4.18)  is  identical  to  the  limiting  motion  of  the  slow 
subsystem  (2.28)  of  the  high  gain  feedback  system  (2.1),  (2.7)  as  gain 


Proof:  By  applying  the  transformation  M introduced  in  (2.9)  in  section 
2.1,  system  (4.17)  with  A = Aq,  B = Bq  and  f = 0 becomes  (2.51),  (2.52). 
In  the  new  coordinates,  the  equivalent  control 


u = - (K2B°)  1[(K1A11  + K2A21)x^  + (kxA12  + K2A22)x2] 


(4.46) 


and  the  switching  planes  are 


s = K xj.  + K2x2  = 0 


(4.47) 


By  substituting  into  (2.51),  (2.52),  the  equations  of  sliding  mode  are 
given  by 


*'  _ »°  « + A°  x' 

1 A11  1 12X2 


(4.48) 


s = 0 


(4.49) 


and  (4.47).  Since  K^exists,  we  can  eliminate  x2  in  (4.48)  which  gives 


*1-  (AU  - A12K21k1)xI 


(4.50) 


i 


1 

I 


Hence,  the  equations  of  sliding  mode  can  be  represented  by  an 
(n-m)4^  order  system  (4.50).  We  call  this  the  equivalent  control  system. 

If  all  the  eigenvalues  of  KBq  have  negative  real  parts,  then  (4.50)  coin- 
cides with  the  slow  subsystem  (2.28)  of  the  high  gain  feedback  system 
(2.1)  (2.7)  as  g -*•  00 . For  the  case  when  KB^  has  unstable  eigenvalues, 

Lemma  4.4  shows  that  the  equivalent  control  system  is  invariant  with  re- 
spect to  switching  plane  transformation.  Since  KBq  is  nonsinguiar,  there 
exists  a nonsingular  matrix  Hg  in  (4.31)  such  that  ReX(HgKBo)  < 0.  This 
completes  the  proof. 

This  theorem  establishes  the  relationship  between  VSS  and  high 
gain  feedback  systems  based  upon  the  existence  of  sliding  model  on  one 
hand  and  the  asymptotic  stability  of  the  fast  motions  on  the  other  hand.  As 
gain  g -*■  00 , their  motions  in  the  null  space  of  matrix  K are  identical 
whereas  outside  this  null  space,  the  motion  of  the  high  gain  feedback  sys- 
tem occurs  infinitely  fast  in  the  range  of  B while  the  motion  of  VSS  de- 

o 

pend  on  the  switching  gains  r^  in  (4.21).  The  theorem  also  allows  us  to 
interpret  the  nonzero  eigenvalues  of  the  equation  of  sliding  mode  as  the 
transmission  zeros  of  the  open  loop  system  (4.17)  with  "output"  s.  Con- 
sequently, by  application  of  the  results  of  Section  2.3  and  2.4,  the  syn- 
thesis of  switching  planes  becomes  an  n-m  order  pole  placement  problem 
or  alternatively  a quadratic  minimization  problems  (problem  "s"). 

In  the  presence  of  system  parameter  variations,  we  see  from 
Section  3.1  that  the  stability  of  the  fast  motion  is  only  robust  in  the 
sense  of  (3.20).  From  (4.41),  the  invertibility  of  the  matrix  K (B  + <5B) 


I 


is  required  for  the  equivalent  control  to  be  unique  and  consequently  the 
equations  of  sliding  mode  to  be  unique.  From  Lemma  4.5  the  sensitivities 


of  sliding  mode  with  respect  to  parameter  variations  and  disturbances  de- 
pend on  range  conditions  of  5A,  6B  and  Eq  analogous  to  the  conditions  for 
high  gain  feedback  systems  given  in  Section  3.1.  In  view  of  the  equiva- 
lence of  sliding  mode  and  the  slow  motion  in  high  gain  systems,  we  obtain 
the  following  new  conditions  of  invariancy  for  sliding  mode. 


Theorem  4.3: 

In  the  VSS  given  by  (4.17),  (4.18)  let  A = A + <5A,  B = B + <$B 

o o 

and  suppose  sliding  mode  exists  on  s = Kx  = 0.  Then  if 


«Eo>  = «V 


R(<5A)  c R(b  ) 
— o 


R(SB)  c R(Bo) 


then  the  equivalent  control  system  is 


X1  (A11  ~ A12K2  Ki)xi 


which  is  identical  to  (4.50). 


(4.51) 


(4.52) 


(4.53) 


(4.54) 


Proof:  If  (4.51)— (4.53)  holds,  then  using  the  same  transformation  M as 


in  theproofof  Theorem  4.2  (4.17)  is  transformed  to 


Y ^ s A ^ y * 4>  y * 

1 11  1 12X2 

X2  = A21X1  + A22X2  + ^B2  + fiB2^u  + E2f 


where  MEq  = [0  E°]  . The  equivalent  control 

Ueq  = “ (K2(B2  + b2)>-1^KiA11  + K2A21)X1  + (K1A21  + K2A22)x2 


(4.55) 


(4.56) 


+ K2E2f], 


(4.57) 


By  substituting  u into  (4.56),  we  obtain  s = 0.  Eliminating  in 
(4.55)  (4.54)  is  obtained. 


Theorem  4.3  states  the  invariancy  conditions  in  terms  of  the 
range  of  Bo  while  in  Lemma  4.5,  the  condition  on  range  of  6B  is  not 
examined.  We  see  that  the  range  conditions  (4. 51)- (4. 53)  for  the 
invariancy  of  sliding  mode  are  identical  to  conditions  (3.17)— (3. 19) 
for  high  gain  state  feedback  systems.  Therefore,  we  conclude  that 
VSS  with  sliding  mode  and  high  gain  state  feedback  systems  reject  the 
same  class  of  disturbances  and  are  insensitive  to  the  same  class  of  para- 
meter variations. 


4.3.  VSS  with  Luenburger  Observer  Feedback 

In  single  input  VSS  in  phase  canonical  form,  the  variable 
structure  control  is  a function  of  the  output  and  its  derivative. 

There  have  been  studies  on  the  construction  of  estimators  for  variable 
structure  feedback  [42,43].  Also,  attempts  to  design  variable  structure 
estimators  have  been  made  [44,45,46,47].  Unfortunately,  extension  of 
these  results  to  multivariable  VSS  is  not  plausible  yet  we  need  them  in 
the  case  when  some  of  the  states  are  inaccessible.  Our  approach  is 
to  construct  the  well  known  full  order  Luenburger  observer  and  examine  the 
behavior  of  the  VSS  with  observer  feedback. 

The  svstem  considered  is 


x ■ A x + B u 
o o 


y - C x 
o 


(4.58) 


(4.59) 


91) 


whore  x t KU,  u t k"'  and  the  measured  outputs  v t K1  as  in  Chapter  3. 

Wo  assume  that  A , B and  C aro  known  matrices.  Then  a full  order 
o o o 

Luenburger  observer  la  of  the  form  1481 

X A A 

x * A x HI  u - l!(v  - v)  (4.60) 

o o * 

v - C x (4.61) 

o 

a n a m 

where  \ t K Is  the  observer  states  aiul  y ^ K Is  the  observer  outputs, 
lot  matrix  H bo  designed  such  that  the  error  aya tom 

o'  - (A  + HC  )e'  (4. 6.’) 

o o 

» x 

with  e x - x (4.bJ) 


ho  aaymptot leal  ly stable  and  possesses  desirable  characteristics.  lot 
tin*  variable  s true  tore  control  be 


u(x) 


i'°x  + $°y 


(4 . os  l 


whore 


V"  " Vdlagfagn  $j,...,sKn 


, ' . . . t h .a 

a tul  v ilenotes  the  i component  or  x. 


* - v 


, s(ix)  > 0 


.1-1 


v-rt  • Vx) 


V s . o S 1 


1 4 . no) 


v s . o 7 ) 


CO  > 


■ ■ *"  


—nr — ’ T" 


X 

1 

J 


I 


91 


and  s^Cx)  ■ 0 are  the  switching  planes.  We  define  an  m-vector  s(x)  whose 
,th 


i component  ia  denoted  as  s^(x)  by 


s (x)  * [s. (x) , . . . ,s  (x)  ]T  - Kx  - 0 
1 m 


(4.68) 


Similarly, 


<J>  - $ diag(sgn  y^ , . . . , sgn  yf) 


(4.69) 


and  y^  denotes  the  iC^  component  of  y. 


4>  — ) 

1 m 


r Pi  . ^ 


, §,(«)  > 0 


, i » 1 , . . . ,m 


(4.70) 


(4.71) 


-P4  , st(x)  < 0 


In  (4.66),  (4.67),  ^ and  r^  are  n-vectors  and  in  (4.70),  (4.71),  4^ 


and 


are  r-vectors. 


Lemma  4.6: 


Let  u be  a variable  structure  control  law  given  by  (4.64), 
then  there  exists  vectors  r^  and  p^,  i * l,...,m  such  that  the  tra- 
jectories x(t)  of  the  VSS  (4.60),  (4.64)  move  towards  the  intersec- 
tion of  the  switching  planes  s *■  K5fc  * 0 and  sliding  mode  exists  on 
“ 0. 

Proof : (4.60),  (4.68)  gives 


S - K(Aq  + HC^)x  + KBqu  - KHy . 


(4.72) 


By  considering  y as  measurable  disturbances,  we  construct  a variable 
structure  control  law  of  the  form  of  (4.64)  where  y is  "fed  forward" 
with  switching  gains.  For  (4.64),  the  condition  (4.35)  of  the  hier- 
archy of  controls  method  becomes 


* r /\  T I ^ I T | , , < . r T~  A I'  T i+1 

a [r  lx|  + p I y I ] - - min[p  x + t y + q u 
ii  l i+1 

u 


(4.73) 


where  the  minimization  in  (4.73)  can  be  performed  component  wise  analogous 
to  (4.36).  This  concludes  the  proof. 

Since  the  switching  planes  (4.68)  are  defined  in  the  observer 
state  space  x,  we  have  to  examine  the  behavior  of  x(t)  when  sliding  mode 
occurs  on  s(x)  = 0. 


Lemma  4 . 7 : 

Let  t > 0 be  the  time  instant  when  sliding  mode  occurs  on  s = 0 

r 

and  let  t ^ > 0 be  the  time  instant  such  that  for  some  positive  number  A 


le'(t) 


< A 


irw 


(4.74) 


for  all  t [t^,  00 ) , then  the  motion  x(t)  of  (4.58)  satisfies 


| K.x  ( t ) | ^ A 


(4.75) 


for  all  t e [max(t^,tr),  °°) 


Proof:  Let  e be  defined  as 

* ""  s 


- „ I 

e = no  . 
s 


(4.76) 


"HI11*  W '"R**\* 


93 


Assume  that  t > t . Then  for  t € [t  Kx(t)  - e (t).  Since  for  all 

6 r r s 

t,  | (t)  | * ||  K||  | e(t)  | , (4.74)  implies  |eg(t)|  5 A for  all 

t € [t  ,<*>) . Hence  (4.75)  is  obtained.  For  the  case  when  t > t , (4.76) 
e re 

implies  that  Kx(t)  ■ eg(t)  + KX(t).  By  the  triangular  inequality  and 


| Kx(t) | - A + | Kx(t) | , for  t € [te,°°)  (4.77) 

Since  by  definition  |Kx(t)|  - 0,  t 6 [tr,ro),  (4.77)  results  in  (4.75). 

The  Lemma  is  proved. 

By  defining  s - Kx  as  in  (4.22),  this  lemma  shows  that  x(t)  is 
confined  to  a A-neighborhood  of  s - 0 after  sliding  mode  occurs  on  s * 0. 
Thus,  the  behavior  of  VSS  with  observer  feedback  is  analogous  to  that  of 
VSS  with  switching  nonidealities.  The  following  theorem  is  a direct 
application  of  Theorem  4.1. 

Theorem  4.4: 

Let  t , t and  A be  as  defined  in  Lemma  4.7  and  let  x(t) 
e r 

be  the  motion  of  the  VSS  (4.58),  (4.64)  and  x(t)  is  the  motion  of  the 
observer  system  (4 . 60) , (4. 61) . Suppose  that  sliding  mode  exists  on 
s ■ 0,  then  there  exists  a positive  number  y such  that 

|x(t)  - ft(t)|  - Y A (4.78) 


for  all  t € [max(tg,  tr),  <*>) . Moreover, 

lim  |x(t)  - x ( t) | = 0 

t-«» 


(4.79) 
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Proof : (4.78)  ts  obtained  by  Theorem  4.1.  Since  lira  e(t)  - 0, 

£ ou 

hence  A *■  0 as  t ■*  00  and  this  concludes  the  proof. 

Clearly,  if  lira  $(t)  m 0,  then  the  VSS  (4.58),  (4.64)  is 

£ oo 

asymptotically  stable. 


Lemma  4.8: 

If  sliding  mode  occurs  on  $ *■  0,  then  the  equ.»t  Ions  of  sliding 
mode  Is  given  by 


x - (I  - H (KB  >-1K)[A  x - HC  e' ) 
o o * O o 


Kx  - 0. 


(4.80) 


(4.81 ) 


Proof:  Bv  setting  the  right  hand  side  of  (4.72)  to  rero  and  solving 

for  u , 


u - - (KB  )~l[KA  x + KtlC  e' 
eq  o o o 


substituting  (4.82)  in  (4.60),  we  have  (4.80),  (4.81). 

To  eliminate  m variables  in  Kx  • 0,  we  Introduce  the 
crausformat ion  M of  (2.-*)  as  in  the  proof  v> t Theorem  4.2.  let 

k 


(4.82) 


MH 


1 

U, 


(4.81) 


where  11^  is  an  (n-m)  x r matrix.  The  next  theorem  shows  the  "separa- 
tion property"  of  VSS  with  observer  feedback. 

Theorem  4.5: 

If  sliding  mode  occurs  on  s($)  » 0 (4.t>8),  then  the  equivalent 


control  s vs  tom  Is 


5;  - <*„  - Ai2K21|ti)*i  + HiV 


e'  - (A  + HC  )e’ 
o o 


(4.84) 


(4.85) 


J A 

where  Mx  - [x^,  x^  ] and  x^  is  an  (n-m)  vector  and  M is  as  defined  in  (2.9), 

Proof : Applying  transformation  M on  (4.60)  and  following  the  proof  of 
Theorem  4.2,  we  obtain  (4.84). 

From  the  block  triangular  structure  of  system  (4.84),  (4.85),  we 

see  that  its  eigenvalues  are  A(A°  - A°  kA.)  and  \(A  + HC  ). 

ll  iz  z L o o 

Furthermore,  the  n-m  eigenvalues  \(A°^  - A^K^Kj.)  are  identi- 
cal to  the  eigenvalues  of  the  equivalent  control  system  (4.50) 
when  the  variable  structure  control  feeds  back  all  the  states  and  switches 
on  s = Kx  « 0.  Therefore,  we  can  design  matrix  K as  if  all  the  states 
are  accessible,  separately  from  the  design  of  matrix  H.  This  property 
is  analogous  to  the  separation  property  of  linear  observer  feedback 
systems  [48].  We  note  that,  in  general,  sliding  mode  does  not  occur  on 
s « Kx  ■ 0.  Hence,  the  motion  x(t)  is  close  to  x(t)  In  the  sense  of 
(4. 78) , (4. 79)  when  sliding  mode  occurs  on  s * 0.  But  if  A^  - A°,,K,A^ 
is  a stable  matrix,  then  the  motion  x(t)  approaches  the  subspace 
s • Kx  ■ 0 and  simultaneously  moves  coward  the  origin  asymptotically  as 
t -*•  What  remains  is  to  show  the  behavior  of  x(t)  before  x(t)  reaches 
s - 0. 


Lemma  4.9: 

Let  t be  as  defined  in  Lemma  4.7  and  let  a = max  Re\(A  + HC  ) 
r oo 

and  3 ” max  Re\(Al>  ~ A^.,K.,^K,),  then  there  exists  positive  numbers  o 
c-  and  c..  which  depend  on  initial  conditions  x(0)  and  x(0)  such  that 


96 


|x(t)|  ^ cl  + c2eClt  + c3oBt  (4.86) 

for  all  t c [0,  t ). 

Proof : Applying  transformation  M on  (4.60)  and  if  sliding  mode  does 
not  occur  on  s = 0, 


A 

X 


1 


(A11  " Ai2K21ki)51  + A12K51S  + HlCoe' 


(4.87) 


Since  the  variable  structure  control  (4.64)  is  designed  so  that  the 

motion  x(t)  noves  towards  s(x)  = Kx  = 0,  there  exists  a positive  number 

c^  such  that  [s(t)|  - cj.  Consequently,  (4.S7)  implies  that  there 

exists  positive  numbers  e'  c'  and  c'  such  that 

- J 4 


| x(t) j £ c'  + c^e 


at 


C4e 


St 


(4.SS) 


(4.85)  implies  that  there  exists  a positive  number  such  that 


| e ' ( t ) | * c'e3t  (4.89) 

Since  x = x + e'  and  x = P^x^  + tor  S0lr-e  matrices  P^  and  POJ  by  the 
triangular  inequality  and  (4.88),  (4.81),  (4.86)  is  obtained 

We  see  that  since  matrices  A°  ^ and  AQ  + HC^  are 

stable,  for  bounded  initial  conditions,  x(0)  and  x(0),  x(t)[  - c^ 

bounded  by  an  exponentially  decay  time  function. 

4.4.  Variable  Struct  ire  Observer  Design 

The  relationship  of  high  gain  feedback  systems  and  VSS  with 


sliding  mode  suggests  the  idea  of  an  observer  with  a variable  structure. 
Such  an  observer  will  possess  similar  insensitivity  property  of  the 
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high  gain  observer  when  sliding  mode  occurs  on  appropriately  chosen 
switching  planes.  The  system  we  consider  is  the  system  (3.1),  (3.2) 


x - (A  + i$A)x  + (B  + SB)u  + E f 
o o o 


(4.90) 


y - (Cq  + 5C)x. 


(4.91) 


Recall  that  A , B and  C denote  the  nominal  values  and  5A,  5B,  SC 
o o o 

represent  the  parameter  variations  and  f is  the  disturbance  which  is  not 
measurable.  Let  a variable  structure  observer  be  of  the  form 


x-AS+Bu-Hv 
o o 


(4.92) 


where  v is  a variable  structure  feedback  discontinuous  on  the  switching 
planes 


s « C (x  - x)  = C e' 
e o o 


(4.93) 


that  is , 


\ (S,y) 


V[(a,y) 


^(S.y) 


s (e’)  > 0 
e 


Sg(e')  < 0 


» 1 “ 1 , . . . , r 


(4.94) 


where  v^  and  s*  are  the  component  of  the  r- vectors  v and  s^,  respec- 
tively. We  note  that  unless  <5C  = 0 in  (4.91),  we  cannot  determine  the 

signs  of  s*(e')  from  the  available  measurements  v(t)  and  $(t) 
e 

Although  for  a different  reason.  Theorem  3.3  shows  that  5C  = 0 
is  also  required  in  high  gain  observer  design.  Thus,  for  the  design 


i 


variable  structure  observer,  wo  make  the  assumption  that  the  measure 


ment  matrix  is  known,  C = C , that  is, 

o 


SC  - 0. 


(4.95) 


Moreover,  we  assume  that  the  matrix  C^H  is  nonsingular.  This  assures 

the  uniqueness  of  the  equation  of  sliding  mode.  We  now  examine  the 

behavior  of  system  (4.92)-(4.94)  when  sliding  mode  occurs  on 

s » C e ' « 0 . 
e o 

Theorem  4.5: 

Let  >1  be  the  transformation  defined  in  (3.38)  and  e = 

7 

[c,  e,]  as  in  (3.37).  Suppose  sliding  mode  occurs  on  s = C e «=  0, 

l ‘ e o 

then  the  equivalent  control  system  is 


~ a -1 

• . ,o  ,o  „o  ,0, 

°1  ~ ^Ail  ■ A12C2  Cl°l 


(4.  %) 


if  and  onlv  if 


R(«A) 


(4.97) 


* ( -B) 


(4 . 98) 


R(E  )C  R(H) 
o — 


(4.99) 


Proof : From  (4.63),  (4.90)  and  (4.92),  we  obtain 


o'  = A e'  + IN’  + <$Ax  + jBu  + E f 
o o 


(4.100) 


Applying  the  transformation  M to  this  equation,  we  obtain  (3.391,  (3.40) 


Solving  for  the  equivalent  control  v 


and  substituting  into  (3.40),  the 


equivalent  control  system  Is 


• l - (A^  - ApC°  C^)el  + 6Atx  + vSB^a  + E^f  (4.101) 

A N 

The  conditions  (4.97)  - i4.99)«re  equivalent  to  $Aj  » 0,  >58^  - 0 and 
- 0.  respectively.  This  concludes  the  proof. 

We  note  that  (4. 97)- (A . '59')  are  the  same  range  conditions  iJ.47)- 
(3.49)  for  insensitivity  in  high  gain  observer.  Moreover,  the  eigen- 
values of  system  (4.9b)  coincide  with  the  eigenvalues  of  the  slow  sub- 
system of  the  high  gain  feedback  system  (,3.36)  as  g -*■  *.  Theorem  3.2 
gives  the  next  lemma. 


I.emiTu^  4.10 : 

Let  (A  ,0  ) be  an  observable  pair,  then  there  exists  a matrix 
'o°  °'o  o'1 

H such  that  \(A^  - A^  ,C\  C^)  are  placed  to  q^,  j - 1 n-r 

where  q ^ are  the  prescribed  locations  of  the  eigenvalues  of  the  equiva- 
lent control  system  of  the  VSS  (4. 92)-(4.94) , 

It  remains  to  determine  the  feedback  functions  v*.  in  the 
variable  structure  feedback  law  (4.94).  Using  again  the  transformation 


s - F’fe,  + F"'s  + C Hv  + C ,SAx  + C «5Bu  + C E f 
el!2eo  o o oo 


(4.102) 


where 


'o  .o)o  , ,o)o  'o.o 
F1  * C1A11  + "2A21  “ F2C1 


t4. 103) 


F2  ’ (C1A12  + C2A22^  C2 


^*4 . 1 0 *4  ^ 
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Since  only  and  u are  accessible,  we  have  to  consider  e^ , x and  f as 
unmeasurable  disturbances.  Then  a feasible  variable  structure  control 
is 


v(s  ,u)  = '{'  S + II  u + TT, 
e v g v b 


(4.105) 


where 


¥°  = V diag(sgn  s\...,sgn  sr) 
v v ° e e' 


(4.106) 


H = Hdiag(sgn  u ,...,sgn  uj 


II  = (tt  ,...,irr) 


(4.107) 

(4.108) 


and  u. ,s ^ denote  the  i^1  component  of  m-vector  u and  r-vector  s , 
l e e 


respectively , 


0 , s (e’)  > 0 

i e 


-a1  ,s*(e')  < 0 


, i ~ 1 , . . . , r , b. 


(4.109) 


In  (4.108),  (4.109),  tt  and  a , i = 1,  — ,r  are  m-vectors  and  tt,  and  o. 


are  r-vectors. 


Lemma  4.11 : 

Let  v be  a variable  structure  control  given  by  (4.105),  then 

there  exists  switching  gain  matrices  and  vectors  0^,  i = l,...,r,b  such 

that  the  motion  e'(t)  moves  towards  the  intersection  of  the  switching 

planes  s = C e'  = 0 (4.93)  and  sliding  mode  exists  on  s =0. 
r n C 


T 


»»  1 
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Proof:  Hv  considering  u ns  measurable  disturbances  and  e^,  x and  f as 

unmeasurable  dint  urbane**  and  apply  lug  the  hierarchy  of  controls  method. 
Inequalities  similar  to  (4.35)  and  (4.73)  are  obtained  for  the  variable 
structure  control  given  by  (4.105). 

Hv  this  lemma,  we  have  concluded  the  design  of  variable  struc- 
ture observer. 

4.5.  VSS  With  Variable  Structure  Observer  Feedback 

Thus  far,  we  have  shown  that  when  all  the  states  are  accessible, 
VSS  with  sliding  mode  reject  the  same  class  of  disturbances  and  parameter 
variations  as  In  high  gain  state  feedback  systems.  The  same  conclusion 
Is  made  regarding  variable  structure  observers  and  high  gain  observers. 
Therefore,  It  Is  expected  that  VSS  with  variable  structure  observer 
feedback,  that  Is,  a VSS  with  two  variable  structure  loops,  possesses  the 
same  Insensitivity  property  as  high  gain  system  with  high  gain  observer 
feedback  studied  In  Section  3.4.  in  this  class  of  variable  sturcture 
feedback  systems,  u Is  a variable  structure  control  tot  the  svstem 
t4.90),  (4. 913  of  the  form, 

(u^(x,y)  If  s^tx)  > 0 

, i - 1, . . . ,m  (4.110) 

u^(x,y)  if  Sj($)  < 0 

where  t ho  switching  planes  are  as  defined  (n  (4 . t>S) , 

s (x)  - Kx  - 0 (4.111) 

and  x Is  the  state  of  the  variable  structure  observer  (4 . 92)- (4. 94) . Kor 


the  same  reason  given  in  the  last  section,  we  assume  SC  - 0 in  (4.91). 
We  note  that  there  are  two  sets  of  switching  planes.  The  set  s(.<)  - 0 
is  defined  in  the  space  of  x and  the  set  se(e')  - 0 is  defined  in  the 
space  of  e’.  We  now  examine  the  behavior  of  the  variable  structure 
feedback  system  (4. 90)- (4. 94) , (4.110),  (4.111)  when  sliding  mode 
occurs  on  s ■ 0 and  se  ■ 0.  We  first  determine  the  equivalent  control. 


Lemma  4.12: 

If  the  (nri-r)  * (nH-r)  matrix 


l = 


KB  - KH 

o 


C SB  C H 
o o 


(4.112) 


is  nonsingular,  then  the  equivalent  control  for  sliding  mode  occurs 
on  s * 0 and  sg  ■ 0 is  unique  and  is  given  by 


u 

KB  - KH' 

-1 

KA  0 

r» 

' 0 

, 

o 

0 

m — 

+ 

V 

CiB  CH 

C 6A  C (A  + 6A) 

e’ 

C E 

eq 

0 0 

o 

o 

o 

o o 

*■*  -* 

. 

L J 

L . 

f (4.113) 


Proof : If  y is  nonsingular,  then  the  equivalent  (4.113)  can  be  solved 

uniquely  in  the  simultaneous  algebraic  equations  a - 0 and  se  ■ 0. 

The  corresponding  conditions  for  (4.112)  in  high  gain  observer 
feedback  systems  is  (3.82). 


Theorem  4.7: 

A 

Let  M and  M be  the  transformations  defined  in  (2.9)  and  (3.38) 

T T 

respectively  and  e ■ [e^  e,]  as  in  (3.37).  Denote  [x^  ] ■ M$.  Let 
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I 

I 


i r>, 


be  as  defined  in  (4.83)  and  N as  in  (2.19).  Suppose  that  £ in  (4.112) 
is  invertible  and  sliding  mode  occurs  on  s ■ 0 and  sg  • 0,  then  the 
equivalent  control  system  is  represented  by 


e - (A°  - A°  C°_1C°)e 

1 U11  A12  2 V 1 


x,  - (A°  - A°  + H.  (C  H)_1C  A Ne 

1 11  12  211  lo  ool 


(4.114) 


(4.115) 


if  and  only  if 


R(E  ) c R(B  ),  R(E  ) C R(H) 
o — o o — 


(4.116) 


R(6A)  c R(B  ),  R(6A)  c R(H) 
— o - 


R(5B)  ^ R(Bq),  R(5B)  5 R(H) 


Proof:  Let 


ft  = (KB  + KH(C  H)  C 6B) , 
o o o 


(4.117) 


(4.118) 


(4.119) 


then  (4.113)  cqn  be  written  as 


u - n_1[(KH(C  H)_1C  5A  - KA  )x 
eq  ooo 

- KH(C  H)-1  (C  (A  + SA)e'  + C E f)]  = U.x  + U.e'  + U-f  (4.120) 
O oo  00  1 i 0 

v - - (C  H)-1  { [C  5Bfi_1(Kll(C  H)-1C  6A  - KA  ) + C 6A]x 
eq  o o o o o o 

- [C  $Bn_1(I  - KH(C  H)_1C  ] [ (A  + <SA)e'  + E f]  = V x + V0e'  + V-f  (4.121) 

Q is  nonsingular  by  the  assumption  that  £ is  nonsingular.  Substituting 
u^  and  v for  u and  v respectively  in  (4.92)  and  (4.100),  the  equation 
of  sliding  mode  is  obtained. 


x - (A  + B u - HV  )x  MB  u,  - HV„)e ' + (B  U~  - HV.)f  (4.122) 

O O I 1 O - Z Oj  J 

Kx  - 0 (4.123) 

4’  - (6BU  + HV.  )x  + (A  + fiBU,  + HV.,)e'  + (E  + 5BU  + HV,)f  (4.124) 
11  o Z Z oJJ 

C e’  - 0.  (4.125) 

o 


The  conditions  r(6B)  C r (h) , R(5A)  c r(h)  and  r(Eq)  c r (H) 
are  equivalent  to  the  identities  H(CoH)  ^Co<SB  - iSB,  H(CqH)  ^Cq6A  - 5A 
and  H(C^H)  C^E^  “ E • respectively  which  simplifies  (4.122), 

(4.124)  to 

x , 

x - [A  + (B  + 5B)U'  + 6Alx  + [ (B  + $B)U'  + H(C  H)  C A + 5A]e’ 
OO  1 o Z o 0 o 

+ [ (B  + 6B)U  ' + E ]f  (4.126) 

o Jo 

e - (5BU  *]x  + [A  + 6BU ' - H(C  H)_1C  A - 6A]e  ' 

1 o Z o o o 

+ [2E  + <5BU  ' ] f (4.127) 

o 3 

where 

U,'  - (K(B  + 5B)]"1(K5A  - KA  ) (4.128) 

1 o o 

u!  - - [K(B  + 6B)]-1(KH(C  H)_1C  A + K5A)  (4.129) 

Z o o o o 

Ui  - - [K(B  + SB) ]_1KE  . (4.130) 

Jo  o 

Applying  the  transformation  M on  (4.126)  and  eliminating  s,  we  obtain 


| 

i 
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" (A11  " AlZK2lKl)*l  * + <5A)Qx^ 

+ M H(C  H)_1C  A e'  + M (<$BU'  + 6A)e'  (4.131) 

1 O O O 1 L 

+ M, (6BU*  + E ) f 
1 J o 

where  M^Bq  - 0 and  KQ  » 0 as  defined  In  (2.21)  and  (3.53),  respectively. 
Applying  the  transformation  M on  (4.127)  and  eliminating  se>  we  obtain 
-i 

e.  - (A°  - A°  C°  C°)e.  + M.  (SBU' - 6A)Ne.  + M.SBU’x  + M (2E  + 5BU!)f 

i 11  1 i 1 1 l ill  lo  3 

(4.132) 

whre  M^H  * 0 and  CqN  » 0 as  in  (3.94)  and  (2.19),  respectively.  The  con- 
ditions (4. 116)- (4. 118)  are  equivalent  to  M,E  * 0,  M,E  - 0;  M,  *5A  **  0, 

1 o 1 o 1 

M^6A  - 0;  and  M^6B  = 0,  M^SB  = 0 respectively.  By  eliminating  sg  from  e' 
in  (4.131),  the  equivalent  control  system  (4.114),  (4.115)  is  obtained. 
This  concludes  the  proof. 

We  note  that  conditions  (4. 116)- (4. 118)  coincide  with  the 

range  conditions  (3. 78)- (3. 80)  for  insensitivity  in  high  gain  observer 

feedback  systems.  Since  system  (4. 114) , (4. 115)  is  in  block  triangular 

form,  its  eigenvalues  are  \(A°^  “ A°?C°  C°)  and  A(A°^“  A°.-,K0‘''K^).  From 

(4.54)  we  see  that  \(A°^  - A^Kj^K^)  are  the  eigenvalues  of  the  equivalent 

control  system  when  the  switching  planes  are  s a Kx  * 0 whereas  from 

-1 

A © AQ  Q Q 

(4.96),  X(A^  “ A19^2  ^1^  are  t^ie  el-8enva^ue3  c^e  variable  structure 

observer  in  sliding  mode.  This  property  is  analogous  to  the  separation 
property  of  VSS  with  Luenburger  observer  feedback. 


In  general,  sliding  mode  does  not  exist  on  s(x)  - Kx  ■ 0. 
However,  it  can  be  shown  that  the  motion  x(t)  of  (4.90)  is  close  to 
x(t)  when  sliding  mode  occurs  on  s * 0 and  sg  * 0. 

Lemma  4.13: 

Let  tr  > 0 be  the  time  instant  when  sliding  mode  occurs  on  § * 0 

and  s » 0 and  let  t > 0 be  the  time  instant  such  that  given  a positive 
e e 

number  A 

|ei(t)|  < I!  ttir1*  (4.133) 

for  all  t 6 [te,  °°)  where  e^  in  the  vector  defined  in  Theorem  4.7  and 

matrix  N satisfies  C N * 0 as  in  (3.94).  Then  the  motion  x(t)  of  the 

o 

VSS  with  variable  structure  feedback  satisfies 

I Rx(t) I - A (A. 134) 

for  all  t € [max(te,  tr)  ,<») . 

Proof : Analogous  to  the  decomposition  in  (2.18),  e'  can  be  written  as 

e’  - Me-  + H(C  H)_1s  . (4.135) 

1 o e 

Assume  that  t > t , then  for  t t [t  , 00  ) Kx(t)  ■ KNe.,(t).  Since  for 
e r r i 

all  t,  | Kx(t)  | - ||  KN||  |e  (t)|,  (4.133)  implies  (4.134).  If 


» 


" 


I 
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t > t , then  for  t€[t,  “») , Kx(t)  * Oe.  (t)  + KH(C  H)  *s  + Kx(t) . 
re  e i o e 

By  the  triangular  inequality  and  (4,133) 

| K.x(t)  | 4 A + ||  KH(Coa)_1!|  |se(t)|  + |a(t)  | (4.136) 

for  t c It  , *).  Since  |s  (t) | - 0 and  |s(t)|  - 0 for  t 6 [t  , «) , 

8 c 1 

by  definition,  (4.136)  results  in  (4.134).  The  lemma  is  proved. 

Be  defining  s(x)  - Kx  as  in  (4.22),  we  see  that  x(t)  is  in 
a ^-neighborhood  of  s(x)  - 0 after  sliding  mode  occurs  on  s(x)  - 0 and 
s (e')  " 0.  This  phenomenon  is  analogous  to  that  in  VSS  with  Luenburger 
observer  feedback  and  can  be  considered  to  be  caused  by  switching  non- 
idealities.  Applying  Theorem  4.1,  we  express  x(t)  as  an  approxi- 
mation to  x(t)  in  the  following  Theorem. 

Theorem  4.8: 

Let  t , t and  A be  as  defined  in  Lemma  4.13  and  let  x(t)  be  the 
e r 

motion  of  the  VSS  with  variable  structure  feedback  given  by  (4. 90)- (4. 94) , 
(4.110)-(4.111)  and  let  x(t)  be  the  motion  of  the  variable  structure  ob- 
server system  (4.90)-(4.94).  Assume  that  Re\(A^  - A^.,C°  C^)  < 0. 
Suppose  that  sliding  mode  exists  on  s - Kx  - 0 and  s(  ■ C e'  ■ 0,  then 
there  exists  a positive  number  y such  that 

| x ( t ) - x ( t ) | - y A (4.137) 

for  all  t [max(t^,  tf) ,°°) . Moreover, 

lim  | x ( t ) - x ( t ) | - 0 (4.138) 

£•*00 


...  ,r, 


I 


Proof : (4.137)  la  a direct  consequence  of  (4.134)  and  Theorem  4.1. 

-1 

Since  the  matrix  A°^  - A°.,C°  is  stable,  lim  e^(t)  - 0.  Hence  A ■*  0 

t "►  00 

as  t + ® and  (4.137)  becomes  (4.138). 

From  the  equivalent  control  system  (4.114),  (4.115)  we  con- 
clude that  of  A^  - A^,,K9  is  a stable  matrix.  Thus  lim  x(t)  ■ 0 and 
Theorem  4.8  implies  that  the  motion  x(t)  goes  to  zero  asymptotically  as 
t ■*  00 . If  Theorem  4.7  holds,  then  x^(t),  the  motion  in  the  subspace 
Kx  ■ 0 is  invariant  with  respect  to  parameter  variations  and  disturbances. 
Since  x(t)  is  0(A)  close  to  $(t)  (in  the  sense  of  (4.137))  and  is  in  the  0(A) 
neighborhood  of  s(x)  - 0 (in  the  sense  of  (4.134)),  we  conclude  that  for 
t fe  [max(te,  tf)  ,°<» ),  the  sensitivities  of  x(t)  with  respect  to  parameter 
variations  and  disturbances  are  0(A). 

Analogous  to  Lemma  4.9,  we  now  show  the  boundedness  of  the 
motion  x(t)  and  e '(t)  before  sliding  mode  occurs  on  s • 0 and  sfi  ■ 0. 

Lemma  4.14; 

Let  t be  as  defined  as  Lemma  4.13  and  let 
r -1 

a - max  Re\(A^  - A^C^  c° ) and  8 - max  Re.\(A^1  - A^.,K?  ^K.^).  Suppose 
that  conditions  (4. 116)-(4. 118)  hold  then  there  exist  positive  numbers 
c^,  and  Cj  which  depend  on  initial  conditions  x(0)  and  $(0)  such  that 

| x( t ) | - + c^e0*1"  + c^e^C  (4.139) 

for  all  t € (0,  tr) . 

Proof:  Applying  transformation  M and  as  defined  in  Theorem  4.7  to 


(4.90),  (4.100).  Under  the  conditions  of  this  Lemma,  if  sliding  mode 
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does  not  occur  on  3 ■ 0 and  sg  - 0, 


A | ,.0  .0  1,,  \ A l . .0  ^"1^  ,, 

X1  " (A11  " A12k2  K1)X1  A12K2  S " H1V 


-1 


V <An  - ?UC(lc>  1 + AnC2  Se 


(4.140) 

(4.141) 


Since  the  variable  structure  control  u and  v are  designed  such  that  the 

motion  x(t)  moves  towards  s(x)  ■ Kx  ■ 0 and  the  motion  e'  (t)  moves  toward 

s (e')  - C e'  -0,  there  exist  positive  numbers  o'  and  o'  such  that 
e o 12 


s(t)|  ^ cj  and  | s ( t ) | ^ c^-  Moreover,  since  v(t)  is  bounded,  consequently. 


from  (4.140) 


|x’(t)|  - c^+  c^t 


(4.142) 


and  from  (4.141) 


|e(t)|*  c5  + V 


'o«t 


(4.143) 


Since  x * x + e'  and  x * P.x'  + P_s,  e'  ■ P-,e,  + P,s  for  some  matrices 

11  2 3 1 •*  e 

P^,  P.,,  P^  and  P^,  by  the  triangular  inequality  and  (4.142),  (4.143), 
we  obtain  (4.139). 

A „ -1 
o o „o  „o 

We  note  that  since  the  matrices  - A^.,C,  and 
A11  " A1'iK'’1ki  are  stable>  _C1  is  bounded  by  an  exponentially 

decaying  time  function  for  bounded  initial  conditions  x(0)  and  x(0). 

In  the  last  section,  the  variable  structure  control  v of  (4.105) 
in  the  variable  structure  observer  includes  feedforward  terms  of  u.  If 
u is  also  a variable  structure  control,  these  terms  are  not  necessary. 


Since  the  accessible  variables  include  x,  sg,  feasible  variable  structure 


control  schemes  for  u and  v are 


1 


I 


v($,S  ) - 4»  « + 4Ts  + TT, 
e v v e b 


(4.144) 


and 


u(x,s  ) ■ f x + 4's  + <p, 

e u e b 


(4.145) 


where  4'°,  ¥°  and  are  as  in  (4.106),  (4.65)  and  (4.105),  respectively 


and  4’v,  4*°  and  are  defined  as  follows. 


4*  - 4'  diag(sgn  x.,...,sgn  x ) 

v in 


^ m 


(4.146) 

(4.147) 


( 5,  , si(e’ ) > 0 


*1  " 


, s1 (e ' ) < O’ 
e 


i - 1, . . . ,r 


(4.148) 


4'°  ■ 4*  diag(sgn  s\...,sgn  sr  ) 


u u 


(4.149) 


and 


Si(x)  > 0 


A - /N  . . 

s1(x)  < 0 


.1-1, 


,m. 


(4.150) 


In  the  above  equations,  and  5^,  i - l,...,r  are  n-vectors,  and  6^ 
are  m-vectors. 


Lemma  4.15: 

Let  v and  u be  variable  structure  controls  given  by  (4.144) 


(4.145),  then  there  exists  switching  gain  matrices  ¥ ,4*v  and  vectors  r^. 


V V “b  Md  V 1 ' l' 


J 


m,  j - 1 


r,  such  that  e'  (t)  and  x(t) 
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move  towards  the  intersections  of  the  switching  planes  s - C e ' **  0 

e o 

and  s ■ Kx  ■»  0,  respectively.  Moreover,  sliding  mode  exists  on  * 0 
and  s » 0. 

Proof : By  considering  e^,  x and  f as  unmeasurable  disturbances  and 

applying  the  hierarchy  of  controls  method,  inequalities  similar  to  (4.35) 
and  (4.73)  are  obtained  for  the  variable  structure  controls  given  by 
(4.144)  and  (4.145). 

4.6.  Regulation  by  VSS 

We  conclude  our  study  on  VSS  with  sliding  mode  by  considering 
the  problem  of  regulation  of  some  specified  variables  z = D^x  in  the 
presence  of  system  parameter  variations  and  disturbances  using  variable 
structure  controls.  We  show  that  the  conditions  for  the  regulation  of 
z in  VSS  are  the  same  as  when  high  gain  feedback  control  is  employed. 

The  system  we  consider  is  system  (4.90),  (4.91),  the  same  one  in  (3.3). 

Let  z ■ Dqx  be  the  variables  to  be  regulated.  We  first  assume  that 
all  the  states  are  accessible  and  let  u be  the  variable  structure  control 
given  by  (4.18)  such  that  x(t)  moves  toward  s(x)  » Kx  « 0 (4.22)  and 
sliding  mode  exists  on  s - 0.  Then  the  condition  for  regulation  of  z(t) 
is  given  in  the  following  theorem. 

Theorem  4.9: 

Let  t > 0 be  the  time  instant  when  sliding  mode  occurs  on 
s(t)  * Kx  ■ 0 and  let  K » [K^K^  K^]  thenz(t)  * 0 for  c - cr  if  and  only  if 

Re\(Z^  + 52  + 5B1L)  < 0 (4.151) 


1 12 


i— 


■T 


!! 


' 


and 

<2?  + 5Z.  + 5B.L|R(M.  [I  - 5B(K(B  + 5B))-1K]E  )>  c sj(D  Q)  (4.152) 
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where  the  notations  used  are  the  same  as  in  Theorem  3.5. 

Proof : Applying  the  transformation  M of  (2.9),  system  (4.90)  is  of  the 
form  of  (3.9),  (3.10).  The  resulting  equivalent  control  system  for 
sliding  mode  on  s ■ 0 is 

x,’  - (Z?  + $Z  + 5B.  L)x'  + (M.  [I  - 5B(K(B  + 5B))-1K]E  )f  (4.153) 
llllil  o o 

Since  by  definition  s - u),  from  Theorem  3.5, 

z - DqQx^  + D^B^ (EB^)-1s  (4.154) 

Hence  z(t)  - DoQx|(C),  t - tr.  Using  Wonham's  result  [33],  the  condition 
(4.152)  is  obtained.  Condition  (4.151)  assures  the  VSS  in  sliding  mode 
is  asymptotically  stable.  This  proves  the  Theorem. 

We  note  that  conditions  (4.151)  and  (4.152)  are  the  same  condi- 
tions (3.33),  (3.55)  for  the  regulation  of  z(t)  using  high  gain  state 
feedback.  For  the  same  reasoning  given  in  Section  3.5,  (3.152)  is 
replaced  by 

rank  KT  - rank[DT  (4.155) 

o 

or  equivalently  there  exists  matrix  P such  that 

D - PK  (4.156) 

o 

An  interpretation  of  (-*.155)  is  that  regulation  of  z(t)  is  achieved  once 
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slide  mode  occurs  on  s = 0 which  can  be  seen  by  multiplying  on  both  sides 
of  (4.156)  with  x(t). 

In  the  case  when  only  y(t)  of  (4.91)  is  accessible,  we  assume 
that  <$C  = 0 and  the  matrices  6A,  5B  and  Eq  satisfy  (4.97)  - (4.99).  Let 
u and  v be  variable  structure  controls  given  by  (4.144),  (4.145)  such  that 
x(t)  moves  toward  s = 0 (4.68)  and  e'(t)  moves  toward  sfi  ■ 0 (4.93)  and 
sliding  mode  exists  on  s = 0 and  sg  = 0.  Theorem  4.6  shows  that  under 
the  assumptions  we  have  made,  the  motion  e'  (t)  is  invariant  with  respect 
to  parameter  variations  and  disturbances  in  the  subspace  sg  = 0.  We  now 
assume  that  this  sliding  mode  in  = 0 is  asymptotically  stable. 

Theorem  4.10: 

Let  t > 0 be  the  time  instant  when  sliding  mode  occurs  on 

a = Kx  = 0 and  s = C e * 0 and  let  a be  as  defined  in  Lemma  4.14, 
e o 

Z = A°x  - + Ml(6BU^  + 6A)Q  (4.157) 

and  B » max  X(Z),  then  there  exist  positive  numbers  c^  and  c9  such  that 

[z(t)|  S C]Leat  + c2e8t  (4.158) 

for  all  t € [t  , °°)  and  lim  z(t)  = 0 if  and  only  if 

t -*•  °o 

ReX(A°1  - A^K'/K,  + M1(6BU;[  + 5A)Q)  < 0 (4.159) 

and 

c(A°i  - + MX(6BU^  + 6A)Q)|R[M1(I-6B(K(Bq  + 6B))-1K)Eo]> 

c V(D  Q)  (4.160) 

— o 


i 


i^a)  Tanar*)  irun  mm 


where  the  notations  used  are  the  same  as  in  (4. 128)- (4. 131) , 


Proof : System  (4.131),  (4.132)  represents  the  equivalent  control  system 

of  the  VSS  with  switching  planes  s » 0 and  se  - 0.  The  conditions (4.97)- 

(4.99)  are  equivalent  to  M.SB  - 0,  M,  6A  ■ 0 and  M, E = 0.  Hence  (4.132) 

1 1 1 o 

is  reduced  to  (4.96)  and  there  exists  positive  such  that 
Ct  t! 

^ c^e  » c ^ tcr*  °°) . Since  x can  be  decomposed  analogously 
as  in  (2.18)  into 


x = Qx'  + B (KB  )_1s 
1 o o 


(4.161) 


and  e'  can  be  expressed  as  in  (4.135),  hence 


z(t)  = DQNe  (t)  + DoQ»* (t) 


(4.162) 


From  (4.131),  the  solution  x^(t)  is  given  by 


x'  (t)  = x'  (t)  + x1u (t) 


(4.163) 


where  xj^(t)  is  the  component  due  to  f(t).  Wonham's  result  [33]  shows 
that  x^(t)  = 0 if  and  only  if  (4.160)  holds.  The  component  x^(t)  is  due 
to  e^(t)  and  "initial  condition"  x^(tr)  and  there  exist  positive  and 
and  such  that 

|*1K(t)|  * cleat  + cle6t 


Applying  the  triangular  inequality  to  (4.162),  (4.158)  is  obtained. 


By  the  stability  assumption  on  the  observer  and  (4.159),  ct  < 0 and  8 < 0, 
(4.158)  implies  lim  z(t)  = 0.  The  theorem  is  proved. 

t-x» 


We  note  that  conditions  (4.159),  (4.160)  are  equivalent  to 
the  conditions  (3.101),  (3.102)  in  Theorem  3.11  for  the  regulation 
of  z(t)  using  high  gain  observer  feedback.  For  the  same  reasoning 
given  in  Section  3.5,  (4.160)  is  replaced  by  (4.155).  If  Dq  satisfies 
(4.155)  or  equivalently  (4.156),  then  from  (4.134), 

| z ( t ) 1 - ||  P ||  A (4.164) 


for  t > max(t  , t ) and  in  view  of  (4.138),  lim  z(t)  ■ 0.  By  comparing 
e r t -*■  00 

Theorem  4.9  and  4.10,  we  see  that  if  y(t)  is  the  only  information  avail- 
able, z(t)  is  regulated  asymptotically  whereas  it  is  zeroed  as  sliding 
mode  occurs  on  s(x)  *■  * 0 in  the  case  when  all  the  states  are  accessi- 

ble. 

Finally,  we  observe  that  the  synthesis  problems  K and  H dis- 
cussed in  Section  3.6  are  encountered  again  in  VSS  with  sliding  mode. 

They  become  the  synthesis  of  the  switching  planes  s * Kx  *=  0 and 
— T—  — 

sg  » H e = 0.  The  switching  planes  sg  = 0 are  constructed  for  the  dual 
of  the  VSS  (4.100)  with  <$A  = 0,  5B  - 0 and  f = 0,  that  is, 

e * A^e  + C^v  (4.165) 

o o 


where  the  components  v^  of  the  variable  structure  control  v is  dis- 
continuous on  the  components  s*(e)  * 0 of  sg(e)  = 0.  We  note  that  for 
VSS,  it  is  not  required  that  the  matrices  KBq  and  CqH  be  stable  as  in 

the  high  gain  feedback  case.  Instead,  they  should  be  invertible.  This 
change  in  specification,  however,  does  not  change  the  solvability  of 
the  synthesis  problems  K and  H of  Section  3.6. 
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In  this  chapter,  we  have  established  the  relationship  between 
VSS  and  high  gain  feedback  systems.  It  is  shown  that  the  insensitivity 
property  of  sliding  mode  in  VSS  and  the  slow  motion  in  high  gain  systems 
depends  on  the  same  range  conditions.  In  other  words,  these  two  types  of 
feedback  systems  reject  the  same  class  of  disturbances  and  parameter 
variations.  Furthermore,  sliding  mode  coincides  with  the  limiting  motion 
of  the  slow  subsystem  as  the  loop  gains  tend  to  infinity.  We  emphasize 
that  the  behaviors  of  VSS  and  high  gain  system  are  similar  only  in  the 
subspace  defined  by  the  intersections  of  the  switching  planes  which  coin- 
cides with  the  null  space  of  the  high  gain  feedback  matrix.  In  high 
gain  system,  the  motion  outside  this  subspace  is  extremely  fast  due  to 
the  large  gain  factor  and  its  behavior  depends  on  the  constant  feedback 
matrix.  On  the  other  hand,  the  motion  before  the  occurrence  of  sliding 
mode  is  not  necessarily  fast  and  its  behavior  is  determined  by  the  vari- 
able structure  control.  In  the  case  when  reduced  order  models  are  used 
in  the  design  of  high  gain  systems  and  VSS,  the  capability  to  reach  this 
subspace  is  susceptible  to  the  neglected  small  time  constants  if  the 
reaching  motion  is  forced  to  be  extremely  fast  by  the  use  of  large  feed- 
back gains.  However,  since  high  gain  is  not  necessary  in  the  design  of 
variable  structure  control  in  VSS,  reaching  can  be  guaranteed.  These 
aspects  regarding  robustness  property  of  high  gain  systems  and  VSS  with 
respect  to  actuator  and  sensor  dynamics  are  examined  in  the  next  chapter. 
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CHAPTER  5 
APPLICATIONS 

5 . 1 Systems  with  Actuator  and  Sensor  Dynamics 
< An  important  consideration  in  the  realization  of  high  gain 

feedback  systems  is  the  effects  of  actuator  and  sensor  dynamics.  In 
many  control  applications,  these  dynamics  are  neglected  in  order  to  reduce 
the  complexity  of  the  system  model.  Since  the  insensitivity  property 
of  high  gain  feedback  system  depends  solely  on  the  system  structure,  it 
is  necessary  to  examine  the  effects  on  the  actual  system  if  a reduced 
order  model  is  employed  in  the  high  gain  feedback  design.  In  this  sec- 
tion, we  shall  examine  only  the  effects  of  actuator  and  sensor  dynamics. 

We  assume  all  the  states  are  measurable  and  there  are  no  system  parameter 
variations  and  disturbances.  The  system 

x - Ax  + B u (5.1) 

o o 

is  considered  where  x t Rn  and  u c Rm.  This  is  the  system  model  we  have  been 
using  in  the  previous  chapters  and  the  high  gain  feedback  is  written  as 

u ■ gKx.  (5.2) 

In  the  presence  of  actuator  and  sensor  dynamics,  (5.2)  becomes 

u - gK-X  (5.3) 


where  n t Rn  Is  the  state  vector  of  the  sensor  system 


and  the  control  variables  are  the  states  of  the  actuator  system. 


e„u  ■ P u + Q U (5.5) 

2 o o 

We  call  x»  the  sensor  variables  and  U,  the  actuator  variables.  In  (5.4), 

(5.5),  £,  and  are  small  positive  scalars  and  matrices  R and  P are 
12  o o 

stable  matrices.  Without  loss  of  generality,  we  assume  that  P and 

o o 

are  invertible.  If  the  feedback  loop  (5.3)  is  opened,  then  for  an 
arbitrary  continuous  function  0(t),  (5.5)  gives 

u ( t ) - - P-1  Q u(t)  + 0(e,)  (5.6) 

o o 2 

for  t - t.,  where  t0  - |max\(P  >1^  2 ^n  and 

X(t)  - - R"1Sox(t)  + 0(e1)  (5.7) 

for  t 2 t where  t^  - |maxA(R^)|  ^e^tnc^.  Thus  the  parameters  and  e0 

determine  the  transient  periods  of  the  open  loop  behaviors  of  the  sensor 

and  actuator.  In  particular,  if  -P  * I and  -R  - I , then  u(t) 

o ^o  m o o n 

and  X(t)  are  0(ei)  and  O(e^)  "close"  to  l)(t)  and  x(t),  respectively. 

We  consider  now  the  behavior  of  the  high  gain  feedback  system 
with  actuator  and  sensor  dynamics,  (5.1),  (5. 3)— (5.5) . First,  we  assume 
that  » g and  >>g,  the  results  are  stated  in  the  following  theorem. 


Theorem  5.1: 

Suppose  and  g are  such  that. 


for  some  positive  number  a 


lim  “ 0, 
e -K)  C1 


lim  - 0 

e,-*0 

r 


(5.8) 


— 


- - ~ 


• w..  — 


119 


and 


- a 


(5.9) 


Let 


'as 


gQ  K a_1P 

o 


and  denote  |j  P KR  *1 
o o o 


a.  B 


b and  | ] ] ] - e , then  if 


(5.10) 


1 + — — a(E,+b)  - 1 

e . c 1 

^ 1 — 
8"  2aabe , ( 1+e , ) = gc 


(5.11) 


then  the  solution  of  system  (5. 1) , (5. 3)- (5 . 5)  can  be  approximated  by 


P t 

r t 

_ i 

‘x(t)‘ 

^*as  e 

'x(o)' 

^■as  e. 

R S 
o o 

■ e 

+ e 

-1  -1 

u(t) 

r _-i_ 

u(0)^ 

1 

-goP  Q KR  S 

O 0 O 0 

x(0) 


-R  S 

o o 

gctP“1Q  KR_1S 
° o o o o 


x(t)  + 0(ei) 


(5.12) 


and 


[A  +gaB  P_1Q  KR_1S  ] 

x(t)  -e°  ooooo  x(Q)  + 0(e^ 


(5.13) 


Proof : Condition  (5.8)  allows  us  to  consider  system  (5.1),  (5. 3)— (5. 5) 
as  a singularly  perturbed  system  in  the  form  of  (C  9),  (C  10)  with  parasitic 
parameter  E. . Thus,  the  block  triangularization  procedure  in  [24,25]  is 
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applicable.  The  remainder  of  the  proof  follows  that  of  Theorem  2.1  which 
is  given  in  the  Appendix  B. 


Corollary  5.1: 

Suppose 


S 

o 


and  let  t 

as 


jraaxX(£as) 


i-l 


e^ne^ 


then 


(5.14) 


X(t)  - x ( t ) + 0(6^  (5.15) 

u(t)  - gaKx(t)  + 0^)  (5.16) 

for  t - t and 

V.S 

[A  +gaB  K]t 

x(t)  - e ° ° x(0)  + 0(GX) . (5.17) 

Thus,  if  ■ £2,  that  is,  a ■ 1 and  the  sensor  and  actuator  systems  are 
such  that  (5.14)  is  satisfied,  then  the  control  u(t)  of  (5.5)  is  O(e^) 
close  to  the  high  gain  feedback  (5.2)  after  an  initial  transient  period  of 
O(e^).  Moreover,  the  motion  x(t)  of  the  high  gain  feedback  system  with 
sensor  and  actuator  dynamics  (5.1) , (5.3)— (5.5)  is  O(e^)  close  to  the  motion 
x(t)  of  (5.1), (5. 2).  If  ct  + 1,  the  approximations  still  hold  but  with  the 
gain  factor  g multiplied  by  a.  In  Theorem  5.1,  an  upper  bound  for  the 
gain  g is  specified  for  which  the  approximations  (5. 12) , (5. 13)  hold.  In 
the  case  when  (5.14)  does  not  hold,  the  next  theorem  examines  the  stability 
of  the  motion  x(t). 


- — — 
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Theorem  5.2: 

Lee  T be  the  transformation  as  in  (2.17)  and  let 

x - [x.  x0]T  = Tx.  Let  P = aP  KR  ''S  B and  suppose  (5.18) 
12'  o 'o  o o o 

ReA(D  < 0 (5.19) 

then  under  the  same  conditions  in  Theorem  5.1  the  motion  x(t)  of  (5.13) 
can  be  approximated  by 

Xj(t)  - eZtx1(0)  + 0(g~l)  + 0(£1)  (5.20) 

and 

x2(t)  - e8rtx2(0)  + Otg'1)  + 0 (ex)  (5.21) 

where  \(Z)  are  the  transmission  zeros  of  the  triple  (KR  ,A  ,B  ). 

o o o o 

Furthermore,  if  g is  sufficiently  large  and  Re\(Z)  < 0,  then 
lim  x (t)  ■ 0 and  lim  x^(t)  * 0. 

fc-*CO  £-+oo  *■ 

Proof : An  O(e^)  approximation  to  x(t)  in  (5.13)  is  the  solution  of  the 
system 

x - (A  + gctB  P_1Q  KR-1S  ) x , x (0)  - x(0).  (5.22) 

a o oo  oooa  * 

Applying  the  transformation  T to  (5.22),  we  obtain 


• 

X - 

al 

(A11  “ A12K2aKla)xaL  + A12K2a  xa^ 

(5.23) 

Ux 

a2 

x + (UH  + Dxa 

I - *2 

(5.24) 
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T 

where  [x  x ] = Tx  , 

al  a2  a 


KR_1S  M_1  - [K.  K.  ] 
o o la  2a 

(5.25) 

and  K„  is  an  m x m nonsingular  matrix  since  by  assumption 
2a 

(5.18)  K2aB°  is 

nonsingular.  Following  Theorem  2.1,  and  letting 

2 - A11  - 

(5.26) 

we  have 

x (t)  = eZtx.(0)  + 0(g  1) 
al  i 

(5.27) 

x (t)  = esrtx9(0)  + 0(g  1) 

(5.28) 

2 


By  Theorem  2.2,  A(Z)  are  the  transmission  zeros  of  the  triple 

(KR~^S  , A , B ).  This  concludes  the  proof, 
o o o o 

Corollary  5.2; 

If  R-1S  = I , then  matrix  Z of  (5.26)  coincides  with  matrix 
o o n 

F, , of  (2.28).  If  in  addition  P_1C)  * I , then  the  0(g)  eigenvalues  of 

11  oo  m 

(5.1) , (5. 3)- (5. 5)  are  the  0(g)  eigenvalues  of  (5.1), (5. 2)  multiplying  by 
the  factor  a. 

Theorem  5.2  shows  that  the  fast  transients  in  high  gain  feedback 
system  are  influenced  by  the  presence  of  sensor  and  actuator  dynamics. 
Neverthelsss , if  (5.19)  is  satisfied,  then  the  fast  transients  decay 
exponentially  in  the  time  scale  gt  to  the  neighborhood  of  Kx  = 0.  The 
presnce  of  sensor  dynamics  alone  alters  the  finite  eigenvalues  of  the  high 
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11  1 ...  i.  , V 


gain  feedback  system  thus  it  affects  the  dynamical  characteristics  of 
x(t)  after  the  fast  transient  diminishes. 

If  there  are  disturbances  and  parameter  variations  in  (5.1),  then 
from  the  above  theorems,  we  see  that  high  gain  feedback  (5.3)  rejects  the 
same  class  of  disturbances  and  parameter  variations  as  when  (5.2)  is  used, 
that  is  when  there  are  no  sensor  and  actuator  dynamics.  Thus,  if  the  gain 
factor  g is  not  "exceedingly  large",  for  example,  if  g ^ g^  as  in  (5.11), 
and  the  sensor  and  actuator  dynamics  do  not  destroy  the  stability  of 
the  system  in  the  sense  of  Theorem  5.2,  then  we  can  proceed  to  design  high 
g3in  feedback  by  assuming  that  the  sensors  and  actuators  do  not  exist. 

We  note  that  if  (5.8)  is  satisfied,  it  is  immaterial  that  and 
I C<,  differ  in  orders  of  magnitudes,  our  conclusion  made  above  is  still 

valid.  The  results  for  case  when  g ■ an<*  8 ” ^2^'’  w^*re  ar>d 

are  finite  positive  constants  will  not  be  presented  in  this  study.  Prelim- 
inary investigations  show  that  the  interaction  of  the  fast  modes  due  to  the 
use  of  high  gain  feedback  3nd  the  natural  fast  modes  in  the  sensor  and 
actuator  gives  rise  to  high  frequency  oscillation. 

Finally,  we  remark  that  the  effect  of  actuator  and  sensor 
dynamics  on  VSS  with  sliding  mode  is  that  nonideal  sliding  mode  occurs, 
that  is,  after  some  time  t^,  the  trajectories  of  the  system  are  confined 
to  some  A-neighborhood  of  the  switching  planes  s(x)  ■ 0 as  in  (4.39). 


Theorem  4.1  shows  that  nonideal  sliding  mode  is  close  to  the  ideal  sliding 
mode,  that  is,  when  there  are  no  sensor  and  actuator  dynamics,  in  the 


While  sliding  inode  is  robust  with  respect  to  sensor  and  actuator 

dynamics,  the  design  of  switching  gains  in  the  variable  structure  control 

has  to  take  into  account  the  matrices  P , Q , R and  S in  (,5.4),  (5.5) 

o o o o 

such  that  the  trajectories  reach  s^x)  -0.  Let  u be  a variable  structure 


control 


o - f°X 


(5.29) 


where  'f  denote  the  switching  gain  matrix  in  (4.19).  Due  to  the  actuator 
dynamics  (5.5),  u is  a continuous  function  and  therefore  sliding  mode  does 
not  occur  on  s(x)-0.  The  trajectories  x(t)  of  the  discontinuous  feedback 
system  (5.1),  (5.4),  (5.5),  (5.29)  are  continuous.  Thus,  the  switching  of 
the  gains  f*"  of  (5 . 2 9 ">  is  regular.  Since  the  feedback  gains  in  '*lt  are  bounded, 
in  the  region  of  the  state  space  of  x where  is  a fixed  gain  matrix,  if 

we  denote,  • K and  let  g-l,  Theorem  5.1  can  be  used  to  show 


u(t)  - 0P  lQ  KR* ' S x(t'  0(C  ) 
O O O O L 


\(t'  - -R  S x 1 1 ) ♦ 0(C) 

O O l 


(5.30) 


(5.3n 


After  some  time  t which  is  of  On,'.  If  (5.14)  is  satisfied,  then 
as  l 


U(,t)  “ O Kx  1 1 ' « 0 (£  ^ ) 


1,5.32' 


and  \ (t)  as  in  (5.15).  If  t is  designed  assuming  that  there  are  no  sensor 
and  actuator  dynamics,  then  if  a > 1 and  |x(t)|  > c^i ^ for  some  positive 
number  c^,  the  trajectories  move  towards  s (x>  - 0.  When  (5.14'  does  not 
hold,  ignoring  the  presence  of  sensor  and  actuator  dynamics  may  lead  to  a 
design  of  n ° which  does  not  guarantee  s(x)  - 0 be  reached  and  hence  nonideal 
sliding  mode  would  not  exist. 


I 


125 


As  we  have  mentioned  in  Section  A. 2 that  in  the  vicinity  of  the 
switching  planes,  the  variable  structure  control  u in  (5.29)  can  be 
approximated  by  U » gKy  as  in  (5.3)  where  Ky  “ 0 approximates  the  switch- 
ing planes  Kx  - 0.  If  Kx  * 0 is  reached  then  the  trajectories  of  the  VSS 
are  in  some  neighborhood  of  Kx  ■ 0,  that  is,  nonideal  sliding  mode  occurs. 

We  note  that  the  reaching  of  Kx  = 0 depends  only  on  a,  the  ratio  between 
and  If  the  nonidealities  in  the  switching  devices  can  be  approxi- 

mated by  (5.3),  then  in  the  vicinity  of  Kx  ■ 0,  nonideal  sliding  mode  is 
the  motion  of  the  high  gain  feedback  system  (5. 1) , (5 . 3)- (5. 5) . The  major 
difference  between  VSS  in  sliding  mode  and  high  gain  feedback  systems  with 
sensor  and  actuator  dynamics  is  their  motions  from  the  initial  states  to  the 
vicinity  of  Kx  •>  0.  In  VSS,  reaching  of  Kx  - 0 guarantees  the  existence 
of  nonideal  sliding  mode  even  in  the  case  when  U * H'°x  has  the  representa- 
tion u - gK\  near  Kx  * 0.  In  contrast,  reaching  of  Kx  = 0 in  high  gain  sys- 
tems depends  on  the  stability  of  the  "fast"  transient  which  is  determined  by  the 
relative  .magnitudes  of  and  g . Thus  we  conclude  that  VSS  in  sliding  mode 

is  more  robust  than  high  gain  feedback  systems  in  the  presence  of  sensor  and 
actuator  dynamics  which  are  neglected  in  the  design  stages. 

5.2  Decentralized  Control 

High  gain  feedback  and  variable  structure  control  can  be  applied 
to  decouple  interconnected  systems.  For  simplicity,  we  shall  limit  our 
consideration  to  systems  with  only  two  interconnected  subsystems.  Let 
S^,  be  the  two  subsystems  given  by 
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Sx:  i1  = A1x1  + B1U1  + E^x2  (5.33) 

•2  22  22  21 

S2:  x - Ax  + B u + E^x  (5.34) 

where  x*,  an  n ^-vector  and  u^,  and  m^-vector  are  the  states  and  control  of 
the  subsystem  S^,  1 * 1,2,  respectively.  Suppose  that  associated  with  S^, 
there  is  a measurement  vector  y^,  an  r^-vector,  of  some  of  the  "local" 
states. 


1 „1  1 
y ■ C x 

(5.35) 

2 .2  2 
y = C x 

(5.36) 

and  the  objective  is  to  regulate  some  "local"  variables,  by  " 

back,  that  is,  the  variables 

local"  feed- 

1 

z = D x , 

(5.37) 

2 _2  2 
z = D x 

(5.38) 

are  to  be  regulated  in  the  subsystem  and  S^,  respectively, 

back  controls  are  of  the  form 

and  the  feed- 

1 ul, 
u = h (y  ), 

(5.39) 

u2  - h2^2) 

(5.40) 

where  h*  denote  some  unspecified  feedback  function.  We  note  that  the 
structure  of  the  interconnected  system  given  by  (5.33),  (5.34)  are  less 
general  than  a linear  multivariable  system  with  an  (n^  + ^^dimensional 


state  space  and  (m^  + m2)  control  variables,  namely,  the  control  vari- 

? 

ables  u do  not  appear  in  and  vice  versa.  The  regulation  problem  we 
pose  for  and  S2  is  a typical  one  in  a decentralized  control  frame- 
work: each  subsystem  has  access  to  measurements  of  its  own  states  and  its 
controller  is  designed  based  on  the  knowledge  of  its  own  dynamics.  We 
recognize  that  if  the  interconnections  E*x^  for  are  considered  as 
unmeasurable  disturbances,  then  using  high  gain  feedback  control  or 
variable  structure  control,  it  is  possible  to  achieve  regulation  of  z~  . 

Our  approach  is  to  construct  "disturbance  rejecting"  observers  for  each 
subsystem  using  the  results  in  Section  3.2  and  4.4,  then  local  high 
gain  observer  feedback  are  designed  accordingly.  For  ease  of  illustra- 
tion, we  first  assume  all  the  subsystem  states  can  be  measured  locally, 

that  is,  » I and  = I and  the  system  matrices  A^,  and  E^  are 

n2  1 J 

known  and  fixed  at  their  nominal  values. 


Lemma  5.1: 


Suppose  local  high  gain  feedback 


u1  = g^x1,  i = 1,2 


(5.41) 


are  applied  to  (5.33),  (5.34)  and  suppose  there  exists  a finite  positive 

i 

constant  a = — and  matrix  K is  designed  such  that  the  0(g)  eigenvalues 
81  1 
of  the  high  gain  system  (5.33),  (5.41)  and  the  transmission  zeros  of  the 

triple  (K*,  A*,  Bi)  are  placed  at  prescribed  locations.  Then  there 

1 2 

exists  a finite  positive  number  t^  such  that  z (t)  and  z (t)  are  reduced  to 
O(g^)  quantities  for  all  t - t^,  that  is,  |zi(t)|  = O(g^),  i » 1,  2 if  and 


i 


) all 


only  if  the  transmission  zeros  of 


IP 

lie  in  the  open  left  half  complex  plane. 


ol 


Proof : By  Theorem  2.1  and  2.2,  the  finite  eigenvalues  of  the  composite 

high  gain  feedback  system  (5.33),  (5.34),  (5.41)  are  the  transmission 

zeros  of  the  composite  system  (5.33),  (5.34)  with  "output"  matrix 

. Since  the  O(g^)  eigenvalues  of  this  system  are  the  union  of 

the  O(g^)  eigenvalues  of  the  high  gain  systems  (5.33),  (5.41)  and 

1 2 

(5.34),  (5.41),  the  regulation  of  z (t)  and  z (t)  results  from  the 
asymptotic  stability  of  the  slow  motion  in  the  composite  system. 

In  order  to  guarantee  that  the  transmission  zeros  of  the  compos- 
ite system  to  lie  in  the  open  left  half  complex  plane,  the  matrices  K^ 

2 

and  K have  to  be  designed  "centrally,"  that  is  for  the  system  (5.33), 
(5.34), 


K1  0 

0 K2 


x„  = A x + B u 
C C c C C 


\l  „l" 

- 1 

-i 

where  A 

_ 

2 2 

and  B * 

B 

0 

c 

• E1  A- 

c 

_0 

B2. 

control  of 

t 

\e  form 

_ 

(5.42) 

, we  design  a high  gain  feedback 


uc  = «l 


x 5 g K x 
c 1 c c 


(5.43) 


K 0 
_0  aK2 

such  that  for  g1  sufficiently  large,  the  finite  eigenvalues  of  (5.42),  (5.43) 
have  negative  real  parts. 

Lemma  5.2 


Let  K={K|K=  block  diag  (K^,^)}  and  let  the  set  of  fixed  modes  of 
(Ac,Bc)  with  respect  to  K be  defined  as 


Then  there  exist 


A(Ac,Bc,K) 


matrices  K and  K 


■H  A (A  + B K). 

K€K  C 

> 

such  that  the  transmission  zeros  of 


t 


f K 0 ^ 

\ 0 ^2  , Ac,  B I lie  in  the  open  left  half  complex  plane  if  and  only 

if  A(Ac,Bc,K)  lie  in  the  open  left  half  complex  plane. 


Proof : The  matrix  g^K£  of  (5.43)  belongs  to  the  set  K.  If  there  exists 
a complex  number  a *=  A(Ac,Bc,K),  then  by  Theorem  2.2  o is  a transmission 
zero  of  (block  diag  (K^,K0),  Ac>Bc)  since  the  O(g^)  eigenvalues  of  (5.42), 
(5.43)  can  be  placed  at  prescribed  locations.  This  proves  the  lemma. 

An  algorithm  to  find  fixed  modes  is  given  in  [49].  By  applying 
the  techniques  in  Section  2.1  for  separating  the  slow  and  fast  modes  of 
high  gain  feedback  systems,  it  can  be  shown  that  the  transmission  zeros 
of  the  composite  systems  are  the  eigenvalues  of  the  matrix 


F A*" 

112  *12 


.2.-1..2 


(5.44) 


where 


M1Ai(M1)-1 


mV  - J1 

J E1 
3 2 


,1*3 


J2  ,mV 


(5.45) 


(5.46) 


EjV«J)'1  ■ 1Ej\l  Ejk2 ^ k ■ !-2-  1 * J 


(5.47) 


K1(M1)_1  - [kJ  K*] 


(5.48) 


for  i.j  * 1,2.  Due  to  the  structure  of  the  matrix  in  (5.44),  algorithm 
such  as  the  one  given  in  [49]  for  finding  decentralized  stabilizing  control 
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can  be  used.  Although  we  manage  to  find  a high  gain  decentralized  control, 
the  design  of  It  requires  the  knowledge  of  the  whole  composite  system. 

In  this  design,  we  have  not  taken  advantage  of  high  gain  feedback;  namely, 
the  insensitivity  property  of  the  slow  motion  with  respect  to  disturbances. 
In  the  following  theorem,  this  property  is  exploited  to  decentralize  the 
design  process. 


Theorem  5.3: 


Under  the  same  assumptions  of  Lemma  5.1,  if  either 


R(E*)E  RCB1) 


(5.49) 


R(E J)c  K(B2) 


(5.50) 


then  there  exists  a finite  positive  number  t^  such  that  |z^(t)|  - O(g^) 
“1 

and  J z — ( t ) J ■ 0(g  ) for  all  t - t. . Moreover,  the  transmission  zeros 

1 1 1 

K 0 

of  „ , A , B are  the  union  of  the  set  of  transmission  zeros  of 

0 K2  C C 
1111'"  2 2 2 

(K  ,A  ,B  ) and  (K  ,A  ,B  ). 

Proof : If  (5.49)  holds,  then  - 0 and  E^0  ■ 0 in  (5.44).  The 

transmission  zeros  of  the  composite  systems  are  the  eigenvalues  of  the 
matrix 


A1  - A1  (kW 
11  12 1 2'  1 

F2  - F2 

Elll  112(V  1 


. 2 o _ l 2 

A",  - Aj2(K‘)  \ 


(5.51) 


If  (5.50)  holds,  E“  ■ 0,  E',,  ■ 0 and  (5.44)  becomes 


AU  - A12CK2>'lKi  EU1  - En2<K2>'lK; 


11  A12^  2;  K1 


(5.52) 


In  both  cases,  the  transmission  zeros  are  the  eigenvalues  of  the 

matrices  on  the  main  diagonal  which  are  the  transmission  zeros  of 

1112  2 2 l9 

(K  ,A  ,B  ) and  (K  ,A“,B  ).  Since  they  are  placed  separately  by  K and  K“ 

to  prescribed  location  on  the  open  left  half  complex  plane,  regulation  of 
1 2 

z (t)  and  z (t)  is  achieved. 

Thus,  if  either  (5.49)  or  (5.50)  holds,  then  the  design  of  the 
high  gain  decentralized  control  (5.43)  is  decomposed  into  two  design 
problems:  each  subsystem  solves  for  its  own  high  gain  feedback  control 

(5.41)  knowing  only  its  own  dynamics.  The  conditions  (5.49),  (5.50) 
mean  that  the  slow  subsystems  of  and  S0  are  connected  in  cascades. 

For  interconnected  systems  with  more  than  two  subsystems,  if  the  slow 
subsystems  of  its  subsystems  have  a multi-level  structure  [50],  that  is, 
the  "A"  matrix  in  the  "slow" composite  system  has  upper  or  lower  triangular 
form,  then  the  design  process  of  high  gain  feedback  (5.43)  is  decentralized. 
This  is  summarized  in  the  next  Theorem. 


Theorem  5.4: 

Suppose  there  are  K subsystems,  Sj , j » 1,...K.  Under  the 
assumptions  of  Lemma  5.1,  if  there  exists  a set  of  integers 
(i^,...,i^  ^)  such  that 


i i 

R(f^j)  c R(B  °) 


(5.53) 


I 


for  k £ {i^ i(j-l} * a " *»••••*  “ 1 with  iQ  = 0,  then  there  exists 

a finite  positive  number  t^  such  that  |z^(t)|  » O(g^),  j - 1,...,k  for 

1 K 

t i t^.  Moreover  the  transmission  zeros  of  (block  diag  (K  , ...,K  ),  , 

1 |C 

block  diag  (B  ))  are  the  union  of  the  sets  of  transmission  zeros  of 

(K^,A^,B^),  j * 1,...,K  where  A^  and  are  defined  analogously  as  in 


(5.42). 


We  now  consider  the  problem  of  designing  high  gain  observers 


for  the  subsystems.  Again  for  simplicity,  we  deal  with  the  system 
(5. 34)- (5. 36) . For  subsystem  S^,  given  the  local  measurements  y*(t),  we 
construct  an  observer  to  estimate  xi(t).  Let  the  systems 


a1  ✓ » i A ..i„i.  *i  . _ i i 

x ■ (A  - g^H  C )x  + B u 


(5.54) 


1 ? 

represent  the  decentralized  high  gain  observers  where  H and  H“  are  the 
design  matrices.  We  assume  that  the  pairs  (A*,C*)  are  detectable. 
Suppose  high  gain  observer  feedback 


u - gjh  8 


(5.55) 


is  applied  to  (5. 33) , (5. 34) . 


Theorem  5.5: 

Suppose  that  matrix  K*  is  designed  such  that  A(KiBi)  and  the 
transmission  zeros  of  (K^.A^.B*)  are  placed  at  prescribed  locations.  The 
matrix  H*  is  designed  such  that  A(C^H^)  and  the  transmission  zeros  of 
(C*,A*,H*)  are  placed  at  prescribed  locations.  Then 
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lim  x*(t)  = 0 and  lim  (x*(t)  - xi(t))  = 0,  i = 1,2 


(5.56) 


t-t-QO 


if  and  only  if  the  set  of  fixed  modes  A(Ac,Bc,K)  defined  in  Lemma  5.2  lie 
in  the  open  left  half  complex  plane. 


Proof : Using  the  results  in  Section  2.1,  it  can  be  shown  that  the  large 
eigenvalues  are  proportional  to  X(KiBi)  and  \(CiHi)  for  sufficiently 
large  gain.  Moreover,  the  transmission  zeros  of  the  triple 
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(5.57) 


are  the  finite  eigenvalues  of  the  decentralized  high  gain  feedback  system 
(5. 33)-(5 . 36) , (5.54),  (5.55)  which  are  the  eigenvalues  of  the  matrix 
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E^r1  - [eJu  £Jk2],  k - 1,2,  1 t j 


Ci(Mi)~1  - [C*  C*] 


(5.60) 

(5.61) 

(5.62) 


for  i « 1,2.  Thus,  the  transmission  zeros  are  the  eigenvalues  of  the 

matrices  - A^2^C2^  ^Cl’  A11  ~ A12^CP  ^C1  an(*  C^e  ^ower  block  diagonal 

matrix  in  (5.58).  The  eigenvalues  of  the  first  two  matrices  are  the  trans- 
111  222 

mission  zeros  of  (C  , A , H X and  (C  , A , H ) , respectively.  The  third 
matrix  is  the  same  as  (5.44)  hence  by  Lemma  5.2,  the  theorem  is  proved. 

The  implication  of  this  theorem  is  that  if  there  exists  decen- 
tralized stabilizing  high  gain  state  feedback,  then  by  substituting  the 
local  states  by  the  states  of  the  decentralized  high  gain  observers  for 
feedback,  the  closed  loop  dynamical  behavior  of  the  interconnected  system 
is  close  to  that  when  full  state  feedback  is  used.  Thus  the  separation 
property  in  linear  observed  feedback  system  holds  in  this  case.  We  note 
that,  in  general,  the  separation  property  does  not  hold  when  regular  gain 

decentralized  observer  feedback  is  applied.  We  observe  that  as  in  full 

1 2 

state  feedback,  the  design  K and  K is  centralized  although  the  designs 

1 2 

of  H and  H are  comp lede ly  independent  of  each  other.  From  the  structure 
of  the  matrix  in  (5.58),  we  see  that  if  either  condition,  (5.49),  (5.50), 
holds  then  again  we  have  a totally  decentralized  design  process.  For 
composite  system  with  more  than  two  subsystems,  the  condition  (5.53)  in 
Theorem  5.4  is  applicable. 
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5 . 3 Design  Examples — High  Gain  Feedback 

The  synthesis  procedures  of  high  gain  feedback  control  developed 
in  Section  3.6  are  now  illustrated  in  the  design  of  control  systems  for  a 
distillation  column  and  a synchronous  machine. 

5.3.1  Distillation  Column 

Precise  control  of  product  quality  in  distillation  columns 
is  considered  to  be  an  important  application  of  automatic  control  in  the 
chemical  and  petroleum  industries.  For  a detail  description  of  the  nature 
of  and  the  difficulties  in  controlling  distillation  processes,  the  readers 
are  referred  to  [51].  Briefly,  distillation  is  used  in  many  chemical  pro- 
cesses for  separating  feed  streams  and  for  purification  of  final  and 
intermediate  product  streams.  In  actuality,  most  columns  handle  multi- 
component  feeds  but  many  can  be  approximated  by  binary  mixtures  [52].  We 
will  consider  only  binary  distillation  columns.  A schematic  diagram  is 
given  in  Figure  5.1.  The  distillation  unit  consists  of  a distillation 
column  which  has,  in  this  case,  eight  plates  (or  trays),  a condenser  and 
a reboiler  unit.  A single  liquid  feed  stream  is  fed  onto  the  feed  plate 
which  is  the  no.  4 plate.  The  vapor  is  condensed  in  a condenser  and  flows 
into  the  reflux  drum.  Reflux  is  pumped  back  to  the  top  plate  which  is  the 
no.  1 plate  of  the  column.  At  the  base  of  the  column,  vapor  boilup  is 
generated  in  a reboiler.  In  [53],  an  approximate  linearized  mathematical 
model  of  this  distillation  unit  for  small  deviations  in  compositions  and 
pressure  about  their  operating  points  is  derived.  The  essential  assumptions 
made  in  this  model  are  that 
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(i)  the  hydraulic  delay  occurring  in  the  liquid  flows  is 

negligible. 

(ii)  the  pressure  variation  is  assumed  constant  throughout 
the  whole  column,  implying  that  the  pressure  drops  across  the  plates  are 
negligible  with  respect  to  the  total  average  change  of  pressure,  and 

(iii)  the  equation  relating  vapor  composition  to  liquid  compo- 
sition on  each  plate  is  linearized  about  the  operation  point. 

The  model  is  described  by  a 11th  order  linear  time-invariant  system  of  the 

• T 

form  x=Ax+Bu+Ef  with  state  vector  x = (x  , x, , . . . ,x 
o o o c l 

x^,  i = 1,...,8  are  the  liquid  composition  on  i*"*1  plate  (mole  fraction  of 
more-volatile  component);  X£  and  x^  are  the  liquid  compositions  in  the  con- 
denser and  reboiler,  respectively  and  p is  the  pressure.  The  control 
variables  are  the  reflux  (u^),  the  steam  input  temperature  to  reboiler 
(Uj)  and  the  cooling  water  temperature  to  condenser  (u^) . These  three 
variables  and  the  feed  input  are  the  disturbance  inputs,  f^fjjf-j  an^ 
respectively.  The  matrices  Aq,  Bq  and  Eq  are  given  in  Table  5.1.  For 
convenience,  we  have  incorporated  the  feed  flow  FF  and  feed  composition 
FCOM  into  a single  disturbance  input  f^,  namely,  f^  = .018FF  + .OllFCOM. 

In  a typical  distillation  control  system,  the  steam  and  cooling 
water  temperatures,  Uj  and  u^,  are  the  only  control  variables.  There  are 
two  proportional  loops,  the  steam  temperature  is  proportional  to  the  tem- 
perature near  the  bottom  of  the  column  and  the  coolant  temperature  is 
proportional  to  the  pressure.  In  [53],  a new  feedback  control  structure 
is  found  by  applying  Rosenbrock's  modal  analysis  procedure  [54].  The  objec- 
tive of  their  design  is  to  minimize  the  steady  state  error  in  output 


I 


139 


product  compositions,  that  is,  xc  and  x^,  caused  by  a constant  disturbance, 

for  example,  the  feed  input.  Thus,  there  are  two  variables  to  be  regulated, 
T 

that  is,  z » [z  z ] and  z = x and  z * v . The  matrix  D where  z = D x 
1 Z lc  z o o o 


F1000000000  0!  a 631 

o " [_0  0 0 0 0 0 0 0 0 1 OJ*  U 

We  note  that  the  first  row  of  D is  orthogonal  to  the  columns  of  B . Hence, 

o o 

T T 

R(D  ; c N (B  ) which  violates  the  necessary  condition  (3.112)  for  the 
o o 

T T 

existence  of  a matrix  K that  satisfies  R(D  ) c R(K  ).  Alternatively  let 

o — 

the  temperatures  of  the  condenser  and  reboiler  Tc  and  T^,  be  the  regulated 
variables.  These  temperature  variables  are  linear  combinations  of  the 


pressure  and  the  respective  composition,  that  is, 

[55  0 0 0 0 0 0 0 0 0 20] 

o = L 0 0 0 0 0 0 0 0 0 55  20J  ' 


(5.64) 


This  matrix  Dq  satisfies  (3.112).  We  verify  that  the  pair  (A^,  A^°IT) 
in  Theorem  3.13K  is  controllable.  Therefore,  there  exists  a high  gain 
feedback  control  u = gKx  such  that  z(t)  is  zeroed. 

From  Table  5.1,  we  see  that  the  first  three  columns  of  B and  E 

oo 

are  identical.  By  Theorem  3.1,  the  sensitivity  of  p(t)  with  respect  to 
reflux,  steam  input  and  cooling  water  temperatures  is  inversely  proportional 
to  the  gain  factor  g.  Since  Tc(t)  and  T^(e)  are  being  zeroed,  the  steady 
state  errors  in  x£(t)  and  ^(t)  (assuming  that  the  feed  input  is  constant) 
are  proportional  to  the  steady  state  error  in  p(t)  caused  by  the  feed 
input.  It  is  of  interest  to  note  that  it  is  concluded  in  [53]  that  steady 
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state  errors  due  to  feed  input  cannot  be  eliminated  which  agrees  with  our 
observation. 

In  [53],  it  is  suggested  that  the  variables  X2>x^  and  p be 

measured.  Since  there  are  four  disturbance  inputs.  Theorem  3.3  demands  at 

least  four  measurements.  If  information  on  x^,  the  composition  on  the  feed 

plate  (no.  4),  is  not  contained  in  the  measurements,  that  is,  the  fifth 

column  of  the  measurement  matrix  C is  a null  vector,  then  the  last  column 

o 

of  Eq  which  is  associated  with  the  feed  input  lies  in  the  null  space  of 
Cq.  Thus  in  order  to  satisfy  the  necessary  condition  R(Eq)  H V (Cq)  = 0 
in  Theorem  3.13H,  we  choose 

00100000000 
00001000000 
00000001000 
00000000001 

T 

that  is,  y * [y^,  y^,  y^,  y^]  and  y^  = y ^ = x^,  v^  = xy»  y^  = P* 

To  check  the  existence  of  matrix  H in  the  high  gain  observer  (3.34),  (3.35) 
which  satisfies  the  specifications  in  Problem  H,  we  proceed  according  to 
Theorem  3.13H.  Since  the  number  of  measurements  and  number  of  disturbance 
inputs  are  the  same,  that  is,  r = p =*  4 the  structural  constraint  in 
Problem  H becomes  R(Eq)  = R(H)  and  by  Corollary  3.9H,  we  need  to  check  the 
transmission  zeros  of  (Cq,  Aq,  Eq).  We  found  that  they  are  the  real  num- 
bers {-.0047,  -.0075,  -.0203,  -.0220,  0.0,  .0176,  .0289).  There  are  three 
transmission  zeros  in  the  closed  right  half  complex  plane.  Hence  the 
resulting  high  gain  observer  which  satisfies  R(Eq)  ■ R(H)  is  unstable.  If 


we  consider  only  the  disturbances  ^2*^3  an<*  ^4*  t*ien  r > P*  According  to 
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Theorem  3.13H,  we  need  Co  check  Che  concrollabilicy  of  (A^ , A^II).  We 
found  ChaC  ic  is  concrollable , hence  Che  eigenvalues  of  Che  slow  sub- 
syscem  of  che  high  gain  observer  can  be  placed  ac  prescribed  locacions. 
However,  chere  is  sCeadv  sCaCe  observacion  error  due  Co  discurbance  f^, 
che  reflux,  (assuming  ChaC  ic  is  conscanc).  If  we  also  measure  x. , chat 
is  r * 5,  y ^ chrough  y^  are  as  before  and  y^  *»  x^,  Chen  ic  is  possible 
co  rejecc  all  disCurbances . We  found  chac  che  pair  (A.^,  A^^)  is 
scabilizable : chere  is  an  unconcrollable  mode  -.0198  in  A^.  Wich  che 
excepcion  of  chis  mode,  che  eigenvalues  of  che  slow  subsysCera  can  be 
placed  and  che  observacion  errors  are  insensicive  Co  all  che  discurbances. 


5.3.2  Synchronous  Machine 

A linearized  mac'nemacical  model  of  a salienc  pole  synchornous 

machine  connecCed  Co  an  infinice  bus  is  given  in  [55].  The  model  is 

described  by  a 9ch  order  linear  Cime-invarianc  syscem  of  Che  form 

x=Ax+Bu  + E f wich  Cwo  concrol  variables  and  Chree  disCurbance 
ooo 

variables.  The  sCace  variables  are  Che  flux  linkages  of  che  field  armacure 

and  damper  windings  (x^,  Xg,  x^,  x^),  che  cerminal  volcage  (x^),  che  field 

volcage  (x^),  Che  speed  of  rocacion  (xg)  and  che  power  angle  (x^).  For 

a chermal  machine,  x^  is  che  mechanical  Corque  T . For  a hydraulic  machine, 

x_  is  defined  Co  be  T + 2T,  where  T.  is  che  gaCe  posicion.  The 

conCrol  variables  are  che  reference  signal  co  che  excicer  (u^)  and  che 

gace  posicion  (u^).  These  Cwo  variables  and  che  discurbance  corque  co  Che 

machine  shafc  are  che  discurbance  inpuCs,  f ^ , f^  and  fg,  respeccively . The 

macrices  A , 3 and  E for  a chermal  machine  are  given  in  Table  5.2.  For  a 
ooo 

hydraulic  machine,  che  A macrix  is  che  same  as  given  in  Che  cable  excepc 

o 

for  che  following  elemencs: 
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“ -26.3,  a^g  “ -.0756,  a^  " -4.87,  a-jjm  —1*24. 


(5.66) 


The  matrices  B and  E are  the  same  as  In  a thermal  machine  except  for 
o o 

the  following  elements: 


b„, y *■  12.42,  b^g 


-110. 


(5.67) 


We  note  that  the  ranges  of  the  input  matrix  for  a thermal  machine  and 
a hydraulic  machine  are  not  the  same.  It  is  expected  that  their  insensi- 
tivity properties  with  respect  to  the  disturbances  will  be  different. 

Since  rank  E ■ 4 1 rank  B - 2,  the  maximum  number  of  regulated  variables 
o o 

is  two.  We  select  the  field  voltage  and  the  speed  of  rotation  to  be  the 
regulated  variables,  that  is,  z^  - x^,  z0  » Xg  and 


(5.68) 


For  a thermal  machine,  the  second  row  of  is  orthogonal  to  the  columns  of 

B and  hence  violates  the  necessarv  condition  (3.112).  If  we  design  a 
o 

high  gain  feedback  u - gKx  ignoring  the  structural  constraint 

R(D  ) S R(K  ),  then  the  steady  state  errors  in  x, (t)  and  xQ(t)  (assuming 

O DO 

constant  disturbances)  are  due  to  f^,  the  shaft  disturbance  torque  since  the 

first  two  columns  of  E coincides  with  B . For  a hydraulic  machine, 

o o 

R(DT)  ^ V(B^)  “ 0.  Since  the  number  of  regulated  variables  is  the  same  as 
o o 

the  number  of  control  variables,  that  is,  S,  - m,  the  structural  constraint 
T T 

becomes  R(D^)  ■ R(K  ).  By  Corollary  3.9K,  we  need  to  check  the  transmis- 
sion zeros  of  (D  , A , B ).  We  found  that  thev  are  the  real  numbers 
o o o 

{-100.,  -60.52,  -34.03,  -26.3,  -19.16,  .0011,  5.02}.  There  are  two 
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transmission  zeros  in  the  right  half  complex  plane.  Hence  the  resulting 

T T 

high  gain  feedback  system  which  satisfies  R(Dq)  " R(K  ) is  unstable.  If 

we  consider  only  the  regulation  of  the  field  voltage,  x, , then  m > £.. 

o 

According  to  Theorem  3.13K,  we  need  to  check  the  controllability  of 

<A11’  ^l?  We  f°un<*  that  it  is  controllable.  Thus,  the  eigenvalues 

of  the  slow  subsystem  of  the  high  gain  feedback  system  can  be  placed 

and  x^(t)  is  regulated.  This  is  also  the  case  for  a thermal  machine. 

On  the  other  hand,  if  we  consider  only  the  regulation  of  the  speed,  Xg, 

then  the  corresponding  pair  (A^'°,  A^'^II)  has  an  uncontrollable  mode 
-2 

1.06  * 10  . Hence  the  resulting  high  gain  feedback  system  which  satisfies 

T T 

R(d‘)  c R(k  ) is  unstable.  If  this  structural  constraint  is  replaced 
by  the  constraint  as  if  x^(t)  is  to  be  regulated,  then  x^(t)  is  zeroed  and 
the  steady  state  error  in  Xg(t)  is  due  to  fg,  the  shaft  disturbance 
torque  only. 

In  [7],  it  is  suggested  that  the  variables  x^,  x^,  Xg  and  x^ 
be  measured,  that  is, 


000010000 

000001000 

000000010 

000000001 


(5.69) 


yl  ■ [yL,  y2>  y3»  y4l  and  y1  - Xg,  y,,  - Xg,  y3  - Xg,  y4  - xg.  For  a 

thermal  machine,  the  second  column  of  E lies  in  the  null  space  of  C . 

o o 

Hence  R(E  ) c M (C  ) which  violates  the  necessary  condition  (3.142).  Thus, 
o o 

there  are  steady  state  observation  errors  due  to  constant  disturbances  f.,. 


the  gate  position.  If  we  consider  only  the  first  and  third  column  of  E^, 


TO 


TO 


then  r > p and  we  found  that  the  corresponding  pair  (A^,  A^^II)  in  Theorem 
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3.13H  Is  stabllizable:  there  Is  an  uncontrollable  mode  -1.24.  With 

the  exception  of  this  mode,  the  eigenvalues  of  the  slow  subsystem  can 

be  placed  and  the  observation  errors  are  insensitive  to  f^  and  f^,  the 

exciter  reference  signal  and  the  shaft  disturbance  torque.  Since  the 

first  and  third  column  of  E in  a hvdraulic  machine  are  identical  to  that 

o 

of  a thermal  machine,  these  conclusions  also  hold.  We  note  that  in  a 
hydraulic  machine,  the  second  column  of  E^  lie  outside  the  null  space  of 
C but  the  span  of  the  second  and  third  column  of  E lie  inside  the  null 
space  of  C^.  If  we  consider  only  the  first  two  columns  of  E^,  we  found 
that  the  pair  ( ) has  an  uncontrollable  mode  5.02.  The  resulting 
high  gain  observer  is  unstable.  We  conclude  that  in  either  a thermal  or 
hydraulic  machine,  the  observation  error  cannot  be  made  insensitive  with 
respect  to  f,,,  the  gate  position  disturbance. 


5 . 4 Design  Example — Variable  Structure  Feedback 

We  illustrate  the  synthesis  procedures  of  variable  structure 
feedback  controllers  in  the  design  of  a model  following  control  system. 

The  idea  of  model  following  is  to  use  a model,  which  specifies  the  design 
objectives,  as  a part  of  the  control  system.  The  objective  of  the  con- 
troller synthesis  is  then  to  minimize  the  errors  between  the  states  of  the 
model  and  the  plant  (the  system  to  be  controlled).  In  the  presence  of 
system  parameter  variations,  the  so-called  "adaptive  model  following  con- 
trol system"  are  developed.  Basically,  the  gains  of  the  controller  in 
these  systems  are  adjusted  so  that  the  plant  states  track  the  model  states. 
Let  the  model  be  given  by 


I 
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A 


where  x^  ^ Rn  is  the  model  state  vector  and  u^  € R1''  is  the  model  input. 


Suppose  the  plant  is  described  by 


x - A x + B u 

p p p p p 


(5.71) 


with  x € Rn  and  u € Rm  are  the  state  vector  and  control  vector  of  the 
P P 

plant,  respectively.  Then  the  error  between  the  model  and  plant  states, 


(5.72) 


is  governed  by 


e - v + (am  - y xP  + Vm  - w 


(5.73) 


In  general,  the  matrices  and  B^  are  known  but  there  are  parameter 

variations  in  A^  and  Bp.  The  objective  is  then  to  design  a controller 

Up  such  that  the  error  e(t)  goes  to  zero  asymptotically.  Our  approach  is 

to  design  a variable  structure  control  u considering  x and  u . as 

P P M 

"disturbances."  Since  Xp  and  u^  are  accessible,  we  let 


u - ¥ x + 4*  e + fii. 
p p p e Mm 


(5.74) 


The  matrices  V , 4*  and  4'  are  defined  as  follows, 
p e M 


4*  - ¥ diag(sgn  x\...,sgn  xn) , xT  - (x\...,xn) 

P P P P P P P 


(5.75) 


4*  ” 4*  diag(sgn  e.,...,sgn  e ),  e - (e. e ) 

e e ± n in 


(5.76) 
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Tm-*m  dlag(sgn  V-“*sgn  “5>*  “m  = (V‘  • ' ,UM)  (5.77) 


where 


«y 


, m.  i i 

••V’  \ 


k for  s.  (e)  > 0 
P i 

-k1  for  s (e)  < 0 
P i 


, i “ 1, 


(5.78) 


^ k^  for  s . (e)  > 0 
e l 

-k1"  for  s (e)  < 0 
e 1 


*S-  <-  t 


^ kjjj  for  s^Ce)  > 0 


^ i 1 j • • • >ni 


y 1 1 ^ ) m 


(5.79) 


(5.80) 


for  s . (e)  < 0 


i-kM  ‘w‘  "i 

In  the  above  equations,  if/*,  k*,  tp*  and  k * are  n * 1 vectors  and  </)*,  kjjj 
are  K x 1 vectors.  The  switching  plane  s^(e)  = 0 is  the  ith  component  of 


s(e)  = Ge  * 0. 


(5.81) 


By  Lemma  4.5,  if 


rank  [ B ! - A ] = rank  B 

Pi  m p P 


(5.82) 


and 


rank  [Bpj B^ ] - rank  B^ 

then  the  VSS  (5.73),  (5.74)  in  sliding  mode  is  described  by 


(5.83) 


-1, 


[I  - Bp(GBp)  G]^e 


(5.84) 


ing  of  the  error  transient  after  sliding  mode  occurs  on  s(e)  ■ 0 is 
achieved. 


We  now  design  the  variable  structure  control  (5.74)  for  an  air- 
craft control  problem.  The  plant  in  this  case  represents  the  three 
degrees-of- freedom  linearized  longitudinal  state  equations  of  a conventional 


[ 


Ki 


subsonic  aircraft,  a Convair  C-131B.  The  model  is  chosen  to  be  the 
estimated  dynamics  of  a large  supersonic  aircraft.  This  is  a benchmark 
problem  considered  in  various  model  following  papers  [56,57,58].  The 
model  and  plant  matrices  are  given  as  follows  [57,58]. 


0.0 

0.0 

-6 

-0.1717 

7.451x10 

-5 

-0.0238 

-7.783x10  J 
-3 

0.0 

3.685x10 

3.0 

1.0 

0.0 

0.0 

1.401x10" 4 

a 

-1.9513 

0.0133 

-2.505xl0"4 

1.0 

-1.3239 

-0.0238 

-0.561 

0.0 

0.358 

-0.0279 

0.0 

0.0 

0.0 

1 

-5.3307 

6. 44 7x10" 3 

-0.2669 

-0.16 

-1. 155xl0~2 

-0.2511 

0.0 

0.106 

0.0862 

0.0 

1.0 

0.0 

0.0 

5. 318x10" 7 

-0.4179 

-0.1202 

2. 319xl0~3 

-4. 6l9xl0-9 

1.0 

-0.7523 

-2 . 387xl0-2 
-2 

(5.85) 

-0.5614 

0.0 

0.3002 

-1.743x10 

(5.86) 


(5.87) 


(5.88) 
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12  3 

The  plant  control  variables  u , u and  u are  the  elevator  command  deflec- 

P P P 

1 2 

tion,  throttle  control  and  flap  command  deflection  where  u^,  are  the 

elevator  and  throttle  inputs,  respectively.  The  plant  state  variables 
12  3 4 

Xp,  x“,  x and  x^  are  the  pitch  angle,  pitch  rate,  angle  of  attack  and 

12  3 4 

air  speed,  respectively.  The  model  state  variables  x^,  x^,  x^  and  x^ 
denote  similar  quantities  as  the  plant  state  variables.  The  plant  matrices 
in  (5.87),  (5.88)  contain  only  one  parameter  a which  varies  between 
-.558  and  -3.558.  We  see  that  conditions  (5.82)  and  (5.83)  are  satisfied. 

By  letting 


(5.89) 


the  equivalent  control  system  which  governs  the  dynamics  of  the  VSS  when 
sliding  mode  occurs  on  the  intersections  of  the  switching  planes  s - Ge  = 0 
is  given  by  a first  order  system 

e^  - -ae^  (5.90) 

where  e^  is  the  first  element  of  the  error  vector  e.  We  conclude  that 
sliding  mode  Is  insensitive  to  the  parameter  O and  by  choosing  a to  be 
positive,  we  guarantee  the  error  goes  to  zero  asymptotically.  Furthermore, 
-a  is  in  fact  the  eigenvalue  of  (5.90),  thus  the  transient  behavior  is 
determined  by  the  selection  of  the  switching  plane  s^(e)  - ae^  + e,,  » 0. 

We  choose  a * 10.  The  hierarchy  of  controls  method  is  used  to  design  the 
gains  in  the  variable  structure  control  law.  We  choose  s0  -*•  s^  -*■  s^ 
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as  the  hierarchy  of  switching  planes  where 

3l(e)  - lOe  + e2  - 0 (5.91) 

s2(e)  - e3  - 0 (5.92) 

s3(e)  - e4  « 0.  (5.93) 

To  conform  with  the  labeling  convention  of  the  hierarchy  of  controls  made 

in  Step  1 of  the  hierarchy  of  controls  method,  we  relabel  sj  » s,,, 
s2  ■ S3  and  S3  “ s^  and  let  U*  be  discontinuous  on  the  switching  plane 
sj!(e)  - 0,  that  is. 


[(kV|e| 


«£T|*pl 


04)T|uMl)8gn  s^. 


1,... ,3. (5.94) 


The  following  conditions  on  the  gains  k^,  k*  and  k^  are  obtained  from  the 
hierarchy  of  controls  method. 
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where  [a^  b^]  > [a9  b9]  means  a^  > a9  and  > b9.  If  these  conditions 
are  satisfied,  then  sliding  mode  will  occur  on  s9(e)  = 0,  then  on  the 
intersection  of  82(e)  = 0 and  s^Ce)  » 0 and  finally  on  the  intersections 
of  all  the  switching  planes.  We  let  the  variable  structure  control 


(5.94)  be 


1 1 
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"3 
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- 
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j 

(5.98) 


We  have  simulated  the  aircraft  dynamics  (5.87),  (5.88)  the 

reference  model  (5.85),  (5.86)  and  the  variable  structure  control  law 

using  the  Runge-Kutta  fourth  order  integration  subroutine  with  automatic 

step  size  adjustments  in  the  IBM  Scientific  Subroutine  Package  [59]. 

Double  precision  arithmetic  is  used  and  the  error  tolerance  for  numerical 
-3 

integration  is  10  . The  parameter  a in  A is  varied  such  that  it  is 

P 

uniformly  distributed  between  -0.558  and  -3.558.  The  reference  model 
1 2 

inputs  Uy(t)  and  u^j(t)  are  taken  to  be  the  same  as  in  [9],  that  is. 


^(t) 


0 < t ^ 20  sec.  and  40  sec.  ^ t ^ 60  sec. 


20  sec.  < t < 40  sec. 


(5.99) 
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The  model  trajectories  which  are  to  be  tracked  by  the  plant  trajectories 
are  illustrated  in  Figure  5.2. 

The  simulation  results  we  have  obtained  indicate  that  sliding 
mode  occurs  as  designed.  Among  these  results,  we  focus  on  the  result 
corresponding  to  the  initial  plant  model  errors  e^(0)  - -1,  e,(0)  - 
e^(0)  » e^(0)  “ 0*  The  error  trajectories  are  illustrated  in  Figure  5.3. 
In  [57,58],  error  responses  for  the  same  initial  errors  using  linear 
model  following  control  and  adaptive  model  following  control  are  obtained. 
A comparison  of  their  responses  to  those  in  Figure  5.3  shows  that 
significant  improvements  on  the  transient  behavior  are  made  using  the 
variable  structure  model  following  control.  We  note  that  sliding  mode 
occurs  on  s, (e)  = 0 and  s^(e)  =*  0 instantaneously  since 
e^(0)  = e^(0)  " 0.  The  trajectory  in  the  phase  plane  (e^,e,)  is  illus- 
trated in  Figure  5.4.  It  is  found  that  sliding  mode  occurs  on  the  inter- 
sections of  all  the  switching  planes  at  t » . 3 sec.  The  time  constant  of 
the  trajectory  e^(t)  for  t 2 .3  sec.  is  approximately  .1  which  confirms 
that  sliding  mode  is  indeed  described  by  (5.40)  with  a - 10. 
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CHAPTER  6 


CONCLUSIONS 

In  this  work,  we  have  developed  to  use  singular  perturbation 

methodology  to  analyze  high  gain  feedback  systems  by  decomposing  the 

motion  into  slow  and  fast  motions.  We  have  shown  that  the  effects  of 

parameter  variations  and  disturbances  on  the  fast  motion  are  attenuated 

by  a factor  of  g ^ where  g is  the  loop  gain.  Their  effects  on  the  slow 

motion  depend  on  the  structure  of  the  parameter  variations  matrices 

($A  and  6B  and  the  disturbance  matrix  Eq.  If  the  ranges  of  6A,  <SB 

and  E lie  in  the  range  of  B , the  input  matrix,  then  these 
o o 

effects  are  also  attenuated  by  a factor  of  g This  new 

methodology  has  not  only  clarified  such  geometric  properties  but  also  made 
possible  a semi-quantitative  analysis  of  the  robustness  properties.  The 
preservation  of  the  stability  under  parameter  variations  is  analyzed 
in  terms  of  preservations  of  the  stability  of  the  slow  and  fast  motions. 
The  stability  of  the  fast  motion  is  affected  by  6B  alone  whereas  the 
stability  of  the  slow  motion  is  affected  by  both  6A  and  5B  and  its  robust- 
ness property  is  the  same  as  the  transmission  zeros.  Using  singular  per- 
turbation techniques,  we  have  obtained  the  lower  bounds  for  the  gain 
factor  that  validate  our  results.  These  bounds  indicate  that  the  gain 
factor  need  not  be  extremely  large.  It  need  only  to  be  large  enough  such 
that  a separation  of  time  scales  takes  place. 

These  analysis  results  are  then  used  to  develop  new  control 
structures.  In  the  case  when  only  some  of  the  state  variables  are 
accessible,  we  have  developed  an  observer  with  high  gain  feedback  loops. 

It  is  the  dual  of  a high  gain  state  feedback  system.  The  high  gain 


feedback  matrix,  in  this  case,  is  the  measurement  matrix  Cq  and  the 
"input"  matrix  H is  a design  parameter.  The  class  of  parameter  varia- 
tions and  disturbances  which  can  be  rejected  are  those  that  lie  in  the 
range  of  H.  This  observation  leads  to  a design  procedure  for  matrix  H 

such  that  the  range  of  H encompasses  the  ranges  of  <$A,  5B  and  E while 

o 

maintaining  the  two  time  scale  property  and  stability  of  the  high  gain 
observer.  For  the  feedback  system  with  two  high  gain  loops:  high  gain 
feedback  from  a high  gain  observer,  we  have  shown  that  the  effects  of 
parameter  variations  and  disturbances  on  the  fast  motion  are  the  same 
as  when  high  gain  state  feedback  is  applied.  The  effects  on  the  slow 
motion  depends  on  the  structures  of  <5A,  <5B  and  Eq.  If  the  ranges  of 
6A,  (SB  and  Eq  lie  in  the  range  of  Bq  and  also  in  the  range  of  H,  then 
the  effects  are  attenuated  by  a factor  of  g Furthermore,  we  have 
shown  that  the  separation  property  of  linear  observer  feedback  systems 
does  not  always  hold.  This  property  holds  only  when  the  observer  loop 
gain  g ig  relatively  larger  than  g and  the  slow  motion  is  0(g  sensi- 
tive to  the  parameter  variations  and  disturbances.  If  only  the  sensiti- 
vity of  the  regulated  variables  are  of  concern,  then  the  structural  con- 
straints on  6A,  6B  and  Eq  are  relaxed  to  a condition  regarding  the 
"controllability"  of  the  slow  motion  with  respect  to  the  disturbance 
inputs.  If  the  motion  of  the  regulated  variables  contains  only  those 
slow  modes  which  are  not  excited  by  the  disturbances,  then  they  are 
0(g  sensitive  to  the  disturbances.  When  the  information  concerning 
the  ranges  of  6A,  6B  and  Eq  is  difficult  to  ascertain,  then  by  assuming 
the  possibility  of  the  worst  case,  this  condition  resolves  to  simply  that 
the  row  space  of  the  regulation  matrix  Dq  lies  in  the  row  space  of  the  high 
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gain  feedback  matrix  K.  This  "worst  case"  condition  is  valid  in  both  cases  when  the 
high  gain  state  feedback  and  high  gain  feedback  from  high  gain  observers 
are  applied.  Moreover,  it  is  the  dual  of  the  insensitivity  condition 
for  high  gain  observer,  that  is,  the  range  of  Eq  lies  in  the  range  of 
H.  Thus,  taking  into  account  the  "worst  case"  structural  constraint, 
the  design  procedure  for  matrix  H developed  earlier  is  applicable  to  the 
design  of  matrix  K. 

A comparative  analysis  of  VSS  and  high  gain  feedback  systems 
has  been  made.  By  showing  that  sliding  mode  in  VSS  is  identical  to  the 
limiting  slow  motion  of  high  gain  state  feedback  systems  as  gain  g tends 
to  inf inity,  we  have  established  that  VSS  and  high  gain  state  feedback  sys- 
tems reject  the  same  class  of  parameter  variations  and  disturbances. 

These  geometric  properties  of  high  gain  feedback  systems  yield  new 
invar iancy  conditions  for  sliding  mode.  By  drawing  an  analogy  between 
VSS  and  high  gain  feedback  systems,  new  control  structures  for  VSS  are 
developed.  We  have  shown  that  Luenburger  observers  remain  effective 
for  variable  structure  feedback.  In  this  case,  the  switching  planes 
are  defined  in  the  observer  state  space  and  the  motion  of  the  resulting 
observer  feedback  system  is  close  to  the  motion  of  VSS  with  switching 
planes  defined  in  the  system  state  space.  In  parallel  to  the  development 
of  high  gain  observers,  we  have  designed  observers  with  variable  struc- 
tures. In  this  case,  the  switching  planes  are  defined  in  the  observation 
error  space  and  we  have  shown  that  variable  structure  observers  reject 
the  same  class  of  disturbances  and  parameter  variations  as  high  gain 
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observers.  Furthermore,  the  insensitivity  property  of  VSS  with  varia- 
ble structure  observers  is  identical  to  the  "two-high-gain-loops"  systems. 
Finally,  we  establish  that  the  sensitivity  of  the  regulated  variables  in 
VSS  is  as  in  high  gain  feedback  systems. 

We  have  demonstrated  that  for  the  purpose  of  sensitivity  reduc- 
tion and  disturbance  rejection,  variable  structure  feedback  is  an 
alternative  structure  to  high  gain  feedback.  Nevertheless  a fundamentally 
important  difference  between  these  two  classes  of  feedback  systems  is  that 
their  motions  are  significantly  different  outside  the  intersections  of 
the  switching  planes  or  equivalently  the  null  space  of  the  high  gain 
feedback  matrix.  In  high  gain  feedback  systems,  this  is  the  fast  motion 
which  is  rapidly  descending  to  this  null  space  or  escapes  to  infinity. 

In  VSS,  the  characteristics  of  this  motion  depend  on  the  design  method 
used  for  the  variable  structure  control.  If  the  hierarchy  of  controls 
method  is  employed,  then  its  trajectory  descends  to  a one  dimensional 
subspace  of  this  null  space  (which  is  a switching  plane),  then  to  the 
intersections  of  two  switching  planes  and  finally  to  the  intersectons 
of  all  the  switching  planes  (the  whole  null  space).  Moreover,  this  motion 
is  not  necessarily  fast.  Our  investigation  on  the  robustness  property  of 
high  gain  feedback  systems  and  VSS  with  respect  to  the  neglection  of 
actuator  and  sensor  dynamics  reveals  that  the  capability  to  reach  this 
null  space  is  susceptible  to  the  neglected  small  time  constants  if  the 
reaching  motion  is  forced  to  be  fast  by  the  use  of  large  feedback  gains. 
However,  since  high  gain  is  not  necessary  in  the  design  of  variable 
structure  feedback,  reaching  can  be  guaranteed  in  VSS.  Hence,  VSS  is 
more  robust  with  respect  to  sensor  and  actuator  dynamics. 


In  conclusion,  we  emphasize  that,  in  this  study,  we  have  not 
considered  the  nature  of  the  disturbance  inputs.  If  we  know  the  type 
of  disturbance  inputs,  then  compensator  feedback  can  be  incorporated 
together  with  high  gain  or  variable  structure  feedback.  For  example, 
in  the  case  of  constant  disturbances,  we  can  apply  high  gain  feedback  to 
reduce  the  steady  state  errors  as  much  as  possible.  If,  owing  to  the 
system  structure,  high  gain  feedback  cannot  make  the  system  motion 
totally  insensitive  to  disturbances,  the  next  step  is  to  incorporate 
integral  feedback  to  remove  the  remaining  steady  state  errors.  This 
amounts  to  a combination  of  fast  and  slow  control  actions  where  the 
sole  purpose  of  the  slow  control  action  (integral  feedback)  is  to  improve 
the  low  frequency  response  and  the  steady  state  characteristics.  We  com- 
ment that  combinations  of  high  gain  or  variable  structure  feedback  and 
compensator  feedback  would  significantly  improve  the  transient  responses 
and  we  foresee  it  to  be  a fruitful  approach  in  control  system  designs. 
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APPENDIX  A 


TRANSFORMATION  OF  HIGH  GAIN  FEEDBACK  TO  SINGULARLY 
PERTURBED  SYSTEM 


Consider  the  transformation  x * Tx  as  two  successive 


transformations,  that  is 


Denote 


x'  - Mx,  x - Tx' . 


then  the  system 


x » MA  M x+MBu 
o o 


by  (A2)  and  (2.9)  is  (2.51),  (2.52).  From  (2.17), 


r1  - 


I 0 

n-m 


C_^C 

2 1 C2 


with  C^,  C9  as  in  (2.12).  The  system 


x - rMA  M-1r'-1x  + TMB  u 
o o 


by  (A2)  and  (A4)  and  x - [x.^  x0]A  is  (2.22),  (2.23)  with 


F ■ A°  - A°  C~^C 
11  11  A12C2  C1 


H2  " (Vl2  + C2A22)C2  • 


T T 

where  M^,  are  (n-m)  x n and  M2,  S.,  are  m * n matrices,  then 


S.M  + S„M„  = I . 
11  2 2 n 


From  (2.12), 


c.  ■ c s,,  c,  = c s. 

1 o 1 2 o 2 

B°  - M0B  , 0 = M.B  . 

2 2 o 1 o 

In  terms  of  (A10),  A°^  in  (A2)  becomes 

A11  - VoV 

Substituting  into  (A6)  yields  (3.38).  Using  the  expressions  (All) 
(A13) , x - M~1r_1x  is 

x - (S.  - S-C~^C..  )x.  + B (C  B )_1x0. 
o 1 2 2 11  000  2 

1 

1 

l 


APPENDIX  C 


APPROXIMATE  SOLUTIONS  OF  SINGULARLY  PERTURBED  SYSTEMS 

In  the  following  lemmas,  it  is  assumed  that  the  matrix  is 
a stable  matrix  and  there  exist  positive  constants  and  such  that 
|w(t)  | * ^ and  | <L( t)  | ^ c2> 


Lemma  Cl: 


The  solution  of  the  singularly  perturbed  system 

£ " 111*  + Ii2a  + niw 
7a  + un2w 


(Cl) 

(C2) 


can  be  approximated  by 


Cut 

5(t)  - e 5(0)  + 


Iu(t-r) 


n u(T)dt  + o(u) 


l -r 

a(t)  - e 22wa(0)  + 0(p). 


(C3) 


(C4) 


Proof : By  application  of  the  variation  of  parameter  formula  to  (C2),  we 

have 


F - 

(>22u 

o( t)  - e o(0)  + 


l It zhI 


e '22  U niu)(t)dT. 


(C5) 


Integrating  the  integral  by  parts,  (C5)  becomes 


- ...  — : 


can  be  approximated  by 


r(  , (^11  ' -12^22^21)t:  . (^11 

5(t)  - e 5(0)  + e 


IlJ22^21)(t"T)nia)(T)dT  + 0(H) 


o”(t)  - e 


Proof:  Let 


^ a»(0)  + .‘-2“  ljjz2l«0)  - r2iZ21«t)  ♦ 0(u) 


(Cll) 


(C12) 


a'  * lulif  + a” 


(C13) 


then  (C7),(C8)  becomes 


^ = (^11  “ ^12^22^21^  + ^12°'  niU 

y“22^21(^ll  ‘ ^12^22^21)5  + (’^22^21“12  + -22)a' 

-1 

+ U(^22^2ini  + n2)u)’ 


(C14) 


(C15) 


Introduce  the  variable 


a - cr'  + yfC 


(C16) 


where 


f ’ C^22  + U^21-12^22)  -21(^11  " ^12^22^21)  + Wl  (C17) 

and  matrix  ft  possesses  a power  series  at  U = 0.  Its  coefficients  can  be 
calculated  from  a recursive  formula  similar  to  (S6)-(BS).  3y  [6,7],  we 
obtain  a equivalent  block  triangular  system 


^ * ^11  ' ^12^22^21  “ '^12'  ^ + ^12°  + IIlt0 


(CIS) 


. “I  2 -1 

yo  ■ f Z 22  + '^22^21^12  + ^ ^12^  + ^^2  + ^22^21^1^<l)' 


(C19) 


This  system  is  now  in  the  same  form  as  (Cl),  (C2) . Hence  (Cll)  and 

l - 

a(t)  - e ya(0)  + 0(u) 


(C20) 


are  obtained  by  Lemma  Cl.  From  (C13)  and  (C16), 


a"  (t)  - a(t)  - [22I2 i?<t)  - yrs(t). 


(C21) 


Substituting  (C21)  into  (C20),  we  obtain  (C12).  This  concludes  the  proof 


Lemma  C3: 


The  solution  of  the  singularly  perturbed  system 


s’  - Inr  + ~ r + n1  o 


(C22) 


r * l21v  + ~ r + n2 


(C23) 


can  be  approximated  by 

, <EirW2222l)t  (^ir^l2^2^21)tv  r_1  , 

5'(t)  - e 5'(0)  - e ^12^22^  (0) 

V r-1  SXii~Zi2^22^21^  St— T>  -1 

+ W22*  (t)  + e [nrI12I22n2]co(T)dT 


-1 


W'  m wX2l5'  + + un 


21l12l22  <• 22 


,0) 


(C28) 


which  is  the  same  form  as  (C14) , (C15).  Hence  by  Lemma  C2,  (C25)  and 

€(t)  „ e(^ir^l2^22^21),:c(0)  + [ e(In-^i2^22^21)(t’T) 


-1 

[n1-Ii2I22!I2la)(T)dT  + 0(y) 


are  obtained.  From  (C26) 


5'(t)  - 5(0  + 


(C29) 


(C30) 


Substituting  (C30)  into  (C29) , we  obtain  (C24).  This  completes  the  proof. 


APPENDIX  D 


TRANSFORMATION  FOR  A "TWO-HIGH-GAIN-LOOPS"  SYSTEMS 


In  order  to  investigate  the  behavior  of  the  high  gain  observer 
feedback  system  (3.57),  (3.58),  it  is  convenient  to  transform  it  into 
singularly  perturbed  systems  of  Appendix  C.  Let  M be  an  2n  * 2n  matrix. 


M. 


Mi 


nxn 


M 


(Dl) 


where  $ is  as  defined  in  (3.38),  M.  and  M.  are  nxn  matrices  such  that 

0-  <5 


M.(B  + 5B)  = 
o o 


(n-m) xm 


(D2) 


M,(B  + 6B)  + M5B  = 0.  (D3) 

o o 

By  the  assumption  that  rank  (Bq  + (SB)  = m for  all  6B,  Mj  exists  and  there 

A 

exists  M^  that  satisfies  (D3) 

Mj  = -MS  (D4) 

where  S is  a solution  of  the  matrix  equation, 


We  define 


SB  = 6B. 


(D5) 


B = B + 6B 
o 


(D6) 


for  convenience.  Introducing  the  transformation 


with  x"  € Rn  m,  x£  € Rm,  ^ Rn  r and  z^  £ Rr.  In  these  new  coordinates, 
system  (3.57),  (3.58)  becomes 


x^'  = A^xJ1  + a-j.2x2  + Elf 


x”  = A21x)’  + A22x^’  + g B2[(K1  + K1)x;['  + (K,  + K2)x2' 


-k^e^  - k2e2J  + E2f 


el  Fllxl  + F12x2  + All£l  + A12e2  + Elf 


(DIO) 


'C2  = F21X1  + F22X2  + ^2iei  + ^22e2  + ^[-(Cj-SC^xJ'  - (C2  - <5C2)x£ 


+ c°e1  + c°z2]  + E2f 


(Dll) 


where 


M ,(A  + 6A)mTL 
0 o 0 


A11  A12 


A21  A22 


m5(ao+5a)m“1  + m.am"1  - maqm  lm5m^1  = 


11 

F12 

21 

F22 

(D12) 


(D13) 
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’ E," 

E, 

m.e  - 

6 o 

1 

.E2. 

. (M-  + M)Ert  - 
6 o 

1 

A 

LE2. 

KMj1  - [Kx  K21 

wTHi^1  - 

[Ki 

k2] 

(D15) 

km'1  - [k2  k2] 

Co*r'*ii's  ’ 

[C\ 

CV 

(D16) 

6cm"1  - 

u 

16C1  fiC2]. 

(D17) 

and  A^j , C°  the  H2  are  aa  defined  in  (3. 38)- (3. 41) . 

We  now  Introduce  another  transformation  that  decomposes  the 
system  (D8)-(D11)  into  slow  and  fast  subsystems.  Let 


T - 


2n-m-r 

K , 


1 0 0 0 
n-m 


0 

0 

0 


0 1 0 
n-r 

[ 0 0 

m 

0 0 1 


where 


*L 


K,  + K, 

A 

-K 

n 

i 

< 

+ 

( 

1 

1 1 

1 

, k2  - 

L i L 

-c'l  + {Cl 

-C2  + 6C2  C° 

In  the  new  coordinates  (x^,  e^,  <j>,  n)t 


(D18) 


(D19) 


(D20) 
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Lemma  C2  shows  Chat 


*(t)  gK_6 . t »(0> 

* - • * + 0(g_1) 


(D26) 


*1  (t>  X^<0>  \ (V^^lKL)(t_T) 


( t ) 


G^O) 


+ 0(g_1). 


Exf (T)dT 


(D27) 


providing  that  Is  a stable  matrix.  From  (D18) 


*1  2 

■ + S 


(D28) 


It  can  be  recognized  from  (D26)  that  [$  n]  are  the  fast 
variables  and  from  (D27)  that  (x'^  e^]  are  the  slow  variables.  Our 

procedure  to  decompose  the  fast  and  slow  variables  follows  closely  that 
in  Section  2.1.  The  variables  x and  e can  be  expressed  in  terms  the 


fast  and  slow  variables  as 


187 


(D29) 


where  f is  a 2n  x (n-m-r)  matrix  whose  columns  span  the  null  space  of 


that  Is, 


K 

5C 


(D30) 


r 


L 
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